B.  G.  Teubners  Sammlung  von  Lehrbüchern 
auf  dem  Gebiete  der  Mathematischen  Wissen- 
schaften mit  Einschlufs  ihrer  Anwendungen. 

Im  Teubnerscheii  Verlage  erscheint  unter  obigem  Titel  in 
zwangloser  Folge  eine  längere  Eeüie  von  zusammenfassenden 
Werken  über  die  wichtigsten  Abschnitte  der  Mathematischen 
Wissenschaften  mit  Einschlufs  ihrer  Anwendungen. 
Die  anerkennende  Beurteilung,  welche  der  Plan,  sowie  die  bis  jetzt 
erschienenen  Aufsätze  der  Encyklopädie  der  Mathematischen  Wissen- 
schaften gefunden  haben,  die  allseitige  Zustimmung,  welche  den  von  der 
'Deutschen  Mathematiker -Vereinigung  veranlafsten  und  herausgegebenen 
eingehenden  Eeferaten  über  einzelne  Abschnitte  der  Mathematik  zu  teil 
geworden  ist,  beweisen,  wie  sehr  gerade  jetzt,  wo  man  die  Eesultate  der 
wissenschaftlichen  Arbeit  eines  Jahrhunderts  zu  überblicken  bemüht  ist, 
sich  das  Bedürfnis  nach  zusammenfassenden  Darstellungen  geltend  macht, 
durch  welche  die  mannigfachen  Einzelforschungen  auf  den  verschiedenen 
Gebieten  mathematischen  Wissens  imter  einheitlichen  Gesichtspunkten 
geordnet  imd  einem  weiteren  Kj*eise  zugänglich  gemacht  werden. 

Die  erwähnten  Aufsätze  der  Encyklopädie  ebenso  wie  die  Referate 
in  den  Jahresberichten  der  Deutschen  Mathematiker -Vereinigung  be- 
absichtigen in  diesem  Sinne  in  knapper,  für  eine  rasche  Orientierung 
bestimmter  Form  den  gegenwärtigen  Inhalt  einer  Disciplin  an  gesicherten 
Resultaten  zu  geben,  wie  auch  durch  sorgfältige  Litteraturangaben  die 
historische  Entwickelung  der  Methoden  darzulegen.  Darüber  hinaus  aber 
mufs  auf  eine  eingehende,  mit  Beweisen  versehene  Darstellung,  wie  sie 
zum  selbständigen,  von  umfangreichen  Quellenstudien  unabhängigen  Ein- 
dringen in  die  Disciplin  erforderlich  ist,  auch  bei  den  breiter  angelegten 
Referaten  der  Deutschen  Mathematiker -Vereinigung,  in  welcher  haupt- 
sächlich das  historische  und  teilweise  auch  das  kritische  Element  zur 
Geltung  kommt,  verzichtet  werden.  Eine  solche  ausführliche  Darlegung, 
die  sich  mehr  in  dem  Charakter  eines  auf  geschichtlichen  und  littera- 
rischen Studien  gegründeten  Lehrbuches  bewegt  und  'neben  den  rein 
wissenschaftlichen  auch  pädagogische  Interessen  berücksichtigt,  erscheint 
aber  bei  der  raschen  Entwickelung  und  dem  Umfang  des  zu  einem 
grofsen  Teil  nur  in  Monographien  niedergelegten  Stoffes  durchaus  wichtig, 
zumal,  im  Vergleiche  z.  B.  mit  Frankreich,  bei  ims  in  Deutschland  die 
mathematische  Litteratur  an  Lehrbüchern  über  spezielle  Gebiete  der 
mathematischen  Forschung  nicht  allzu  reich  ist. 

Die  Verlagsbuchhandlung  B.  G.  Teubner  giebt  sich  der  Hoffnung 
hin,  dafs  sich  recht  zahlreiche  Mathematiker,  Physiker  und  Astronomen, 
Geodäten  und  Techniker,  sowohl  des  In-  als  des  Auslandes,  in  deren 
Forschungsgebieten   derartige   Arbeiten   erwünscht   sind,    zur   Mitarbeiter- 


schaffe  an  dem  Unternehmen  entschliefsen  möchten.  Besonders  nahe  liegt 
die  Beteiligung  den  HeiTen  Mitarbeitern  an  der  Encyklopädie  der  Mathe- 
matischen Wissenschaften.  Die  umfangreichen  litterarischen  und  speziell 
fachlichen  Studien,  welche  für  die  Bearbeitung  von  Abschnitten  der 
Encyklopädie  vorzunehmen  waren,  konnten  in  dem  notwendig  eng  be- 
grenzten Rahmen  nicht  vollständig  niedergelegt  werden.  Hier  aber,  bei 
den  Werken  der  gegenwärtigen  Sammlung,  ist  die  Möglichkeit  gegeben, 
den  Stoff  freier  zu  gestalten  und  die  individuelle  Auffassung  und  Richtung 
des  einzelnen  Bearbeiters  in  höherem  Mafse  zur  Geltung  zu  bringen. 
Doch  ist,  wie  gesjtgt,  jede  Arbeit,  die  sich  dem  Plane  der  Sammlimg 
einfügen  läfst,  im  gleichen  Mafse  willkommen. 

Bisher  haben  die  folgenden  Gelehrten  ihre  geschätzte  Mitwirkung 
zugesagt,  während  erfreulicherweise  stetig  neue  Anerbieten  zur  Mitarbeit 
an  der  Sammlung  einlaufen,  worüber  in  meinen  „Mitteilungen"  fortlaufend 
berichtet  werden  wird  (die  bereits  erschienenen  Bände  sind  mit  zwei  **, 
die  unter  der  Presse  befindlichen  mit  einem  *  bezeichnet):  • 
**P.  Bachmann,  niedere  Zahlentheorie.     (Band  X  der  Sammlung.) 

M.    Bocher,    über    die    reellen    Lösungen    der    gewöhnlichen    linearen 
Differentialgleichungen  zweiter  Ordnung. 

G.  Bohlinaun,  Versicherungsmathematik. 

G.  H.  Bryan,  Lehrbuch  der  Thermodynamik. 

G.  Castelnuovo  und  F.  Enriques,  Theorie  der  algebraischen  Flächen. 
**E.  Czuber,  Wahrscheinlichkeitsrechnung  und  ihre  Anwendung  auf  Fehler- 
ausgleichung, Statistik  und  Lebensversicherung.     (Band  IX.) 
**L.  E.  Dickson,  Linear  Groups  with  an  exposition  of  the  Galois  Field 

theory.       [in  engUscher  Sprache.]       (Band    VI.) 

F.  Dingeldey,  Kegelschnitte  und  Kegelschnittsysteme. 

F.  Dingeldey,  Sammlung  von  Aufgaben  zur  Anwendung  der  Differential- 

und  Inteprralrechnunff. 

G.  Eneström   (in  Verbindung   mit  andern  Gelehrigen),   Handbuch  der 

Geschichte  der  Mathematik. 

F.  Enriques,  Prinzipien  der  Geometrie. 

Pli.  Furtwängler,  die  Mechanik  der  einfachsten  physikalischen  Apparate 

und  Versuchsanordnungen. 
**A.    Gleichen,    Optische    Abbildungslehre    u.    Theorie    der    optischen 

Instrumente.     (Band  VIII.) 
M.  Grübler,  Lehrbuch  der  hydraulischen  Motoren. 
J.  Harkness,  elliptische  Funktionen. 
L.  Henneberg,  Lehrbuch  der  graphischen  Statik. 
K.  Heun,  die  kinetischen  Probleme  der  modernen  Maschinenlehre. 

G.  Jung,  Geometrie  der  Massen. 
G.  Kohn,  rationale  Kurven. 

**A.  Krazer,  Handbuch  der  Lehre  von  den  Thetafunktionen.    (Band  YIJ.) 
H.  Lamb,  Akustik. 
R.  V.  Lilienthal,  DiffHVPntialgeometrie. 


A.  Loewy,  Vorlesungen  über  die  Theorie  der  linearen  Substitutionsgruppen. 
**G.  Loria,  spezielle,  algebraische  und  transcendente  Kurven  der  Ebene. 
Theorie  und  Geschichte.     (Band  V.) 
A.  E.  H.  Love,  Lehrbuch  der  Hydrodynamik. 
A.  E.  H.  Love,  Lehrbuch  der  Elasticität. 
R.  Mehmke,   über  graphisches  •  Rechnen  und  über  Eechenmaschinen, 

sowie  über  numerisches  Rechnen. 
W.  Meyerhof  er,  die  mathematischen  Grundlagen  der  Chemie. 
**E.  Netto,  Lehrbuch  der  Combinatorik.     (Band  VII.) 
W.  F.  Osgood,  allgemeine  Funktionentheorie. 
E.  Ovazza,  aus  dem  Gebiete  der  Mechanik. 
**E.  Pascal,   Determinanten.     Theorie   und  Anwendungen.     (Band  III.) 
S.  Pincherle,  Funktional  -  Gleichungen  und  -Operationen. 
Fr.  Poekels,  Krystalloptik. 
^  A.  Pringsheim,  Vorlesungen  über  Zahlen-  und  Funktionenlehre.    (Ele- 
mentare Theorie  der  unendlichen  Algorithmen  und  der  analytischen 
Funktionen  einer  komplexen  Veränderlichen.)    Bd.  I,     Zahlenlehre. 
Bd.  n.    Funktionenlehre.     (Band  I.) 

C.  Segre,   Vorlesungen   über   algebraische   Geometrie,   mit   besonderer 

Berücksichtigung  der  mehrdimensionalen  Räume. 

D.  Seliwanoff,  Differenzenrechnung. 
M.  Simon,  Elementargeometrie. 

P.  Stäckel,  Lehrbuch  der  allgemeinen  Dynamik. 

P.  Stäckel,  Differentialgeometrie  höherer  Mannigfaltigkeiten. 

O.  Staude,  Flächen  und  Flächensysteme  zweiter  Ordnung. 
**0.  Stolz   und  J.  A.  Gmeiner,   theoretische   Arithmetik.     (Band  IV.) 

R.  Sturm,  Theorie  der  geometrischen  Verwandtschaften. 

R.  Sturm,  die  kubische  Raumkurve. 

H.  E.  Timerding,  Theorie  der  Streckensysteme  und  Schrauben. 

K.  Th.  Vahlen,  Geschichte  des  Fundamentalsatzes  der  Algebra. 

K.  Th..  Vahlen,  Geschichte  des  Sturmschen  Satzes. 

A.  Voss,  Prinzipien  der  rationellen  Mechanik. 

A.  Voss,  Abbildung  und  Abwicklung  der  krummen  Flächen. 

J.  G.  "Wallentin,  Lehrbuch  der  theoretischen  Elektrizität. 
**B.  V.  "Weber,  Vorlesungen  über  das  Pfaffsche  Problem  und  die  Theorie 
der  partiellen  Differentialgleichungen  1.  Ordnung.     (Band  IL) 
*A.  G.  Webster,  the  Dynamics  of  Particles,  of  rigid,  elastic,  and  fluid  Bodies 
being  Lectures  on Mathematical  Physics.  pn  englischer  Sprache.]  (Band  XI.) 

A.  "Wiman,  endliche  Gruppen  linearer  Transformationen. 

W.  Wirtinger,  algebraische  Funktionen  und  ihre  Integrale. 

"W.  Wirtinger,  partielle  Differentialgleichungen. 

H.  G.  Zeuthen,  die  abzählenden  Methoden  der  Geometrie. 

Mitteilungen  über  weitere  Bände  werden  baldigst  folgen. 


Leipzig,  PoststraXse  3.  ^^    ^^   Teublier. 
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Vorwort. 


JL)ie  älteren  Werke  über  Differenzenrechnung  besitzen  nicht  die 
Strenge,  die  unsrer  Zeit  eigentümlich  ist.  In  dieser  Beziehung  ist 
befriedigend  das  Buch  von  A.  A.  Mark  off  (Differenzenrechnung, 
St.  Petersburg  1889 — 1891;  deutsch  von  Th.  Friesendorf  und 
E.  Prümm.  Leipzig,  bei  B.  G.  Teubner.  1896).  Ich  habe  versucht 
einige  Verbesserungen  zu  machen.  Vielleicht  ist  es  mir  gelungen 
einige  Beweise  zu  vereinfachen  und  den  Stoff  systematischer  an- 
zuordnen. Um  die  Arbeit  den  Studierenden  zu  erleichtern,  habe  ich 
speziellere  Probleme  beiseite  gelassen  und  mich  hauptsächlich  auf 
elementare  Fragen  beschränkt. 

Worin  nach  meiner  Meinung  das  Wesentliche  in  der  Differenzen- 
rechnung besteht,  habe  ich  in  einem  Referat  dargestellt.  (Encyklo- 
pädie  der  mathematischen  Wissenschaften,  Bd.  I,  p.  918 — 937.)  Dort 
findet  man  hinreichende  Literaturangaben.  In  kleinem  Räume  konnte 
ich  nur  Sätze  aussprechen,  die  Beweise  sollte  der  Leser  selbst  suchen. 
Um  ihm  die  erforderliche  Mühe  und  Zeit  zu  ersparen,  folge  ich 
dem  Vorschlag  des  Herrn  B.  G.  Teubner  und  übergebe  der  Öffentlich- 
keit dieses  Lehrbuch. 

St.  Petersburg,  den  6.  Dezember  1903. 


D.  Seliwanoff. 
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Erster  Teil. 
Differenzen. 


Erstes  Kapitel. 

Hauptsätze  über  Differenzen. 

1.  Definitionen.     Wenn    eine    Funktion    f{x)    gegeben   ist,   so 
heißt  der  Ausdruck  f{x-{-}i)  —fix)  Differenz  von  f{x),  in  Zeichen: 

f{x  +  }i)-f{x)=^Lf{x). 

Wenn   man    diesen   Ausdruck    in    eine    andre  Form  bringt,    so 
entsteht  die  Gleichung 

f{x  +  h)-f{x)  =  ^{x) 
oder 

/^f{x)^q>{x). 

Hier   kann   x   beliebige  Werte    annehmen,    h  soll  aber  konstant 
bleiben. 

Ist  z.  B,  f{x)  =  x^,  h  =  l,  so  hat  man 

/•(0)  =  0,      /'(1)  =  1,      /■(2)  =  8,        /•(3)  =  27,    /•(4)  =  64,... 

A/'(0)=1,   A/'(2)  =  7,  A/-(3)  =  19,    A/-(4)  =  37, . . . 

Die  Differenz  von  ^f(x)  heißt  zweite  Differenz  von  f(x)  oder 
Differenz  zweiter  Ordnung: 

A'f{x)  =  Af{x-\-h)-Af(x). 

Entsprechend  ist 

A'f(x)  =  AV  (x  +  h)  -  A'f{x),      A^f{x)  ==  A'f(x  +  A)  -  A'f(x), . . . 

In  dem  Beispiel  ist 

A7(0)  =  6,     AY(1)  =  12,     AY(2)  =  18, . . . 

AY(0)  =  6,     AY(1)  =  6,... 

AY(0)  =  0,... 

Seliwanoff,  Differenzenieolinung.  1 


2  Erster  Teil.    Differenzen. 

Diese  Resultate  lassen  sich  in  folgender  Tabelle  zusammenstellen: 


X 

f 

A/- 

AY 

AY 

AY 

0 

0 

1 

6 

6 

0 

1 

1 

7 

12 

6 

2 

8 

19 

18 

3 

27 

37 

4 

64 

2.  Differenzen  einfacher  Funktionen.  Unter  einfachen 
Funktionen  versteht  man  diejenigen,  deren  Differenz  direkt  bestimmt 
werden  kann.  Die  Funktion  heißt  zusammengesetzt,  wenn  ihre 
Differenz  durch  die  Differenzen  andrer  Funktionen  ausgedrückt 
werden  kann. 

In  der  Differentialrechnung  ist 

Die  Differenz  von  x^  hat  keine  einfache  Form,  weil 


h*' 


ist.     Die    einfache   Funktion,    die    der   Potenz  entspricht,    ist  in  der 
Differenzenrechnung  die  folgende 


f(x)  =  x(x  —  'h){x  —  2Ä)  . . .  {x  —  n  —  Ih). 

Die  direkte  Rechnung  ergibt 

f(x  -{-h)  —  f(x)  =  (x  -}-  }i)x{x  —  K)...{x  —  n  —  2h) 

—  x(x  —  h)(x  —  2h) . .  .  (x  ~  n  —  Ih). 

Um    diesen   Ausdruck    zu   transformieren,    setzen   wir    vor   die 
Klammern  den  gemeinsamen  Faktor 

x{x  —  h)(x  —  2h)  . . .  (x  —  n  —  2h) 

und  machen  Reduktionen  in  den  Klammern. 
Auf  diese  Weise  findet  man  die  Formel 

A[x(x  —  h)(x  —  2h)  . .  ,{x  —  n  —  \h)'\ 

=  nhx{x  —  h)(x  —  2h)  . .  .{x  — 

Spezielle  Fälle  dieser  Formel  sind 

Lx  =  h,     L.\x{x  —  h)'\  ==  2hx, 

A[x{x  -h)(x-  2h)]  =  ^hx{x  -  h), 


2h). 


Erstes  Kapitel.    Hauptsätze  über  Differenzen, 

Etwas  einfacher  ist  die  Formel  '•;. 

x{x-h){x-2K)...{x-  n^l  h)  _  ,    x{x-h){x-'lh) . .  .{x-n-  2h) 
^  1.2-3...n  1.2-3...(w-l)  ' 

die  ganz  analog  bewiesen  wird. 

Wir  werden  noch  folgende  einfache  Funktionen  angeben: 

>     w*,     smic,     cosa;. 


a;(a;  +  Ä)(a;  +  2Ä) . . .  (a;  +  w  —  1  i 

Hier  ist  m  eine  Konstante. 

Die  direkte  Rechnung  ergibt 

1 


\{x^h){x^2h)  ...{x-\-n-lh)       (aJ  +  Ä)  (a;-f-  2Ä)  (a;  +  3Ä) . . .  (a;  + wÄ) 

1 

x{x ■\-h){x-\-'ih)  ..  .{x -\-n  —  lh) 
1 


{x-\-h){x-\-2h)  .  .  .  (a;  +  w-lÄ) 
1 


— r-^"  -  -1 

i  K)  La;  +  nh       rc  J 


=  ~  nh 


xix  -\-'h){x  ■\-  ^h) . . .  {x  •{■  nh) 
Spezielle  FäUe  dieser  Formel  sind 

^x^~^'x{x  +  hf      ^x{x-\-h)'^~^^'  x(x-\-h){x-\-2hy 
Etwas  symmetrischer  ist  die  Formel 

^  1.2.3...(w-l)  =_/i.  l-2.S...n 


x{x-\-h){x-\-2h)...{x-{-n-lh)  x(x-Jrh)ix  +  2h)...{x-\-nh) 

Eine  sehr  einfache  Funktion  ist  m*,  weil 

Am^  =  m^+''  —  m^  ==  m^(m''^  —  1). 

Um  die  Differenz  dieser  Funktion  zu  bilden,  genügt  es,  m^  mit 
dem  konstanten  Faktor  m''  —  1  zu  multiplizieren. 

Nach  den  bekamiten  Formeln  der  Trigonometrie  ist 

sin  (x  +  Ji)  —  sin  ic  =       ^  sin  —  cos  (x  +  —j  =  2  sin  —  sin  (a;  +  2"  +  -|) ' 

co8(a;  -\-  h)  —  cos x=  —  2 sin  —  sin  ( ic  +  y)  ==  2 sin  —  cos (^  +  ^  +  ^y 
Also  wird 

A  sin a;  =      2  sin-  cos  (x  -\-~\  =  2  sin  —  sin  (x  -\ — y~j' 

Acos  rr  =  —  2  sin  —  sin  (o;  +  —  j  =  2  sin  —  cos  Ix  -\ ^ — j 

1* 
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Die    wiederholte    Anwendung    der    angegebenen    Formeln    liefert 
uns  folgende  Ausdrücke  für  Differenzen  höherer  Ordnung 


.mx{x  —  h)(x  —  2h)...(x  —  n  —  lh)  ,m  x(x—h)...(x—n—ni—lh)„.. 

1  •  2  •  3  .  .  .  n  1  •  2  •  3  .  .  .  (w  —  wi)  ' 


.n  x{x  —  h){x  —  2Ä) . .  .{x  —  n  —  ih) , « 


1-2  •  3.  ..w 


1 .  2  •  3  .  .  .  (w  -  1)  _  ,      ,  .m        1  •  2  •  3  .  .  .  (n  4-  m  -  1) 


x{x  -^h)  {x  -\-  2h)  . . .  {x  -\-  n  —  Ih)  x(x  ■\-'h)  . .  .  {x  -\- n  -\-  m  —  Ih) 

An       X  X  f     h         i  \w 

m  =  m  [m  —  \)  , 

A  sm a;  =  ( 2  sin  —  j  •  sm  I a;  +  w — ^ — )> 

A  cosa;  ==  (2  sm  —  \  •  cos  I a;  +  w  — ^ — I  • 

Das    sind   die    wichtigsten  Formeln,    die  für  das  Folgende   aus- 
reichen. 

3.  Differenzen,  zusammengesetzter  Funktionen.    Wenn 

wobei  C  eine  Konstante  ist,  so  ist 

f{x  +  Ä)  -  fix)  =  \sp{x  +  h)  +  C\-  {(fix)  +  C]  =  cp{x  +  'h)-  (fix), 

also  Af(x)  =  Ag){x). 

Bei  der  Berechnung  der  Differenz  kann  also  der  konstante  Sum- 
mand weggelassen  werden. 

Wenn  ^r  \       n    r  \ 

f{x)  =  C(p{x) 

ist,  so  ist 

Af(x)  =  C(p(x  +  h)-  C(p(x)  =  C[(p(x  +  h)-  9)(ic)] 

oder  A[C(p{x)']  =  CA(p{x). 

Der  konstante    Faktor    kann    also    vor   das    Zeichen  A    gesetzt 
werden. 

Ist  die  gegebene  Funktion  f(x)  Summe  von  mehreren  Funktionen: 

f{x)  =  f,{x)  +  /•^(a;)  -f- . . .  +  f^(x), 
so  findet  man 

^f{^)=[fi{^-^^)-^fi{^+^)  +  ---+U^+h)]-[f,(x)+f,(x)^-^ 

=[f^(x+h)-f,ix)-]-h[Ux  +  h)-f,{x)]  +  ---  +  \f,(x-\-h)-f^{x)l 

Es  besteht  also  die  Formel 

A[/i(a^)  +  f,ix)  +  . . .  +  Ux)]  =  Af,{x)  +  A/aCo;)  +  • . .  +  A^x). 
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Die   bewiesenen    Sätze   lassen    sich   in    einer  Formel  vereinigen. 

A[0  +  CJ,{x)  4-  C,Ux)  4-  •  •  •  +  G,f,{x)-\ 

wobei  C,  C^,  C^,  .  ■  -,  C/ii  konstant  sind. 

4.  Differenzen  ganzer  Funktionen.    Es  sei 

Es  ist  bequemer,  diese  Funktion  zuerst  in  die  Form 

y/  \        A     ,     A  x  —  a  ,    .  (x—a)(x—a—h)  ,    .  (x—a)(x—a—h)(x—a—2h)   , 
f(x)  =  ^0  +  A-J-  +  A^ Tl +-^3  1.2.3 +  •  •  • 

{x  —  a)(x  —  a  —  h)  ...(x  —  a  —  n—  Ih) 


+  X 


1-2-3 


zu  bringen.  Zur  Bestimmung  der  Koeffizienten  bemerken  wir,  daß 
bei  der  Division  von  f(x)  durch  x  —  a  der  Rest  Äq  ist  und  der 
Quotient 

/•  /'^^        A   j_  A   x-a-h         .  {x-a-h)(x-a-2h) 
l\\^)  —  ^1  +  ^2      1T2        '  '^s  1.2.3  '    ■  ■  ■ 


.    {x  —  a  —  'h)...{x  —  a  —  n  —  l'h) 
^  ^"  1-2.3...W 

Es  ist  noch  zu  bemerken,  daß  f^ix)  die  Form  hat 

Dividiert  man  f^  (x)  durch  x  —  a  —  h,  so  findet  man  den  Rest  Ä^ 
und  den  Quotient 

r-  /  \ A     1      i^j  a;— a  —  2Ä  .    (a;  —  a  —  2ä)  . . .  (x  —  a  —  n  —  lh) 

U[X)  —^21:2""^  ^3    1.2.3     +"-  +  ^^  1.2-3...W  ~* 

Das  erste  Glied  von  f^{x)  ist  p^x"*  —  ^. 

Schließlich  ist 

j-        /  \         1  1  ,     A    X  —  a  —  n  —  Ih  , 

fn-ii^)=Ä„-i  •  iT2'...(r.-l)  +-^     l-2.3...n     =^0^  +  ?• 

Hier  ist  g  eine  Zahl,  die  sich  vollständig  bestimmen  läßt. 

Durch   Division  von   fn—i(x)    durch   x  —  a  —  n  —  1   findet  man 

den  Rest  t— ; — —, rr  und  den  Quotient  ,    -    —     Dieser  Quotient  ist 

l-2...(n  — 1)  ^  l-2...n 

aber  p^.     Es  ist  also 

i-2.3...n  ^  ^o- 

Nachdem   die   Koeffizienten  Äq^  A^,  ä^,  .  .  .,  A„  bestimmt    sind, 
erhält  man  mit  Anwendung  der  Formeln  (Nr.  2): 
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oder 

Es  bestellt  also  der  Satz: 

Die   erste   Differenz   einer   ganzen  Funktion  w*®^  Gerades 
ist  vom  (w  —  1)*®°  Grrade  und  die  fi}^  Differenz  ist  konstant. 

Für  X  =  a  ist 

f{a)  =  A,       A/-(a)  =  /^  •  A,       AY(a)  ^h'-A^,... 

Wir   sehen,    daß    die   Koeffizienten   Aq,  A^,  .  .  .,  An   sich   durch 
Differenzen  ausdrücken  lassen: 

Diese  Werte   der  Differenzen  bestimmen    sich   durch  die  Werte 
der  Funktion  f{x)'. 

f(a),    f{a  +  h),    f(a-^U),...,f{a  +  nh). 

Da  aber  f(x)    eine  beliebige  ganze  Funktion  n^^^  Grades  ist,   so 
besteht  die  Formel 

fr^\  _  fi'r.\  _L  ^-^    ^ /"(<»)    I    {x-a)ix-a-h)     A^fja) 

(x  —  a)ix  —  a  —  h)...{x  —  a  —  n  —  lh)     A"f(a) 


'  1-2-3...W  ^w 

für  jede  Funktion  von  der  Form 

f(^x)=PQX''  -\-  p^X'^-^  +  p^X""-^  ^ ypn-lX+Pn, 

wobei  einige  Koeffizienten  gleich  Null  sein  können. 

Das  ist  die  Newtonsche  Interpolationsformel. 
Es  sei  z.  B. 

a  =  0,    h  =  \,    fix)  =  x^. 
Wir  haben  schon  in  Nr.  1  gesehen,  daß 

m  =  0,     A/'(0)  =  1,     AY(0)  =  6,    AY(0)  =  6. 
Folglich  ist 

x^^x^-  Zx{x  -  1)  +  x{x  -  l){x  -  2). 
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Als  zweites  Beispiel  nehmen  wir 

a  =  0,     /^  =  1,     r(^)  = i.2.3...n 

Hier  ist  y  konstant  vorausgesetzt. 

Nach  der  bewiesenen  Formel  (Nr.  2)  ist 

K^f(c\\  _  y(y-^)---(y-^+^»+^) 

^   /W—  i.2.3...(w-m) 

Daraus  folgt,  daß 

(a!  +  y)(a;  +  y-l)...(a;  +  y-n+l)  ^  y(y-l) . . .  (y-w+1) 
1-  2  ■  3  .  .  .  w  1-  2-  3  ..  .n 

_L  ^    y(y-l)...(y-n  +  2)       a;(a;-l)     y(y- 1)  •  •  •  (y-w  +  3) 

'     l'       l-2-3...(«-l)       "•"       1-2       '  1-2. ..(w-2)         "^ 

a;(a;  — 1) . . .  (a;  — w-fl) 
■  ■  '  "^  l-2-3...'n 

Diese  Relation  besteht  für  beliebige  Werte  von  x  und  y\  sie  ist  von 
Cauchy  benutzt  worden,  um  den  binomischen  Lehrsatz  auf  beliebige 
reelle  Exponenten  auszudehnen.* 

5.  Differenz  des  Produktes  zweier  Punktionen.  Wenn 
man  die  Differenzen  zweier  Funktionen  tf{x)  und  ^(ic)  kennt,  so  ist 
es  leicht,  die  Differenz  ihres  Produktes  zu  bilden. 

Wenn  f{x)  =  ff{x)-'^{x)  ist,  so  ist 

t^f{x)  =  ^{x  +  }i)  •  i){x  +  )i)  -<f{x)' ^{x) 

=  Sjfix)  +  t\^>{x)\y\}{x  +  /i)  -  (fix) •  '^{x) 
=  ^){x)\}\>{x  +  /i)  -  ^{x)\  +  ^{x  +  li)  ■  A(p(x). 
Man  findet  also  die  Formel 

A[(p{x)  ■  ^(a;)]  =  (p{x)-  L-ip{x)  -{-  ^ {x -[■  h)  ■  A(p(x). 

Z.  B.  für  Ä  =  1  ist 

|_  2        J  2  i5 

6.  Differenzen  von  sin  u  und  cos  i*.  Wir  wollen  voraussetzen, 
daß  u  eine  Funktion  von  x  ist,  deren  Differenz  bekannt  ist.  Es  ist 
Asin!*  und  Acos?t  zu  berechnen. 

Wenn  man  x  durch  x  -\-  h  ersetzt,  so  geht  u  in  u  -\-  Au  über. 
Deswegen  ist 

A  sin  w  =  sm  (w  +  Am)  —  sm  tt  =  zsm  -^  •  cos  Im  +  -g- 1> 

/■AN  O     •      ^«         •       /         I     ^"\ 

AcosM  =  co8(m  +  Aw)  —  cosM  =  —  2sm—  •  sm  ( w  +  -g-l- 
*  Cauchy,  Cours  d'analyse,  p.  159.  Paris.  1821, 
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Damit   beschließen   wir    die  Sätze   über   Differenzen   zusammen- 
gesetzter Funktionen. 

7.  Berechnung  äquidistanter  Werte  einer  ganzen  Funk- 
tion.    Es  sei  f{x)  eine  ganze  Funktion  n*®'^  Grades,  deren  Werte 

f{a),    f{a  +  h),    f{a  +  2h),  .  .  .,    f{a  +  nh) 

ausgerecbnet  seien. 

Es  sollen  die  Werte 


/"(a  +  w+l/i),    f{a-\-n  +  2h),    f{a  +  n  ^  ^h),  . . . 
und 

f{a-h),    fia-^h),    f{a-^h),... 

bestimmt  werden. 

Diese  Aufgabe  läßt  sich  leicht  mit  Hilfe  der  Differenzen  auf- 
lösen. 

Durch  wiederholte  Subtraktionen  bestimmt  man  die  Reihe  der 
Werte 


A      f{a),     A      f(a-\-}i),     A  f(a  +  2h),  . .  .,     A  f(a  +  n  -  Ih), 
A2     f{a),     A2     f{a  -f  h),     AY(a  +  2h),  . .  .,     A^f{a  +  n-^2h), 

A'^-'fia),     A'^-'fia  +  h), 
A^     f{a). 

Da  aber  A"f(x)  konstant  ist  (Nr.  4),  so  ist 

AY(«  +  ^)  =  AY(4 
Mit  Benutzung  der  Formel 

A'-^f{x  +  h)  =  A'-^f(x)  +  A'f{x) 

und  durch  wiederholte  Additionen  erhält  man  die  Werte 

A"-V(a  +  2Ä)         =  A"-^f(a  +  h)  -f  A''     f(a  +  h), 

A^'-'fia  +  Sh)         =  A"-V(a  +  2h)  +  A"- V(a  +  2h), 


A2    7(a -I- w  -  1/i)  =  A2     f{a-i-n-2h) +A^     f{a  +  n-2h). 


A      f{a  +  nh)         =A      f{a  +  n-lh)  +A^     f{a  +  n-lh), 
/•(«  +  w  +  Ih)  =  f{a  +  nh)         +A      f{a  +  nh). 

In  derselben  Weise  werden  auch  die  Werte 


f{a  +  n  +  2h),     /"(a  +  w  +  2,h), 
berechnet. 
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Um  fia  —  ^0  ^^  bestimmen,  bemerken  wir,  daß 

A7(a  -  h)  =  A"f(a) 

ist,  und  daß  die  Relation 

A'-'f(x)  =  A'-'f(x  +  70  -  A'f(x) 
besteht. 

Wiederholte  Subtraktionen  geben  uns  die  Werte 

A"- Y(a  -  h)  =  A"- V(«)  -  A"     /  (a  -  h), 
A«- V(a  -h)  =  A"-^f{a)  -  A''"  V(«  -  h), 

A'     f{a  -Ji)  =  A^     f{a)  -  A^     f{a  -  h), 
A      f(a-h)  =  A      f{a)  -  A^     f(a  -  h), 
f(a  -  h)  =  f{a)  -  A      f{a-  h). 

Ebenso  findet  man  auch  die  Werte 

f{a--2h),    f{a-Sh),  ... 
Als  Beispiel  nehmen  wir 

a  =  -l,     h=l,    f(x)  =  x^+x^-2x-l. 
Die  Rechnung  gibt 

A-l)  =  l,   /•(0)  =  -i,   /•(i)  =  -i. 

Durch  Subtraktionen  findet  man 

A /■(-!)  =  -2,     A/-(0)  =  0, 
AY(-  1)  =  2. 
Nach  der  bewiesenen  Formel  (Nr.  4)  ist  aber 

A^f{x)  =  6 
für  jedes  x. 

Die  Anwendung  der  Additionen  gibt  uns 

AY(0)  =  8,       A/-(l)  =  8,       /■(2)=7; 

AV(1)  =  14,     A/(2)  =  22,    /(3)  =  29. 

Um  /"(—  2)  und  f(—  3)  zu  bestimmen,  sind  folgende  Rechnungen 
auszuführen: 

AY(-  2)  =  -  4,       Afi-  2)  =  2,      /•(-  2)  =  -  1 ; 
AY(-3)  =  -10,     A/-(-3)  =  12,    /(_3)  =  -13. 
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Diese  Resultate  lassen  sich  in  der  Tabelle  zusammenstellen: 


X 

f 

A/- 

AY 

AY 

-3 

-  13 

12 

-10 

6 

-2 

—  1 

2 

-4 

6 

-1 

-2 

2 

6 

0 

—  1 

0 

8 

6 

1 

—  1 

8 

14 

2 

22 

3 

29 

Aus  den  Vorzeiclien  von  fix)  erkennt  man,  daß  die  Gleichung 
f[pc)  =  0  drei  reelle  Wurzeln  hat,  die  in  den  Intervallen  (—  2,  —  1), 
(-  1,  0)  und  (1,  2)  liegen. 

Diese  Methode  führt  immer  zur  Absonderung  der  Wurzeln  einer 
numerischen  Grieichung,  sobald  der  absolute  Betrag  der  Differenz 
zweier  Wurzeln  kleiner  als  Eins  ist. 

8.  Relationen  zwischen  den  sukzessiven  Werten  und 
Differenzen  einer  Funktion.  Wir  stellen  uns  die  Aufgabe,  die 
Werte 

f{a),     d.f{a),     AY(«),  ...,     AY(a) 
durch 

f{a),    f{a  +  h),    f{a  +  2h),  .  .  .,     f(a  +  nh) 

auszudrücken. 

Um  die  Schreibweise  abzukürzen,  wollen  wir  folgende  Bezeich- 
nungen einführen: 

f(a)  =  UQ,    f(a-{-h)  =  Uj^,    f(a -{- 2)1)  =  u^,  .  . .,    f(a -\- nh)  =  u„. 

Die  gesuchte  Relation  kann  auf  symbolischem  Wege  ge- 
funden werden. 

In  der  Gleichung 

AUk  =  Uk+i—  U/c 

können    symbolisch   die    Indizes    durch  Exponenten    ersetzt    werden. 
Wenn  man  aus  der  Gleichung 

Am*  =  #+^—  # 

den  gemeinsamen  Paktor  m*  wegläßt,  so  wird 

A  =  w  —  1. 
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Daraus  folgt 

Nach  Multiplikation  mit  w*  und  Ersetzen  der  Exponenten  durch 
Indizes  erhält  man 

Diese  Formel  ist  damit  zwar  gefunden,  aber  noch  nicht  bewiesen. 
Das  kann  in  verschiedener  Weise  geschehen. 
Wenn  man  in  der  Gleichung 

AUi:=Uk  +  i—  Uk 

a  durch  a  -\- h  ersetzt,  so  vermehren  sich  alle  Indizes  um  Eins  und 
es  wird 

Aujc+i=Uk+2—tik+v         X  i  — A' <4 

Durch  Subtraktion  findet  man  a^  , 

A^Uk=Uk+2—  2uk+i+  Uk. 
In  derselben  Weise  fährt  man  weiter  fort 

A^Uk  +  l  =  Uk  +  3—2Uk  +  2+Uk  +  i, 

A^Uk        =  (Uk+S  —  2Uk  +  2  +  Uk  +  l)  —  (tlk  +  2  —  2Mit+l  +  Uk) 

==Uk+B  —  3Wi_|_2  +  ^Uk  +  1  —  Uk, 
A^Mi^-l  =  WA;  +  4—  3Wi  +  3+  3«ti+2—  Uk  +  \, 
A^Uk        =(Mi+4— 3%+3+-^Wi-  +  2  — MA-t-l)  — (W;t  +  3— 3Ma.  +  2+3Ma4.i  — M;t) 

Damit  ist  die  Richtigkeit  der  gefundenen  Formel  für 

m  =  l,  2,  3  und  4 

bewiesen.    Durch    vollständige    Induktion   läßt    sich    zeigen,    daß    die 
Formel  allgemein  gültig  ist. 

Wir   werden    aber    einen    andern    Weg    einschlagen.      Aus    der 
früheren  Betrachtung  folgt,  daß  A'^m*  die  Form  hat 

A'^llk  =  A^llk^m-^r  A^Uk  +  m-l  +  A^UkJf.m-2  -\ \- A^Uk, 

weil  wir   immer   mit   linearen  und  homogenen  Funktionen   der  u  zu 
tun   haben.     Die  Koeffizienten  Aq,  A^,  .  .  .,  Am  sind  von   der  Funk- 
tion fix),    von   a   und   %   unabhängig.     Deswegen    genügt    es,    diese 
Koeffizienten  in  einem  speziellen  Falle  zu  bestimmen. 
Wir  setzen 

a  =  0,     }i  =  \,    f{x)  =  e 
und  finden' 
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Uq  =  1,      U^=C,      W2  =  C^,  .  .  .,      Uk^ra  =  C^^"", 

Am  k  /  ■<  \m 

Uk=C  (c  —  1)  , 

Folglich  ist 

Die  Koeffizienten  A   sind  aber  aus  dem   binomischen  Lehrsätze 
bekannt. 

Es  ist  also  allgemein 

Setzt  man  in  der  bewiesenen  Formel 

Ä;  =  0,     m  =  1,  2,  3,  .  .  .,  w, 
so  erhält  man  die  gewünschten  Formeln 
A  Wo=  Wj  —  Uq, 
A^Wq  =  Mg  —  2u^  +  4*0, 
A^Mo  =  W3  —  3^2  +  3%  —  Uq, 


A  «  ,       W  (W  —  1)  I       /  i  \M 

A   Mo  =  M„  —  W  Un-X  +      \.<^       ^n-2 1"  (—  1)    % 

9.  Jetzt  wollen  wir   zur   umgekehrten  Aufgabe   übergehen.     Es 
soUen  die  Werte 

f(a),    f(a-\-h),    f{a-\-2}i),,..,    f{a  +  nh) 
durch 

f(a),     A/-(a),     AY(a),  ...,     AV(«) 
ausgedrückt  werden. 
Aus  der  Formel 

folgt 

Symbolisch  ist 
oder 

Daraus  ergibt  sich 


Uk+i  =  Uk+  Auk. 

w*+i  =  M^-f  Am* 

M  =  l  +  A. 


M-      =l  +  mA4-^t^A2+---+^". 

1     *     A 

^m+fc  _  ^*  _|.  ^  A  M*  +  ^'^P;^^  A^M*  +  •  •  •  +  A'"m*, 

Mm+i  =  t6i  +  m A M* H — \.2     ^^-t  + ^  ^  ^** 
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Um    diese  Formel    zu    beweisen,   beginnen   wir   mit   den  Fällen 
m  =  l,  2,  3. 

Wir  haben  schon  gesehen,  daß 

Daraus  folgt  (Nr.  3) 

Uk+i  =  {Ui  +  Am*)  +  (A  %  +  A*w*)  =  m*  +  2  A  m*  +  A'wjt, 
diUk+2=  Awi+  2A^Mi+  A^Wi, 
Ma+8  =(uk+2Auk-\-  A^Uk)  +  (Am*  +  2  A^w*  +  A^tt*) 
=  Mi  +  3  A  W;t  +  3  A^  %  +  A^M*. 
Wir   haben   hier  immer  lineare  und  homogene  Funktionen  von 
i(k,  Aiik,  A^Uk,  .  .  .,  folglich  ist 

Uk+m=  B^Uk -^  B^AUk+  B^  A^Uk H h  BmA'^Uk, 

wo  die  Koeffizienten  von  der  Wahl  der  Funktion  f(x)  unabhängig  sind. 
Für 

findet  man 

c*+«=  5^c*  +  B.c'^ic  -  1)  +  B^c^(c  -  1)2 -f  ...  +  5^c*(c  -  1)"» 

odpT 

c^=  B,+  B,{c-1)  -\-  B,{c  ~iy-\-  ■•■  +  Bm{c  -ir. 
Setzt  man  c  =  1  -\-  t  ein,  so  wird 

{l-i-t)"^=Bo+Bj  +  B,t^-{----  +  Bj-. 
Die  gesuchten  Koeffizienten  haben  also  die  Werte 

B,  =  l,    B,=  m,     :B2  =  ^'Tr^'-'-'     ^"^l- 
Es  bestehen  daher  die  Formeln 

^2  =  ^0+  2  A  Mo  +  A^  Mo, 

^3  =  Mo+  3Amo+  dA^UQ-\-  A'mq, 


Wn  =  Wo  +  w  A  tto  +  ^Tg-^  A^Mo  H h  A^Mo- 

Alle  diese  Formeln  können  durch  die  Formel 

m.=  Mo+^AMo+-£7^A«Mo+---+  i.2.3...n       ^  **o 

umfaßt  werden.     Hier  kann  2  nur  die  Werte 

,  0  =  0,  1,  2,  3,  ...,  M 

annehmen. 

Damit  schließen  wir  dieses  Kapitel  und  gehen  zu  Anwendungen 

der  bewiesenen  Sätze  über. 
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Zweites  Kapitel. 
Interpolation. 

10.  G-euaue  Interpolation.  Die  einfachste  Aufgabe  ist  die 
folgende:  Es  sind  die  Werte  einer  ganzen  Funktion  f(t)  von  nicht 
höherem  als  n^^^  Grade  gegeben 

f{a),    f{a  +  li),    fia-\-2h),...,    f(a-\-nh). 

Es  soll  der  Wert  der  Funktion  für  t  =  x  berechnet  werden.     Diese 

Aufgabe  war  schon  gelöst  für  den  Fall,  daß  — t —  gleich  einer  ganzen 
Zahl  ist  (Nr.  7). 

Die  bekannten  Werte  können  wir  wieder  mit 

bezeichnen.     Es  ist  also  der  Wert  der  Funktion  u^  gegeben  für 

2  =  0,  1,2,...,  n. 

Nach  der  bewiesenen  Formel  (Nr.  9)  ist 

,       .         ,    z(z—l).^       ,  ,    z(z—l)...(z  —  n-\-l)j.n 

U,=  t^o+^At^o+-j7Y-A'Mo+---+  i.2.3...n        ^  ^0 

für  diese  Werte  von  0. 

Wenn  wir  setzen 

t—a 

so  ist 

j'w=««.+-^"-A«o+"-f';-;-''A'«o+-- 


j^  (t—a)(t  —  a  —  h)...{t  —  a  —  n  —  lh)  .n 
eine  ganze  Funktion,  bei  der  die  Relation 

F{t)=m 

für  n  -\-l  verschiedene  Werte  von  t  besteht.  Diese  Relation  muß 
identisch  erfüllt  sein,  weil  die  algebraische  Grleichung  Fit)  —  f{t)  =  0 
von ,  nicht  höherem  als  w*®"  Grade  nicht  mehr  als  n  verschiedene 
Wurzeln  besitzen  kann.     Wenn  wir  t  =  x  setzen,  so  finden  wir 

f(x)  =  f{a)  +  ^  A  fia)  +  ^"^"^^.^^'r^  AY(«)  +  •  •  • 
{x  —  a){x  —  a  —  h) . . .  {x  —  a  —  n  —  Ih)  j.n  j^f  \ 

Das  ist  die  Newtonsche  Interpolationsformel,  der  wir  schon 
früher  begegneten  (Nr.  4). 
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11.  Andre  Aufgaben  der  Interpolation  wollen  wir  nur  für  ganze 
Funktionen  von  nicht  höherem  als  dem  4*®**  Grade  behandeln. 

Es  sind  Werte  einer  ganzen  Funktion  f{t)  vom  4*®''  Grade  ge- 
geben 

f{a),    fia  +  h),    f(a  +  2h),     f{a+'dh),     f{a-h4.h). 

Man  soll  für  t  =  x  den  Wert  der  Funktion  und  ihrer  Ableitungen 
berechnen. 

Nach  der  New  ton  sehen  Formel  ist 

m = txa)  +  '-^  ■  Af  («) + '■'-f'-^r"^  AYw 

.   it-a){t-a-h)it-a-2h)  A3r/^\  ■  {t-a)it-a-h){t-a-2h){t-a-dh)  *  4/r«\ 
+  1.2- 3Ä»  ^  /  W+  1.2.3.4Ä*  ^  ^W 

für  jeden  Wert  von  t. 

Der  Einfachheit  wegen  setzen  wir  t  =  a  +  112,  die  Funktion  f(t) 
verwandelt  sich  dabei  in  F(z)  und  wir  finden 

F(,)  =  f(a)  +  {Lfia)  +  i^  AY(a)  +  '-^^^^^  AY(«) 

oder 

F(^)  =  f{a)  +  ^ A/-(a)  +  (0^-.)  ^  +  (^«-  3^^+  2^)-^ 

Durch  DifiBrentiation  erhält  man  die  Formeln 

F^   (^)  =  A/-(a)  +  (2.-l).^  +  (3^^-6^  +  2).^ 
+  (4^3- 18^2+ 220 -6). 3^^, 

JP"  (^)  =  AY(a)  +  (^  -  1)  AY(a)  +  (6^^-180  +  11) -^T^' 

i^'"  (0)  =  AY(«)  +  (2^  -  3)  ^^, 

^""(0)  =  AY(a). 

Wenn  z  =  y  für  ^  =  a;  ist,  so  sind  die  gesuchten  Werte 

/'(^)  =  |i^'(2/),     rix)  =  l,F'\y),     f"'(x)=^_lF"'(y), 

Hier  ist  y  durch   — v--  zu  ersetzen. 

Es  seien  Zv  B.  für  eine  ganze  Funktion  f(t)  von  nicht  höherem 
als  S*««*  Grade  die  Werte 
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m  =  l,    /■(2)  =  3,    /•(4)  =  53,    /■(6)  =  199 
gegeben. 

Es  soll  /""(!)  bereclmet  werden. 
In  diesem  Falle  ist 

a  =  0,     h  =  2,     x  =  l,     y  =  Y 

r(l)  =  i[AY(0)-|AV(0)] 
AY(0)  =  53-2-3+  1=48, 
A3 /-(O)  =  199 -3- 53 +  3-3-1  =  48, 

r(l)  =  1(48-24)  =  6. 
12.  Für   eine   ganze  Funktion  4*®*^  Grades   und   für  /i  =  1  seien 

f{0),    A/-(0),    AV(0),     AY(0),    AY(0) 
gegeben.     Es  sollen  die  Werte 

f{^'    /(fo)'     /©'■■'    fiU 
bereclmet  werden. 

Wenn  wir  die  Bezeicbnungen 

öf{a)  =  f(a  +  ^)  -  f{a),     d'f{a)  =  d/-(a  +  ^)  -  (J/^Ca),  .  . . 

einführen,  so  ist  nach  der  Newton  sehen  Formel 

m  (w  -  1)  (m  -  2)  (m  -  3)  ^^    .q. 

Die  Werte  d/-(0),  dY(0),  dY(0)  und  dY(0)  lassen  sich  aus  der 
Formel 

/•(«)  =  AO)  +  <A/-(0)  +  '-fc±>  AY(0)  +  '('-'H»-^)^s^(o) 

bestimmen,  die  für  jedes  t  gültig  ist.     Man  findet 

fQö)  =  m  +  ~  Ar(0)  -  4  AY(0)  +  ^  AY(0)  -  i^  AY(0), 
f{^)  =  f(0)  +  Ä  ^  «0)  -  i^  ^V(O)  +  j^  AY(0)  -  ^  AY(0), 
f  (l^)  =  rtO)  +  A  ArtO)  -  ^  AY(0)  + 1^„  AY(0)  -  S  AY(0), 
^(Ä)  =  «0)  +  ^  A  AO)  -  ^  AY(0)  +  ^  AY(0)  -  ^  AY(0). 
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Die  Reihe  von  Subtraktionen  gibt 

ärn  =  /o  A/-(0)  -  4  AY(0)  +  ^  AY(0)  -  i?^  AY(0), 

*^Q  =  i^  ^/"W  -  2^  ^V(0)  +  ^  AY(0)  -  S  ^'«0), 

*^(rö)  =  ,ö  ^^(0)  -  4  ^ V(0)  +  ^  AY(0)  -  ^  AY(0), 
*^Q  =  fo  ^«0)  -  4  ^V(0)  +  s^  AY(0)  -  ^  AY(0). 

*  V(0)  =  i^o  ^V  (0)  -  j4  a  Y(0)  +  ^  AY(0), 
*V(i^)  =  iJö  ^V(0)  -  j4  aY(0)  +  j?55_  AY(o), 

*V(|,)  =  i^  A Y(0)  -  j^  AY(0)  +  ^  AY^O). 
«V(0)  =  ^AV(0)-Ji,AY(0), 


*'/Q=i4^v(o)-2iSiö^v(o).  ^ 

Die  gesuchten  Formeki  sind  also 

d/(0)  =  i  A^(0)  -  4  AY(0)  +  ^  AV(0)  -  S  AY(0), 

-»^(0)  =  J^  AY(0)  -  -^4  AY(0)  +  j^  AY(0), 

*V(0)  =  j^AY(0)-^AY(0), 

*V(0)  =  iö^AV(0). 

Wir  haben  schon  gesehen  (Nr.  7),  daß  die  Gleichung 
(p(t)  =  t^+t^-2t-  1  =  0 
im  Intervall  (—  2,  —  1)  eine  Wurzel  besitzt.     Die  Gleichung 

f{x)  =  (p{-2-\-x)  =  0 
hat  aber  eine  Wurzel  zwischen  0  und  1. 

Seliwanoff,  DifTerenzenrechnang.  9 
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Nach  der  angegebenen  Methode  können  die  Werte 

berechnet  werden. 

Wir  wollen  voraussetzen,  daß  /"fe)  und  /"(^^r^)  verschiedene 
Vorzeichen  haben.  Dann  hat  die  Gleichung  (p(f)  =  0  eine  Wurzel 
im  Intervall  (-  2  +  ^^  -  2  +  ^^)- 

Jetzt  ist  in  cpif)  die  Substitution 

^^-"^  +  "10" 
anzuwenden.     Wir  erhalten  eine  neue  Funktion 

die  auch  eine  Wurzel  zwischen  0  und  1  besitzt. 

Nach  demselben  Verfahren  findet  man  zwei  Zahlen  %  und  ^-i-, 
/c  \  /e  4- 1\  ^^ 

für    die    /i(^)  imd  /i(    lo    )  "^^rschiedene    Vorzeichen   haben.      Die 

Wurzel  der  Gleichung  g)(t)  =  0  liegt  also  im  Intervall 

{-^  +  ro+-w'    -^  +  ro  +  '4^y 

Auf  diese  Weise  kann  man  die  Wurzel  der  Gleichung  mit  beliebiger 
Genauigkeit  berechnen. 

13.  Wir  wollen  jetzt  folgende  Aufgabe  der  Interpolation  lösen. 
Es  ist  eine  ganze  Fimktion  4*®^  Grades  zu  bilden,  die  für 

M,  f{o),  /■(«  +  h),  f{a  +  h),  f{a  +  2h) 

vorgeschriebene  Werte  besitzt. 
Die  Funktion 

Fix)  =  f{a)  +  ^  A/-(a)  +  (^-«K^-^«-^)  AV(a) 
hat  die  Eigenschaften 

F{a)  =  fip),  F(a  +  h)=  f{a  +  h),  F(a  +  2h)  =  f{a  +  2h). 
Dieselben  Eigenschaften  bestehen  auch  für 

0{x)  =  Fix)  +  {x  —  a){x  —  a  —  h){x  —  a  —  2h)(p{x), 

wie  auch  die  Funktion  (pix)  gewählt  wird. 
Aus  den  Forderungen 

^\a)  =  fip),     0'{a  -f-  Ä)  ==  f'(a  +  h) 
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lassen  sich  die  Werte  (p{a)  und  (p{a  -{-  h)   berechnen.     Deswegen  ist 
die  Funktion 

vollständig  bestimmt. 

Es  soUen  z.  B.  die  Bedingungen 

/(0)  =  1,    r(0)  =  l,    /•(!)  =  2,    /•'(!)  =  3,    /•(2)  =  7 

befriedigt  werden.     In  diesem  Falle  ist 

A/-(0)  =  1,     A/'(l)  =  5,     AY(0)  =  4, 

F{x)  =  l  +  x  +  '^^^^-4c  =  2x^-x-\-l, 

^{x)  =  2x^-x^-l-\-  x(x  -l){x-  2)(p(x) 

=  2x^-x+l-\-{x^-  3a;2  +  2x)(p(x), 

0'{x)  =  4:x-l  +  {x^-  3x^  +  2x)(p'(x)  +  (3a;2  -Qx  +  2)(p{x), 

^'(0)  =  -  1  +  2g)(0)  =  1,     9'(0)  =  1, 

^'(1)  =  3 -9.(1)  =  3,  9^(1)  =  0, 

q)(x)  =  1  -{•  X '  {—  1)  =  —  X  +  1, 

fix)  =  2a;2-  a;  +  1  +  {x^ -  3a;«+  2x)  (-  a;  +  1) 

=  -ic*+4a;3-  3a;2+a;+  1. 

Wir  wollen  noch   ein  Beispiel  behandeln.     Es   soll   eine   ganze 
Funktion  3*®°  Grades  die  Bedingungen  erfüllen 

/•(o)=^i,  /•(i)  =  i,  Ai)  =  i,  r(i)  =  4. 

Nach  dem  Taylor  sehen  Satze  hat  die  Funktion 

F{x)  =  l  +  {x-l)-\-  4^^^  =  2x^-  ^x  +  2 

die  Eigenschaften 

F(l)  =  l,     F'(l)  =  l,     i^"(l)  =  4. 

Dasselbe  gilt  auch  für 

^{x)  =  2x^-^x-\-2-\-{x-  Ifcpix) 

bei  willkürlicher  Wahl  von  (p{x).  Da  aber  nur  eine  Bedingung  zu 
erfüllen  ist,  so  können  wir  cpipc)  konstant  voraussetzen.  Es  soll 
^(0)  =  1  sein,  folglich  ist 

2-qp(0)  =  l,    9)(0)=-l 
"^^  f(x)  =  2x^-  dx  +  2  +  {x-  iy=x^-x'+  1. 

Dieselbe  Aufgabe  könnte  man  in  andrer  Weise  lösen. 
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Den  Bedingungen  f{0)  =  1,  /"(l)  =  1  genügt  die  Funktion 

^(ic)  =  1  +  x(x  —  l)gp(a;)  =  1  +  (ic^  —  x)(p(x). 
Die  Ableitungen  dieser  Funktion 

0'  (x)  =  {x'-x)(p'  (x)  +  {2x  -  l)(pix), 
0"{x)  =  (x^-  x)(p"(x)  +  2(2x  -  \)^>\x)  +  2if{x) 
haben  für  x=\  die  Werte 

^'(1)  =  9.(1)  =  1,     ^"(1)  =  29'(1)  +  29.(1)  =  4. 
Daraus  folgt 

9,(1)  =  1,     ^''(l)-!, 

Kf{x)  =  \  ■\-  {x—V)  =  X, 

fix)  =  1  +  (a;^  —  (£)x  =  x^—  x^-\-  1. 

Ganz  ähnlich  kann  man  ganze  Funktionen  höheren  Grrades  bilden, 
die  vorgeschriebene  Werte 

f{a),  f{a),  r{a),...,f^<'-^\a), 

f{a  +  h),      f'(a-\-h),        f"(a  +  h),...,f^-^-')(a-\-h), 

f(a  +  mh),  f'(a  +  mh),     f"{a  +  mh),  ...,  /"("». -i)(a  +  mh) 
besitzen.     Der  Grad  dieser  Funktion  kann  nicht  die  Zahl 

a  +  tti+ttg-^ h«w— 1 

übersteigen.     (Vgl.  Markoff,  Differenzenrechnung,  S.  1 — 6.) 

In  den  betrachteten  Fällen  der  Interpolation  wurde  vorausgesetzt, 
daß  die  Werte  der  unabhängigen  Variabein  in  äquidistanten  Intervallen 
auftreten,  weil  nur  in  diesem  Falle  die  Differenzen  anzuwenden  sind. 

14.  Angenäherte  Interpolation.  Für  Werte  von  t,  die  in 
einem  Intervall  liegen,  wollen  wir  die  gegebene  Funktion  f{t)  an- 
genähert durch  die  ganze  Funktion  F{{)  ausdrücken,  die  gewissen 
Forderungen  genügt.  Der  Fehler  dieser  Annäherung  soU  abgeschätzt 
werden. 

Z.  B.  soU  die  ganze  Funktion  F(t)  vom  nicht  höheren  als 
M*®"^  Grade  sein  und  die  Gleichung 

m  =  ^(0 

erfüllen  für 

t  =  a,     a  +  hj     a  +  2h,  ...,  a -\- nh. 

Es  ist  uns  schon  bekannt,  wie  man  eine  solche  ganze  Funktion 
bildet  (Nr.  10). 
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Es  sei  X  eine  Zahl,  die  im  gegebenen  Intervalle  liegt.  Wir 
werden  näheriingsweise 

fix)  =  Fix) 

setzen  und  es  ist  der  Fehler  dieser  Annäherung  abzuschätzen. 
Die  Funktion 

i^(4\        wrA        W   (t-  a)  {t-  a-h) .  .  .  {t -  a-nh) 
fit)  -  Fit)  -  K i.2.3...(n  +  l) 

bei  beliebigem  Werte  von  K  verschwindet  für 

t  =  a,     a-\-h,     a-\-2h,  ...,  a -]- nh. 

Wir  wollen  K  so  wählen,  daß  diese  Funktion  noch  für  t  =  x 
verschwindet.  Hier  ist  x  eine  Zahl,  die  von  a,  a-\-h,  a -\-  2h,  ..., 
a  -\-  nh  verschieden  ist. 

Die  Gleichung 

^lO  ==  f{t)  -  Fit)  -  K i.2.3...(n  +  l) " 

hat  dann  mindestens  n  -\-  2  Wurzeln.  Nach  dem  Rolleschen  Satz 
hat  die  (w  +  1)*®  Ableitung  von  ^(i)  mindestens  eine  Wurzel,  die 
wir  mit  |  bezeichnen  woUen,  Diese  Zahl  ist  eine  mittlere  zwischen 
der  größten  und  der  kleinsten  unter  den  Zahlen  a,  a  -\-  nh  und  x. 
Es  ist  aber 

K=fin+x)il), 

Es  besteht  also  die  Formel 

fix)  =  f{a)  +  y  Af(a)  +  <"-f  .<^-r ''  AY(a)  +  ■  - 

{x—a){x—a—li). .  .{x—a—n—lh),n_p/  ^  ,  (x-a){x-a-h). .  .(x-a-nh)  /.„4.i)/fc\ 
^  1.2.S...n.h-  ^  '^""^^  l-2-3...(n  +  l)  ^  ^^^^ 

Das  ist  die  Newtonsche  Interpolationsformel  mit  dem  Rest- 
gliede  von  A.  Cauchy  (vgl.  Cauchy,  Oeuvres,  1  serie,  tome  5, 
p.  422  oder  Comptes  rendus,  tome  11  [1840],  p.  787). 

Besonders  sind  die  Fälle  zu  bemerken 

n  =  0,     fix)  =  fia)  +  {x-a)f'i^\ 

»^  ==  1,     fix)  =  fia)  +  ^  Afia)  +  (^-«)(^-«-^)  ^».(|), 

n  =  2,     fix)  =  fia)  +  ^  Afia)  +  ^""1^.^^;,"" '^  AY(a) 
(x-a)(x-a-h){x-a-2h)  .,,,,^. 

1     *     A    *    O 

Die  erste  dieser  Formeln  ist  die  Lagrangesche. 
Durch  Weglassung   der  Restglieder  entstehen  Näherungsformeln, 
deren  Fehler  durch  die  Restfflieder  ausgedrückt  sind. 


folglich  ist 
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15.  Wir  wollen  einen  besonderen  Fall  der  Interpolation  betrachten. 
Für  t  =  X  soll  die  gegebene  Funktion  f(J)  angenähert  durch  eiue 
solche  ganze  Funktion  3*^"^  Grades  F{t)  dargestellt  werden,  daß  die 
Bedingungen 

-F(- 1)  = /"(- 1),    F{Q)^m,    F'(0)  =  f{0),    i^(l)  =  /•(!) 
erfüllt  seien. 

Die  Funktion  F{t)  hat  die  Form 

m  -  /•(- 1)  +  ^  A/^c- 1)  +  ^^*  AY(- 1)  +  ^'\']l^'-'^  Ä, 

wo  die  Konstante  Ä  sich  aus  der  Bedingung 

•F'(0)  =  /-'(0) 
bestimmen  läßt.     Wir  finden 

F\t)  =  A/-(-l)  +  ^(2t  +  1)AY(-  1)  +  |(3^^-  1)^, 

/■'(0)  =  Ar(-l)  +  4AY(-l)-{^, 

A  =  6A/-(-  1)  +  3AY(-  1)  -  6f'(0), 

A  =  6[/-(0)  -  /•(-  1)]  +  31/(1)  -  2f(0)  +  /•(-  1)]  -  6f'(0). 

Nach  Reduktionen  ergibt  sich 

JL  =  3/-(l)-3/'(-l)-6/"(0). 

Um  die  Interpolationsformel  mit  dem  Restgliede  zu  finden, 
betrachten  wir  die  Funktion 

wo  die  Konstante  K  so  gewählt  ist,  daß  0(x)  gleich  NuU  ist. 
Wenn  wir  die  vier  Zahlen 

-  1,     0,     1,     X 

der  Größe  nach  ordnen,  so  entstehen  solche  vier  Zahlen 

%,      «2?      %;      ^i> 
die  den  Bedingungen  genügen 

»1  <C  «2  *^  ^3  "^  a^. 

Diese  Zahlen  genügen  der  Gleichung  ^(^)==0.  Nach  dem 
Roll  eschen  Satze  hat  die  Gleichung  ^'(^)  ==  0  mindestens  drei  solche 
Wurzeln  \,  h^,  \,  daß 

a^<i\<  a^<\K  a^<})^<  a^. 
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Keine  der  Zahlen  i^,  \  mid  Jg  ist  also  gleich  Null.  Außerdem 
ist  diese  Gleichung  noch  für  ^  ==  0  befriedigt.  Die  Gleichung  ^ '  (^)  =  0 
hat  also  mindestens  vier  verschiedene  Wurzeln.  Nach  dem  Roll  eschen 
Satze  hat  die  dritte  Ableitung  von  ^\t)  oder  ^^^{f)  mindestens  eine 
Wurzel  t  =  ^,  die  zwischen  a^  und  a^  liegt.     Da  aber 

iet,soist  *-(0  =  ^-©-if 

K=r{%). 

Wir  haben  also  die  folgende  Formel  bewiesen 
fix)  =  f(-  1)  +  ^  Af(-  1)  +  ^^^ AY(-  1)  +  ^  (^+^i)^^(^-i) 

I        {X-\-l)  X      {X  —  1)  ^jy  /gv 

"T"         1 .  2  •  3  •  4       '       ^^^' 

Der  Koeffizient  Ä  hat  den  angegebenen  Wert. 

16.  In  andern  Fällen  werden  die  Restglieder  mit  Hilfe  der- 
selben Betrachtungen  gebildet. 

Wenn  z.  B.  die  Näherungsformel 

genau  ist  lur 

t  =  a,    a  +  h,    a  -\-  2h,    a  +  3ä, 

f'(t)  =  F'(t)     für     t  =  a,     a  +  h,    a -\-  2h, 

f"(t)  =  F"(t)    für     t  =  a, 

dann  ist  F(f)  eine  ganze  Funktion  vom  V^^  Grade  und  die  Näherungs- 
formel f(x)  =  F(x)  enthält  den  Fehler 

(x  -ay(x-a-  Ä)*  (x-a-2hy(x-a-Sh) 


1-2-3-4-5-6-7 


■rm- 


Die  Zahl  |  ist  zwischen  der  größten  und  der  kleinsten  unter 
den  Zahlen  a,  a  -{•  dh,  x  enthalten.  (Vgl.  Markoff,  Differenzen- 
rechnung, S.  6 — 8.) 

17.  Anwendung  der  angenäherten  Interpolation  auf  die 
Berechnung    der    Wurzeln    numerischer    G-leichungen.     Wir 

machen  die  Voraussetzung,  daß  die  numerische  Gleichung  fiji)  =  0 
im  Intervall  (a,  h)  nur  eine  Wurzel  besitzt  und  daß  /"  (t)  immer 
dasselbe  Vorzeichen  in  diesem  Intervall  behält. 

Da  die  Vorzeichen  von   f(a)  und  /"(&)  verschieden  sind,   so  ist 

fia)-m<Ö. 

Wir  werden  f(f)  näherungsweise  durch  eine  solche  lineare 
Funktion  J^(^)  ausdrücken,  für  die 

F{a)  =  f{a)  und  F{b)  =  f{b). 
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Hier  ist  h  =  h  —  a  und  deswegen  ist 

F 
Für  t  =  X  wird 


F(t)  =  f{a)^l^Mib)-fia)l 


Den  Wert  x  wählen  wir  so,   daß  F{x)  =  0  wird.     Folglich  ist 

(.-«)(.- 6)  =  '^-;>;r%,f<o. 

Es  liegt  also  a;   zwischen  a   und  &.     Für   diesen  ic   kann   nicht 
/""(!)  =  0  sein.     Die  Gleichung 

^{t)  =  fit)-Fit)  =  0 

wäre  dann  befriedigt  für  t  =  a,  x  und  &.  Nach  dem  Rolleschen 
Satze  hätte  die  Gleichung  ^  (t)  =  0  mindestens  zwei  Wurzeln  zwischen 
a  und  h,  und  die  Funktion 

^"(t)  =  f"{t) 

müßte  in  diesem  Intervall  ihr  Vorzeichen  wechseln,  was  gegen  die 
Voraussetzung  ist. 

Weil  F{x)  =  0  ist,  so  ist 

Der   betrachtete   Wert   von   x   läßt    sich    einfach    durch   Deter- 
minanten ausdrücken.     Es  ist 

m  m  .  fia)  f{b) 

a        b       '      1        1 
Wenn  t  im  Intervall  (a,  h)  liegt  und 
fia)f"{t)>0 
ist,  so  ist  /•(«)/""(!)  >0  und  folglich 

fia)nx)<0, 

da  {x  —  a){x  —  h)  negativ  ist.  Die  Wurzel  der  Gleichung  f(f)  =  0 
liegt  also  zwischen  a  und  x. 

Wenn   wir   a   und   h   als   erste   Annäherungen   der   Wurzel   be- 
trachten, so  sind 

n   -n        h     .bf(a)-am 

»1  -  »;    Ol  -  f(^a)  -  m 

zweite  Annäherungen. 

Es  sei 

mr(t)<o,  oder  mr{t)>o, 

für  Werte  von  t,  die  im  Intervall  (a,  &)  liegen. 
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Dann  ist 

fib)f"(l)>0     und     f(h)f(x)<0. 

Die  Wurzel  von  fit)  =  0  liegt  zwischen  h  und  x.  Die  zweiten 
Annäherungen  sind  in  diesem  Falle 

In  derselben  Weise  kann  man  von  a^  und  h^  zu  den  folgenden 
Annäherungen  «2  und  \  übergehen  usw.  (Vgl.  Weber,  Lehrbuch 
der  Algebra,  I,  §  116.) 

18.  Anwendung  der  angenäherten  Interpolation  auf  die 
Berechnung    der    Logarithmen    und    Antilogarithmen.     Wir 

nehmen  z.  B.  eine  Tabelle  fünfstelliger  gewöhnlicher  Logarithmen 
aller  ganzen  Zahlen  von  1000  bis  einschließlich  9999.  Mit  Hilfe 
der  Werte  log  iV  und  log  (N  +  1)  soll  log  (N  +  x)  berechnet  werden, 
wobei  0  <  ic  <  1  ist. 

Die  Literpolationsformel 

f{x)  =  m  +  X  A/-(0)  +  -^^  /•"(!) 

gibt  die  Näherungsformel 

log  {N-\-x)  =  lo^N-\-x  [log  {N+x)-  log  N^ 

mit  dem  Fehler 

x{l  —x)       log  e 


Es  ist 


Daraus  folgt 


1-2        {N+IY 

e  <  3  <  -/lÖ,     log  e  <  l, 

iV+|>103. 


£<j6^e     oder     8<l-\0-\ 


Dieser  Fehler  hat  keinen  Einfluß  auf  die  fünfstelligen  Logarithmen. 

19.    Die  Berechnung  der  Antilogarithmen  führt  auf  die  Frage: 
Es  sind  die  Werte  gegeben: 

logJ\^==a,     log(iV^+  1)  =  &. 

Aus  der  Gleichung 

log(iV+ ?/)  =  a  +  x 

soll  y  berechnet  werden  für  einen  Wert  von  x,  der  zwischen  0  und 
b  —  a  liegt. 
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Wir  setzen  ^,  .       ^^  _^    , 

f(x)^N  +  y=  10«  +  ^. 

Die  Anwendung  der  Interpolationsformel 

fix) = m  +  ^  A  m  +  ^(^^^  fi^i) 

gibt  die  Näherung 

y  = 


b  —  a 
mit  dem  negativen  Fehler 

Es  ist 

Aus  der  Relation 

e  1 10  ==  10 


110-,-^- 
log  e 


folgt 

Da  aber  10«  +  ^<10*  ist,  so  ist 

-^<i.l0^.(loge)-2.(^)'. 
Wenn  wir  bemerken,  daß 

10*  =  J^^+  1,    &  -  a  =  log  (i\r  +  1)  _  log  i^=  log  (l  +  ^), 
so  finden  wir 

-^<^(J!f+l).(loge)-^[log(l+l)J. 

Dieser  Ausdruck  kann  noch  vereinfacht  werden. 
Für  positive  z  ist 

9  {2}  =  log  (1  +  ^)  —  <s  log  e  <  0, 
weil  gD(0)  =  0  und 

sind.     Deswegen  ist 

log  (l  +  1=)  <  ^  log  e 
und  folglich  ,^  ,  ^ 

Da  aber  JV+  1  ^  10*,  N^  10^,  so  ist  die  positive  Zahl 
-£<2.10-3. 

In  diesem  Kapitel  ist  die  Interpolation  nur  als  Anwendung  der 
Differenzenrechnung  behandelt  worden.  Ausführliche  Mitteilungen 
über  diesen  Gegenstand  findet  man  im  Referat  von  J.  Bauschinger 
(Encykl.  der  Math.Wiss.,  I,  S.  799—820). 

*  1  10  ist  der  natürliclie  Logarithmus  von  10. 
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Drittes  Kapitel. 
Angenäherte  Berechnung  bestimmter  Integrale. 

b 

20.  Einleitende  Bemerkungen.     Um  das  Integral    1  f{x)dx, 

a 

WO    a  <.h    ist,    angenähert   zu    berechnen,    gebraucht   man   eine    der 
Interpolationsformeln  (Nr.  14 — 16) 

fix)  ^  F(x)  +  (pix)f(pK^). 

Hier  sind  F(x)  und  (p(x)  rationale  ganze  Funktionen  und  |  ein 
Mittelwert. 

Nach  Integration  findet  man  die  Näherungsformel 

b  b 

jf{x)dx=  JF{x)dx 

a  a 

mit  dem  Fehler 


=-f(p{x)f^p\^)di 


Es  ist  I  eine  unbekannte  Funktion  von  X]  wir  wissen  nur,  daß 
für  a<x  <b  auch  a<^<b  ist.  Der  Fehler  der  Näherungsformel 
ist  nur  dann  abzuschätzen,  wenn  die  Funktion  q)(x)  für  alle  x  zwischen 
a  und  b  ihr  Vorzeichen  behält.  Nach  einem  Satze  der  Integral- 
rechnung ist  dann 


=  f^p){r})J(p(x)dx, 


wo  a  <Crj  -Cb  ist. 

Die    Aufgabe    ist    also    auf    die    Integration    ganzer    rationaler 
Funktionen    zurückgeführt.      Die    Rechnungen   sind    noch    einfacher, 
wenn  die  Grenzen  des  bestimmten  Integrals  —  1  und  -f  1  sind,  z.  B. 
1 


/' 


2  2 


3-^2    '     5. 
—  1 

Um  die  Grenzen   des    gegebenen  Integrals   iu  —  1    und   +  1    zu 
verwandeln,  setzen  wir 

X=-l\){t) 

und  es  sollen  die  Bedingungen  erfüllt  werden 
t^(-l)  =  a,     T/;(l)  =  &. 
Nach  der  Newtonschen  Interpolationsformel  ist 

a;  =  a  + -X_(&  _  a)  =  — -<  + -— . 
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Vermöge  dieser  Substitution  findet  man 

6  1 

a  —1 

In  den  folgenden  Paragraphen  werden  wir  verschiedene  Inter- 
polationsformehi  anwenden.  Daraus  entstehen  verschiedene  Methoden 
zur  Berechnung  des  gegebenen  Integrals. 

21.  Methode  der  Rechtecke.     In  der  Formel 
b  1 


wo 


a  —1 


setzen  wir  (Nr.  14) 

Daraus  entsteht  die  Näherungsformel 


/• 


1 


mit  dem  Fehler 
1 


f{t  +  l)J"(t^)dt  =  ^'{t)  f(t  +  1)^^  -  2j:'(r). 
—1  —1 

Hier  ist  -  1<  t  <  1. 
Um  die  Funktion  f{x)  einzuführen ,  bemerken  wir,  daß 


0 

ff 


2 
Wir  finden  die  Näherungsformel 

f(x)dx  =  (b  —  a)f(a) 
mit  dem  Fehler 

wo  a  <  I  <  &  ist. 

Geometrisch  ist  das  gegebene  Integral  eine  Fläche,  die  in  der 
Näherungsformel  durch  ein  Rechteck  ersetzt  wird.  Daher  rührt  der 
Name  der  Methode. 

Um  den  Fehler  zu  vermindern,  zerlegen  wir  das  Intervall  {a,  &) 
in  n  gleiche  Intervalle,     Das  Integral  verwandelt  sich  in  die  Summe 

b  a-\ a+2 6 

/»                           /^         n               f*           n                             /» 
f{x)dx  =  /  f{x)dx  +  /  f{x)dx  +  •  •  •  4-  /  f{x)dx. 
U  tJ  b  —  a  *y,  -b  —  a 

a  a  a-\ a-\-n  —  1 
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Nachdem  wir  auf  jeden  Summanden  die  Näherungsformel  angewandt 
haben,  finden  wir 


/ 


f{x)dx  =  ^—^{%+yt  +  2/2  +  •  •  •  +  Vn-l), 


wo  der  Kürze  wegen  die  Bezeiclinung 

eingeführt  ist. 

Der  Fehler  dieser  Näherungsformel  ist 


s  = 


1     (&-a)\., 


[/'(y  +  /"(y +  •••  +  /"(!«)], 


2 

wo  a  +  (k  —  1)  - — -  <  li  <  ö^  +  ^  - — -  ist. 

Diesen  Ausdruck  kann  man  vereinfachen. 

Es  wird  vorausgesetzt,  daß  f'{x)  endlich  und  stetig  im  Intervall 
(a,  h)  ist.  Diese  Funktion  hat  dann  eine  obere  Grenze  A  und  eine 
untere  Grenze  B.* 

Aus  der  Reihe  von  Ungleichheiten 

folgt 

nÄ  ^  /•'(!,)  +  /"(y  +  •  •  •  -f  f'iL)  ^  nB 

"^^  ^  >>  r(gi)+/^(y+---+r(i^o  >  jg 

Nach  einem  Satze  von  Weierstraß  kann  f'{x)  jeden  Wert 
zwischen  A  und  B  annehmen.     Deswegen  gibt  es  eine  solche  Zahl  |, 

Es  ist  also 

^-1^'- /■'«)■ 

Wenn  wir  die  Interpolationsformel 

benutzen  und   dieselben  Rechnungen  wiederholen,    so  finden  wir  die 
Näherungsformel 


ß 


f{x)dx  =  ^(^1  +  2/2  +  •  •  •  +  ^«) 


mit  dem  Fehler  ,  ,,       ,, 


*  Tannery,    Introduction    ä   la    thöorie    des    fonctions    d'une    variable, 
p.  101  et  126. 
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Wenn  f{x)  immer  dasselbe  Vorzeichen  im  Intervall  {a,  h)  hat, 

6 

SO  ist  der  Wert  von    /  f{x)  dx  zwischen 

a 

-^  (^0  +  2/1  +  2/2  +  •  •  •  +  «/n- l) 
und 

enthalten,  weil  e  und  £j  verschiedene  Vorzeichen  haben. 


22.     Um  das  Integral 

6  1 

ff{x)dx  =  ^ß{t)dt 


1 
zu  berechnen,  kann  man  nicht  die  Formel 

gebrauchen,  weil  t,  der  Faktor  bei  J^'(^i),  das  Vorzeichen  im  Intervall 
(—  1,  1)  wechselt  und  der  Fehler  der  Näherungsformel  nicht  ab- 
geschätzt werden  kann.     Es  ist  zu  setzen  („  Markoff  %  S.  53) 

und  man  findet 

1 

—1 

Daraus  folgt  die  Näherungsformel 

b 

fm = (j  -  «)/(^) 

a 

mit  dem  Fehler 

Ä(*-«)r(i)=i?-^'g|'- 

Um  eine  genauere  Formel  zu  bekommen,  zerlegt  man  das  Intervall 
(a,  h)  in  n  gleiche  Intervalle.  Wenn  man  dieselbe  Methode  wie  in 
Nr.  21  benutzt,  so  findet  man  die  Näherungsformel 

b 


ß 


f(x) dx  ==  -^  (y^+  ^^+  •  •  •  +  ysn-i) 
22  3 


mit  dem  Fehler  ^      ,,       .,     „„.., 

JL  .  (&-«)'   r  (g) 

2w**  1-2-3  *   1-2 ■ 

Hier  ist  wieder 

2/,  =  /"(a+i^~)- 


Drittes  Kapitel.    Angenäherte  Berechnung  bestimmter  Integrale.  31 

23.  Methode  der  Trapeze.     Zur  Bereclinung  des  Integrals 

b  1 

a  —1 

werden  wir  die  Funktion  J^(t)  angenähert  durch,  eine  lineare  Funktion 
^(i^)  ausdrücken,  die  den  Bedingungen  genügt 

0i-  1)  =  /(-  1),     ^(1)  =  J^(l). 

Dadurch  entsteht  die  Interpolationsformel 

Ht)  =  J^i-  1)  +  *4i  [i^(l)  -  ß(-  1)]  +  (i+Ml^)  j:''(^J 

Die  Anwendung  derselben  gibt 
1 

J/(0^^  =  ^(-l)+/(l)-|^"(r). 
—  1 

Daraus  folgt  die  Näherungsformel 

6 

fmdx  =  ^-^[fia)  +  f{b)-] 

a 

mit  dem  Fehler 

1-2-3  '   1-2* 

Die  rechte  Seite  der  Näherungsformel  ist  der  Flächeninhalt 
eines  Trapezes. 

Nach  Zerlegung  des  Integrals  in  n  Summanden  findet  man  die 
Näherungsformel 

6 

//•(^)^^  =  ^  [(2/0  +  2/i)  4- (2/x  +  2/2)  +  •  •  •  +  (2/„  _  1  +  2/J], 


oder 


ff{x)dx  =  ^  (i/„  +  2i/,  +  2y,  +  •  •  •  +  22/„_,  +  y„) 


mit  dem  Fehler 

1     (6-a)»    /•"(§) 


h  =  - 


n*     1-2-2      1-2 


Wenn  f'(x)  dasselbe  Vorzeichen  hat  für  alle  x  zwischen  a  und 
b,  so  haben  s^  und  s  (in  Nr.  22)  verschiedene  Vorzeichen.     Deswegen 


ist  das  Integral    /  f(x)dx  zwischen 
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-^  (^1  +  2/1  +  2/1  +  •  •  •  +  y^n-i) 

.  2  2  2  2 

und  ,  _ 

enthalten. 

24.    Simpsonsche    Formel.      Wenn    wir    die    Interpolations- 
formel (Nr.  15) 

/(O  =  /(-  1)  +  4-'  ^  ^^-  1)  +  ^TTY-*^'-^(-  1) 

anwenden,  so  wird 

J/(Oa!^  =  2/(-  1)  +  2A^(-  1)  +  I- A2/(-  1)  -  ^  j:iv(r). 
—1 

Es  ist  aber 
A/(-  1)  =  J^(O)  -  /(-  1),     AV(-  1)  =  /(l)  -  2/(0)  +  ^(-  1); 

Nacli  einigen  Reduktionen  finden  wir  die  Simpsonsche  Formel 

b 

}(x)dx = ^  [f(a) + if{^) + m] 

mit  dem  Fehler  ,,       ,.  /.ttt/j-x 


ff 


1-2-3-4-5     1-2-3-4 

Wenn  f(x)  eine  ganze  Funktion  vom  3*^"^  oder  niedrigeren  Grade 
ist,  so  ist  das  Restglied  gleich  Null  und  die  Simpsonsche  Formel 
ist  ganz  genau.     Z.  B.: 


1 


ß 


(2  a;  +  3)  (2ic  -  1)  (2ic  +  5)  £?a;  =  I  [-  3  •  5  +  5  •  7J  =  ^ 

ö 
Wenn   aber  f(x)   eine    andre   Funktion    ist,    so    kann    man    den 

Fehler  verkleinern,   nachdem  man  das   gegebene  Integral  in  n  Sum- 
manden zerlegt  hat.     Man  findet  dann  die  Näherungsformel 

b 


f(so)dx  =  -^(yQ-{-4yi  +  2^/1 -f  42/3  H 1- 4^2n-i+ 2/») 


a 

mit  dem  Fehler 


1       (p-ay        /-i^d) 

£  = 


n*     1-2-3-4-5     1-2-3-4 

In  diesem  Kapitel  haben  wir  mit  geringen  Änderungen  die 
Darstellung  vom  Herrn  A.  A.  Mark  off  angegeben.  (Differenzen- 
rechnung, S.  50 — 60.) 

Wir  schließen  damit  die  Anwendungen  der  Sätze  über  Differenzen 
und  gehen  zu  den  Summen  über. 


Zweiter  Teil. 
S  n  m  m  e  n. 


Erstes  Kapitel. 

Unbestimmte  und  bestimmte  Summen. 

25.  Definitionen.  Unter  Summe  von  f{x)  versteht  man  eine 
Funktion  (p(x),  deren  Differenz  gleich  f(x)  ist,  also 

(p{x-\-  h)  —  (p(x)  =  fix). 

Diese  Gleichung  bedeutet,  daß  der  Ausdruck  (p{x  -\-li)  —  (p{x) 
so  transformiert  werden  kann,  daß  er  in  f{x)  übergeht.  Dieselbe 
Relation  wird  auch  bezeichnet  durch 

^{x)  =  2:f{x). 

Es  ist  vorausgesetzt,  daß  f{x)  nicht  in  der  Form  (pix  +  K)  —  (p{x) 
gegeben  ist.     Sonst  wäre 

fp{x)  =>  2[(p{x  +  Ä)  —  (p{x)\ 

eine  reine  Identität,  die  keinen  Nutzen  haben  könnte. 

Wir  werden  untersuchen,  ob  es  möglich  ist,  daß  eine  andre 
Funktion  tl>{x)  dieselbe  Eigenschaft 

^{x  +  h)-'4>{x)  =  f{x) 

hat.     Wenn  wir  i^ix)  —  (p(x)  mit  F{x)  bezeichnen,  so  ist 

ij)(x)  =  cp(x)  +  F(x) 
und  nach  Nr.  3 

Ai(f(x)  =  ^(p(x)  +  £^I{x). 
Es  ist  also 

LF{x)  =  0     oder    F{x  +  /^)  =  F{x). 

Beide  Funktionen   i^ix)   und   (p{x)   unterscheiden   sich   also   um 
eine  periodische  Funktion  mit  der  Periode  h. 
Ist  ^{x)  eine  beliebige  Funktion,  für  die 

^{x  4-  h)  =  ^{x)     oder     A  ^{x)  =  0, 
so  ist 

t^[(p{x)  +  ^(x)]  =  A<p{x)-\-^  ^{x)  =  f{x). 

Sellwanoff,  Differenzenreolinang.  3 
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Also  haben  alle  Funktionen,  deren  Differenz  gleich  f(x)  ist,  die  Form 
(fix)  +  ^{x), 

wo  Q{cc)  eine  willkürliche  periodische  Funktion  mit  der  Periode  h  ist. 
Um  die  Aufgabe  zu  vereinfachen,  werden  wir  von  nun  an  voraus- 
setzen, daß  X  nur  Werte  der  Reihe 

a,     a  +  h,    a  -\-  2h,  . . . 

annimmt.     Dann  ist  0(x)  eine  willkürliche  Konstante,  die  wir  mit  C 
bezeichnen  werden.     Es  ist  also 

2Jf(x)  =^  (fix) -i- C. 

Dieser  Ausdruck  heißt  unbestimmte  Summe,  weil  der  Wert 
von  C  ganz  beliebig  ist. 

26.  Summatiou  einfacher  Funktionen.  Aus  jedem  Satz 
über  die  Differenz  (Nr.  2  —  6)  können  wir  einen  Satz  über  die  Summe 
herleiten. 

Aus  der  Relation  Ax  =  h  folgt  A^  =  1   und 


i[71  =  f +  0. 


Weil 


Ax{x  —  h){x  —  2h)  ...{x  —  n—  Ih)  =  nhx(x  —  h)...(x  —  n  —  2h) 
ist,  so  ist 

A  -^x{x  —  h)(x  —  2h)  ...(x  —  n—  Ih)  =  x(x  —  h)...(x~n  —  2h) 
und 

I]x(x-h)  (x-2h)...  (x-n-2h)  =  ^x(x-h).  ..{x  -n-lh)  +  G. 

Damit  man  unter  dem  Summenzeichen  n  Faktoren  hat,   ist   es 
zweckmäßig,  n  durch  w  +  1  zu  ersetzen.     Man  findet 

I]x{x  —  h)(x  —  2h)  ...{x  —  n  —  lh)=     j-iife^C^  —  h)...{x  —  nh)  -f  G. 
Wir  hatten  die  Formel  (Nr.  2) 

=  —  nh 


x{x-\-h)...{x-\-n-lh)  x{x-\-h)...{x  +  nh) 

Im  Nenner   der   rechten   Seite    sind   mindestens    zwei   Faktoren, 
weil  n^l  ist. 

Aus  dieser  Formel  folgt 

2;___L____  =  _  A 1  .   r? 

x{x  +  h)...{x-\-nh)  nh    x{x-\-h) . .  .{x^-n-lh) 

Die  Funktion  unter  dem  Summenzeichen  hat  im  Nenner  mindestens 
zwei  Faktoren. 
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Wenn  wir  n  —  1  statt  n  einsetzen,  so  entsteht  die  Formel 

2 '—= ,-?L^ '      _     +C 

x{x-\-h)...{x-{-n-lh)  (w-l)Ä    x{x  +  h) . .  .(x-\-n-2h) 

Hier   ist   w  >  1.     Dem   Wert  n  =  1    entspricht   eine   besondere 

1 


transzendente  Funktion,  die  wir  mit  U—  bezeichnen 

'  X 

Die  Relation 

Am^=m^(w^— 1)     oder     A -^^ — - 

V  j  w*  —  1 

gibt  die  Formel 


==  m^ 


2;w*  =  — — -  +  G. 
mfi  —  l 


Es  ist,  wie  wir  wissen  (Nr.  6), 


Wir  setzen 


A  cos  w  =  —  2  sin  -^-  sin 

u  ==  ax  -{-  ß 


(«  +  ^) 


,    Au 

u-\-  -j-  =  X 


und  wählen  a  und  ß  so,  daß 

wird.     Aus  der  Bedingungsgleichung 
folgt 

Deswegen  ist 


ax  +  ß  -{-  —  =  x 
a=l,     i3  =  -^- 


Acos 


(^-1) 


2  sin  —  •  sin  a;. 


Daraus  folgt 


und 


hm 


2  sin 


2 

cos 


=  siniC 


Usinx  =  — 


i-l) 


2  sin  — 


c. 


Ganz  analog  findet  man  die  Formel 


27co8rc 


+  a 


2  sin 


27.    Summatiou  zusammengesetzter  Funktionen.     Wenn 
Ä  ein  konstanter  Faktor  ist,  so  ist 

AÄ(p{x)  =  ÄAq)(x) 

2][AAg){x)]  =  Ä(p(x)  +  C. 

3* 


oder 
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Wenn  wir  Acp(x)  mit  f(x)  bezeiclmen,  so  entsteht  die  Formel 

Die    Konstante    C    ist    schon    im    Ausdrucke    2f(x)    enthalten. 
Wenn  wir  eine  Reihe  von  Gleichungen  haben 

so  ist 

Daraus  folgt 

^[fl(.Oc)  +  f,{x)-\----+  fm(x)]  =  (p,(x)  +  g),{x)+---+  (pm{x)  +  C 

oder 

^iftiP^)  +  /"sC«^)  +  •  •  •  +  fm(x)]  =  Zf,{x)  +  27/2(0;)  +•■•  +  2fm(x). 

Mit  Hilfe  dieser  Formel  ist  es  leicht,  ganze  rationale  Funktionen    ^ 
zu  summieren. 

Es  sei 

f(x)  -  A  + A^  +  A^^^i^^^  +  - 

,      j    (x  —  a)(x  —  a  —  h)(x  —  a  —  271) . . .  (x ~  a  —  n  —  Ih) 

Dann  ist 

Zf(x)  =  0 4-  A ^~" 4- Ä  (^~<^)(^-<^-^)  1  ^4  (^-«)(^-«-^)(a?-«-2fe) 

_L  j    (ic  —  a) (a;  —  a  —  Ä.) .  . .  (aj  —  a  —  wÄ) 
"'  1"^  1  •  2  •  3  .  . .  (w  +  1)  •  Ä 

Die  Summe  einer  ganzen  Funktion  n^^^  Grades   ist  also 
eine  ganze  Funktion  vom  (n  +  1)*®"*  Grade. 
Für  h=l  ist  z.  B. 

ic^  =  rr  +  x{x  —  1),    a;'  =  ic  +  3a;(rr  —  1)  +  x{x  —  l){x  —  2) 

und  folglich 

Zx^  =  C  +  |^a;(a;  -  1)  +|a;(a;  -  l)(a;  -  2), 

2;ä;3  =  G-\-\x(x  -  1)  +  rr(a;  -  l){x  -  2)  +  la:(a;-  l)(:»-2)(a;-3). 
Es  ist  auch  leicht,  gebrochene  Funktionen  von  der  Form 

m 

{x-\-a){x-\-a-\-li) . .  .{x-\-a-\-nh) 
zu  summieren,  wo  f{cc)  eine  ganze  Funktion  ist. 

28.  Partielle  Summation.    Aus  der  Formel 

A[(p{x)  •  ■^{x)'\  =  (p{x)  •  AiIj(x)  +  g)(x  +  h)  •  Ag)(x) 

^  ^     2[(p{x)  •  Acp{x)  +  if(x  +  h)  '  A(p(x)]  =  (p(x)  •  T^(a;)  +  C 

oder 

2;9)(a;)  •  A^(«)  +  I]<p{x  +  Ä)  •  A9)(a;)  =  (p(x)  •  t(x)  +  C. 
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Diese  Gleichung  liefert  die  Formel  der  partiellen  Summation 

E(p{x)  •  Af(x)  =  <p(x)  •  i{f(x)  —  Hip^x  +  h)  •  Aq){x). 

Die  willkürliche  Konstante  erscheint,  nachdem  die  Summation 
auf  der  rechten  Seite  ausgeführt  ist.  Mit  Hilfe  dieser  Formel  kann 
man  das  Produkt  (p(x)  •  f{x)  summieren,  wenn  man  einen  Wert  von 
2Jf(x)  kennt,  und  wenn  die  Operation  2]ijj(x  -\-  h)  •  A9(a;)  ausführbar 
ist,  wo  xlf(x)  =  2Jf(x)  ist. 

Wir  wollen  z.  B.  2JS^-  x  für  h  =  1  ausrechnen. 

Da  ein  Wert  von  2^3^  gleich  —3*  ist,  so  ist 

2;3^-  a;  =  2;a; .  A  (I  •  3^)  =  a;  •  ^  •  3*-  2:|  •  3^+^  •  Ax 

Es  wäre  nicht  zweckmäßig,  die  Summation  beim  Faktor  x  anzufangen. 
Da  Ux  =  -^ — -  ist,  so  wäre  dann 

^3*-  X  =  2; 3^-  A ^^-=^  =  3^ .  ^fcl)  _  2:^^+J-)^  ^3x 

A  A  A 

=  ^  S^x(x  -  1)  -  U(x  +  l)x  •  3% 

und  die  gegebene  Summe  würde  auf  eine  kompliziertere  zurückgeführt. 
Noch  ein  Beispiel: 


2x  sin  a;  ==  Ux-  A 


_cos(a;-|) 


2sm-       _, 


=  X 


-cos(x--) 


2  sin  — 
2 


—  cos 


Z— ^ Ax  = j^x  cos  (a^  -  2)  +  -^ ^8  -81110;+  G. 

2sin-  2sin-  (2sin-j 

29.  Eigenschaften  der  bestimmten  Summen.     Wir  haben 

schon  gesehen  (Nr.  25),  daß,  wenn 

(p{x  +  Ä)  -  (p{x)  =  fix),      t(x  +  h)  -  rlf{x)  =  fix) 

und  X  eine  Zahl  der  Reihe 

a,     a  -{-h,     a  +  2h,  . . . 

ist,  (p{x)  und  iIj(x)  sich  durch  eine  Konstante  unterscheiden 

(p(x)  =  t{x)+C. 
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Subtrahiert  man  voneinander  zwei  Werte  dieser  Funktionen  für 
x  =  a-i-n}i  und  x  =  a-{-mh,  so  fällt  C  weg.     Man  findet  die  Grleiehung 

q}{a  -\-  nh)  —  (p{a  +  mh)  =  -il^^a  -f  nh)  —  ■^{a  +  mh). 

Der  Ausdruck  (p(a  -\-  nJi)  —  (p{a  -\-  mJi)  bleibt  also  ungeändert, 
wenn  man  (p{x)  durcb  einen  andern  Wert  von  2f{x)  ersetzt.  Dieser 
Ausdruck  heißt  bestimmte  Summe  und  wird  bezeichnet  durch 

aA-nh 

y^,  f{x)  =  (p{a  +  nh)  —  cp(a  +  mh). 

a-\-mh 

Die  Zahlen  a  -}-  mh  und  a  +  nh  werden  Grenzen  dieser  Summe 
genannt. 

Aus  dieser  Definition  folgt 

a-j-mh  a-{-(m-\-l)fi 

^/■(^)  =  0,  ^f{x)  =  fia  +  mh). 

a-{-mh  a-\-mh 

Wenn  c  eine  Zahl  a  -{-Ich  ist,  die  zwischen  a  und  h  =  a  +  nh 
liegt,  so  ist 

6 


oder 


^f{x)  =  cp{b)  -  q>{a)  =  [cp{h)  -  cp(c)-]  +  [_cp{c)  -  9(a)] 


^m=^m-^^m' 


Wenn  man  diesen  Satz  mehrmals  anwendet,  so  findet  man 

&  a-\-h  a+ih  b 

^m = ^f{^) + ^m + •  •  • + ^i^)' 

Daraus  folgt  eine  wichtige  Formel 
f{a)  +  f{a  +  h)  +  /■(«  +  2ä;  +  ...  +  fQ)  -h)^  (p{h)  -  tp{a). 

Hier  ist  f(x)  die  gegebene  Funktion,  (p(x)  irgend  eine  Funktion, 
deren  Differenz  gleich  f{x)  ist,  — r—  eine  ganze  positive  Zahl.  Diese 
Formel  kann  nur  dann  nützlich  sein,  wenn  (p{cc)  nicht  die  Form  hat 

f{a)  +  /-(a  4-  /O  +  /"(«  +  2  /^)  +  •  •  •  +  f{x  -  h). 

30.  Berechnung  einiger  bestimmten  Summen.     1)  Es  sei 

s  =.  a -\- {a  +  h) -\- {a  +  2h) -\ h  («  +  n  -  ^h). 

Wenn  wir  a  -\-  nh  mit  h  bezeichnen,  so  wird 

s  =^x. 
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Ein  Wert  der  unbestimmten  Summe  2x  ist  (Nr.  26) 

_  x{x-h)^ 

deswegen  ist 

s  =  ^[h{b  -  h)  -  a{a  -  h)] 

oder  nach  einigen  Reduktionen 

n 


s  =  -^[a  +  {a  + n-lh)]. 


13 


2)      1 . 2  •  3  +  2  •  3  •  4  +  •  •  •  +  10  •  11  •  12  =^x{x  -l)(x-  2) 


=  |.13121110. 

23 

3^    J_4.J_4.J__L      _L    ^    _ y     1      =  1  _  Jl  =  H 

>*       l-3"^3.5"^5-7~'"'"     "^21 -23      -^  a;.(a;  +  2)         2        46        23" 

4)  l  +  g  +  g«.+  --4-2"-^.=2V=^- 

0  ^ 

5)  1  +  cos  h  +  cos  2h  -\ +  cos  (w  —  l)h  =  ^,  cos  a; 

0 

.     2w  — 1  ,    ,     .     Ä 
Bin  — - —  Ä  -I-  sin  - 

2  sin  — 
2 

Wenn  wir  von  beiden  Seiten  —  abziehen,  so  finden  wir 

.    2«  — 1, 
sm  — — —  h 
1  2 

—  +  cos  Ä.  4-  cos  27i  +  •■•  4-  cos  (w  —  1)7*  = r — 

2  sin  — 
2 

«A  ^  2w  — 1- 

^  cos  —  —  cos  — - —  h 

2  2 


6)  sin  Ä  +  sin  2J^  H 1-  sin  (n  —  l)h  =^^ sin  x  = 


2  sin  — 
2 


7)  Wir  wollen  jetzt  mit  Anwendung  von  Formeln  in  Nr.  27  die 
Summen  von  Potenzen  natürlicher  Zahlen  berechnen. 

X 

0  +  l  +  2+.+(rr-l)  =^0:  =  ^^. 

0 

X 

0'+V-{-2'+-   +(x-iy=^x'  =  ^x{x-l)-\-^x(x-l)ix-2) 

0 

a;(g-l)(2a?-l) 
""  6 
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03+   13+  23+  .  ..  +  (a;  _   1)3  _2'^3_   1  ^(^  _   1)  ^  ^(^  -!)(«;-  2) 

0 

+  i^(^  - 1)  (^  -  2)  (..  -  3)  =  p^=i)]^ 

Daraus  folgt  die  bemerkenswerte  Relation 
1»  +  2»  +  33  +  . . .  +  (a;  -  1)3  =  [1  +  2  4-  3  +  •  •  •  +  (a;  -  1)]^  * 

Wenn  wir  — —  und  statt  rr*  und  a;^  summieren,  so  finden  wir 

X 

^^X  X^  1  x{x  —  l) 

0 

y»_^  _     a;»      _  £    _^   ,   J^     _  a;(a;-l)[a;-^j 

^1-2        1.2-3         21.2"'~12^~  1.2-3  ' 

0 

"y    a;" ^ £        a;»  ^    _a^  _  a;Va!-l)» 

-^l-2-3  ~  1.2-3-4         2'l.2.3~^12*1.2~1.2-3-4* 
0 

Im  nächsten  Kapitel  werden  wir  die  Funktion  betrachten 

-^l-2-3...(w-l)  ^  ^"^^)' 
0  ^  ^ 

die,  abgesehen  von  einem  konstanten  Faktor,  von  Jacob  Bernoulli 
eingeführt  ist. 

Zweites  Kapitel. 
Die  Jacob  Bernoullische  Funktion. 

31.  Bestimmung  der  Koeffizienten.     Nach  dem  in  Nr.  27 
bewiesenen  Satze  ist  die  Jacob  Bernoullische  Funktion 


(1)  ^1.2.3^... (n:ri)  ==  9'«  W 


eine  ganze  Funktion  n^^^  Grrades,  die  wir  in  der  Form 

(2)  yn(^)=ArT^+Ai.2!'!'.(n-i)+^2i.2!'."(n-2)+---+-^--^^ 

schreiben  werden.  Das  konstante  Glied  fehlt,  weil  (pn(0)  =  0  ist. 
Denn  in  Nr.  29  haben  wir  gesehen,  daß  die  bestimmte  Summe  gleich 
Null  ist,  wenn  ihre  Grenzen  einander  gleich  sind.  Wir  wollen  die 
Koeffizienten  Äq,  A^,  A^,  . .  . ,  Än—i  bestimmen. 

*  Vgl.  Encykl.  der  Math.  Wiss.,  Bd.  I,  p.  579. 
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Da 

9)„(a;+  l)=^i.2...(n-i)=^^l.2...(w-i)+-2;r7 


^   (w-1)      ^l.2...(w-l)   '  -^1.2...(«-l) 

ist,  so  besteht  die  Relation 

/» 1 

(3)  ^)n{x  +  1)  —  (fnix) 


1-2...  (w-1) 
die  nur  für  ganze  positive  x  bewiesen  ist.     Die  ganze  Funktion 

Fix)  =  gp„(a;  +  1)  —  ifnix) 


1  •  2 ...  (w  - 1) 

vom  nicht  höheren  als  yC-^'°-  Grade  verschwindet  für  jeden  ganzen 
positiven  Wert  von  x.  Also  ist  F{x)  identisch  gleich  Null  und  die 
Relation  (3)  besteht  für  jedes  x. 

Wenn  eine  andre  ganze  Funktion  ^(a;)  die  Eigenschaften 

hat,  so  ist  i\}{x)  mit  ^)n{x)  identisch,  weil  beide  Funktionen  vom 
»*®°  Grade  (Nr.  4)  sind  und  außerdem  die  Gleichungen 

^(0)  =  9,(0),     ^(l)  =  9p(l),     ^(2)  =  qp(2),  ...,1^(»^)  =  qp(,^) 

bestehen.  Die  Werte  ^(l),  9? (2),  .  .  . ,  qp(w)  bestimmen  sich  sukzessive 
aus  der  Gleichung  (3). 

Durch  die  Differentiation  von  (3)  findet  man 

^«(^  +  1)  -  ^'n{x)  =  ,.2.''"(^_2)- 

Nach    der    Gleichung    (2)    ist    qp',(0)  =  J.„_i.       Deswegen    ist 
^\{x)  —  An—\  eine  ganze  Funktion,  die  für  ic  =  0  verschwindet  und 

deren  Differenz  gleich  ^7^ 1  _<^\  ist.    Folglich  ist  für  n>  \ 

,  9>«(^)  —  ^-1  =  ^>n-\{x) 

oder 

(4)  tf\,{x)  =  (pn-lix)  +  Än-l- 

Aus  dieser  Relation  schließen  wir,  daß 

ist,  und  daß  die  Koeffizienten  J„,  Ä^,  . .  . ,  An—%  hier  dieselben  sind, 
wie  in  der  Gleichung  (2). 

Wenn  wir  in  (3)  x  =  0  setzen,  so  finden  wir  9P,,(1)  =  0  oder 

^^^    1 .  2  .". .  n  +  1.2...  (n-1)  ^  1-2...  V- 2)  +  *  *  *  +  iV  +  ^'— ^  ^  ^- 
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Diese  Relation  gibt  die  Mögliclikeit,  Än—i  zu  berechnen,  sobald 
die  Koeffizienten  Äq,  A^,  .  .  .,  An—^  bekannt  sind. 
Aus  der  Definition 


.X 

folgt,  daß  ccia^-i)    \^        1 

ist.     Deswegen  ist 

A=  1;      A  =  ~"  2' 

Die  Gleichung  (5)  für  n  =  3,  4  und  5  gibt 
iiT2~3  "^  IT2  +  ^  ^  ^^ 

l-2-3-4-5"^l-2-3-4"^12-3"^l-2"'"4 

Daraus  findet  man 

^2  =  I2'     ^3  =  0,     A^=  —  ^^- 
In  dieser  Weise  kann  man  beliebig  viele  Koeffizienten  berechnen. 

32.  Eigenschaften  der  Koeffizienten.     Wenn   man   in  der 
Relation 

9'2;fc+l(l  +  ^)  -  9C'2,  +  i(^)  =  i.2.'^3...2fc 

X  durch  —  X  ersetzt,  so  findet  man 

(ö)  9'2Ä+i(l  -  ^)  -  ^'ait+iC-  ^)  =  i.2.l...^k' 

Wenn  man  die  Bezeichnung 

einführt,  so  ist 

-  9'2.+i(-  ^)  =  -  ^'g.+Jl  -  (^  +  1)]  =  ^i^  +  !)• 
Die  Gleichung  (6)  nimmt  die  Form  an 

■F(^  + 1)  - -z^i^)  =  rÄ*- 

Da  aber  F(p)  =  0    ist,   so    ist  nach   dem   in  Nr.  31   bewiesenen 
Satze  F(x)  =  gpgxr+iC^)  ^^^^ 
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Vergleicht  man  die  Ableitungen  der  beiden  Seiten  dieser  Gleichung, 
so  findet  man  mit  Berücksichtigung  der  Gleichung  (4) 

-  [9P2i(l  -  a;)  +  ^2*]  =  -  [(p^,(x)  +  Äu] 
oder 

(8)  9>2,(1  -^)  =  9'2i(4 

Für  Z:  =  1  ist  die  Ableitung  von  (p.^f.{x)  keine  Jacob  Bernoullische 
Funktion,  weil  die  Definitionsgleichung  (1)  in  Nr.  31  für  n  =  l  keinen 
Sinn  hat.     Aus  der  Differentiation  von  (8)  folgt  für  Z;  >  1 

-  [^'si-iCl  -  ^)  +  -^ai-i]  =  9'2i-i(^)  +  ^2i-i 
oder  nach  (7) 

—  -0.2*  — 1  =  -a.2A  — 1. 

Es  ist  also 

(9)  A^k-i  =  0  für  Ä;  >  1. 

Wir  haben  schon  in  Nr.  31  gesehen,  daß  A^  nicht  Null  ist,  und 
daß  ^3  =  0  ist.    Jetzt  überzeugen  wir  uns,  daß  A^  =  0,  ^  =  0, . . .  sind. 

33.   Das    Vorzeichen    von    <pn{pc)    im    Intervall    (O,  1). 

Wenn  man  in  (7)   x  =  -    setzt,  so  findet  man 


2 


^2.  +  l(|)=0. 


Wir  wollen  untersuchen,  ob  für  ^>  1  die  Gleichung  ^'gi-i-iC^)'^^ 
noch  eine  Wurzel  im  Intervall  (0,  1)  hat. 

Es  sei  vorausgesetzt,  daß  diese  Gleichung  außer  0  und  1  noch 
die  Wurzeln  %  und  a^  besitzt,  die  den  Ungleichheiten 

0  <  «j  <  «2  <  1 

genügen.  Hier  darf  a^  oder  a«  gleich  —  genommen  werden.  Nach 
dem  Roll  eschen  Satze  hat  ^a^-iiC^)  ==  0  mindestens  drei  Wurzeln 
?)i,  &2  ^^uid  &3,  die  den  Ungleichheiten  genügen 

0<\<a^<\<a^<l.,<  1. 

Die   Gleichung    9'2/t+i(^) ""  ^   ^^^    mindestens  zwei  Wurzeln   c^ 
und  Cg,  die  zwischen  0  und  1  liegen. 
Es  ist  aber 

9'2*  +  l(^)  =  9'2i-l(^). 

weil  ^2it  — 1  =  0  ist.  Die  Gleichung  <p^j^_-^ix)  ^0  hat  außerdem  die 
Wurzeln  0  und  1. 

Aus  der  Voraussetzung,  daß  die  Gleichung  92*+i(^)  ^  ^  zwischen 
0  und  1  zwei  Wurzeln  hat,  folgt  also,  daß  die  Gleichung  9'2i_i(^)  =  0 


44  Zweiter  Teil.    Summen. 

in  demselben  Interyall  auch  zwei  Wurzeln  hat.  Wenn  wir  diesen 
Schluß  mehrmals  anwenden,  so  finden  wir,  daß  (Pq{x)  zwischen  0 
und  1  zwei  Wurzeln  hat.  Das  ist  aber  unmöglich,  weil  die  Funktion 
3*^^  Grades  93  (ic)  nicht  für  vier  verschiedene  Werte  von  x  (0  und  1 
inbegriffen)  verschwinden  kann. 

Es  ist  also  bewiesen,  daß  die  Gleichung  9'2/fc4-i(^)"^^  zwischen 

0  und  1  nur  eine  Wurzel  ^  ==  ^  hat. 

Wir  wollen  jetzt  die  Gleichung  (p^^(x)  =  0  untersuchen.  Wenn 
diese  Gleichung  zwischen  0  und  1  die  Wurzel  a  hat,  so  folgt  nach 
dem  Rolleschen  Satze,  daß  die  Gleichung  qp^^(ic)  =  0  mindestens 
zwei  Wurzeln  zwischen  0  und  1  hat.  Das  ist  aber  unmöglich,  weil 
^2i(^)  "^  9^2^— 1(^)  ^^^  ^^^  ^^^  eben  bewiesen  haben,  daß  <P2k—ii^) 
nicht  für  zwei  Werte  zwischen  0  und  1  verschwinden  kann. 

Es  besteht  also  der  Satz:  Für  alle  Werte  von  x  zwischen 
0  und  1  hat  die  Funktion  ^PjäC^)  dasselbe  Vorzeichen. 

Um  dieses  Vorzeichen  zu  bestimmen,  betrachten  wir  die  Gleichung 
9^2* +i(^)  ^  ^>  ^^®  befriedigt  wird  für  x  =  0  und  x  =  -^'  Nach  dem 
Rolleschen  Satze  hat  die  Gleichung 

"P^k+ii^)  =  9P2i(^)  +  ^2^=  0 
mindestens  eine  Wurzel  x  =  a,  die  zwischen  0  und  —  liegt. 

Aus  der  Relation 

folgt,  daß  Ä2jc  nicht  gleich  Null   sein   kann,    weil   (p^j^ix)   nicht   im 
Intervall  (0,  1)  verschwinden  kann,  und   daß  in  demselben  Intervall 
(f^j^ipc)  das  entgegengesetzte  Vorzeichen  hat,  als  A^k- 
Es  ist  aber 

9'2^(^)  =  1.2...2Z;  +  ^1 1-2...(2Ä-1)  +  ^  l-2...(2A-2) 


-2i  — 4 


a3' 


+  ^M.2...(2Ä-4)+'*"  +  ^^^-2r.2 

Nach   einem   Satze  [der  Algebra   hat   (p^{x)    für   kleine    positive 

.  .  .  x^ 

Werte   von   x   dasselbe   Vorzeichen    wie    das    Glied   Azk— 27-70  ^^^^ 

wie  Aik-s. 

Die  Koeffizienten  A2k—2  und  A^k  haben  also   verschiedene  Vor- 
zeichen.    Da  aber  ^9  =  -^  >  0  ist,  so  ist 

'S         12  ' 

A,<0,    A>o,    ^8<0,... 
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und  für  0  <  a;  <  1  besteht  die  Ungleichheit 

Die  Funktion  (p2k+i(^)  wechselt  ihr  Vorzeichen  im 
Intervall  (0,  1),  weil  x  =  ^  eine  einfache  Wurzel  der  Gleichung 
9?2j,  j(a;)  =  0  ist.  Um  diese  Behauptung  zu  beweisen,  nehmen  wir 
an,  daß 

sei.     Die  Funktion  (p2j^{x)  +  -^2*  verschwindet  für  x  =  cc  und  rc  =  —  • 
Nach  dem  Roll  eschen  Satze  soll 

verschwinden   für    einen  Wert    von   x   zwischen    a  und  — ;  was   un- 
möglich ist. 

34.  Bemoullische  Zahlen.  Die  Koeffizienten  Ä2,A^,Äq,  ... 
werden  oft  durch  die  Zahlen  B^,  JB^,  B^, . ..  ausgedrückt,  die  man 
Bemoullische  Zahlen  nennt.  Es  sind  Zahlen,  die  mit  dem  Vor- 
zeichen +  oder  — ,  als  Koeffizient  von  x  in  der  ersten  Potenz,  in 
den  Ausdrücken  für 

X  X  X 

y^i  ^  >  Xi  ^J  X>  X  ,  .  .  . 

auftreten.  <>  «>  <> 

Es  ist 

X 

^x'=l-2g>,(x)  =  l-2A^x  -  x' +  ^. 
0 

X 

^a^  =  l-2  ^'4<p^{x)  =  12-3'4:Ä^x+-" 
0 

X 

^x'=l-2'3'4:'b'Q(p,{x)  =  l'2-^-4c'6-6Ä^x  +  -" 
0 
Es  wird  daher  gesetzt 

12^  =  ^1,     1 -2.3.4^4  =  -.^2,     1-2.3-4   5-6J«  =  53, ... 

Allgemein  ist 
oder  l'2-..2k-Ä,,  =  i-iy-^B, 

Die  ersten  Bernoullischen  Zahlen  sind 
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Die  Tabelle  der  62  ersten  B  findet  man  bei  J.  C.  Adams  im 
Journal  für  Mathematik,  Bd.  85  (1878),  S.  269—272.  Die  ersten 
20  Zahlen  sind  im  Lehrbuch  von  A.A.  Markoff  abgedruckt  (S.  125). 

Auf  die  arithmetischen  Eigenschaften  der  Bernoulli sehen  Zahlen 
gehen  wir  hier  nicht  ein.* 

36.    Entwicklung  von  (f^^ix')  in  eine  trigonometrische 

Beihe.    Aus  der  Theorie  der  trigonometrischen  Reihe  (L.  Dirichlet, 
Werke,  Bd.  I,  S.  119  —  132)  folgt,  daß  für 

0^a;<l 
die  Entwicklung 

(p^J^x)  =  -^  «0  +  '^i  ^^^  ^^^  +  ^2  COS  4jric  H \-  ün  cos  2n^x  -\ 

+  &i  sin  27tx  +  &2  sin  4^%x  -\ f-  &„  sin  2n%x  -\ 

besteht.     Die  Koeffizienten  haben  die  Werte 

a„  =  2  I  (p^j^ipc)  cos  2'yf:tx-dx,     h»  ==  2  1  (p^f.(x)  sin  2nnx  •  dx. 

0  0 

Wenn  wir  x=^l  —  y  im  Ausdruck  für  &„  einsetzen,  so  wird 
wegen  der  Relation  (8)  in  Nr.  32 


&«  = 

1 
=  -  2  /  (p^j^{y)  sin  2njty  ■  dy. 

Folglich  ist  hn  =  ■ 
Der  Wert  für 

0 

—  tn  oder  hn  =  0. 

1 

«0  =  2  /  (p^j^(x)dx 

0 

läßt  sich  aus  der  Relation 
bestimmen.     Man  findet 


\%=- 1  (p2k(.^)^^  =  -  ^2*. 


Um  ttn  zu  berechnen,  werden  wir  das  Integral 

1  1 

/  9^24  (^)  COS  2n3tx  •  dx     durch       /  (fs^-^i^)  ^os  2nnx  •  dx 

0  0 

ausdrücken.     Zu  diesem  Zweck  gebrauchen  wir  die  Formel 

*  Diesem  Gegenstand  ist  der  dritte  Abschnitt  des  Werkes  von  Herrn 
L.  Saalscliütz  gewidmet  (Vorlesungen  über  die  Bernoullisclien  Zahlen  .  .  . 
Berlin  1893). 
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(11)  I  (pip"dx  =  (fxl;'  —  q)'t -{-  ljljcp"dx, 

die  eine  Folge  der  zweimaligen  partiellen  Integration  ist. 
Setzt  man  in  (11) 

9  =  g?2j  (x),  ^"  =  cos  2nycx, 

so  findet  man  für  Ä;  >  1 


,  f       Sin  2nnx 

if  =  —        f 


ff 


9>2k-^i^)  +  ^2*-2;    t  = 


Da  aber 


ist,  so  wird 


2njt 
cos  2nnx 

{2nny" 


ß 


cos  2n7tx  •  dx  =  0 


0  0 

Wenn  man  diese  Formel   Qc  —  1)  mal  anwendet,   so    erhält   man 


J  "PiM 


cos  2nnx  ■  dx 


± 

(2w 


x)  COS  2n7cx  •  dx. 


Es  sei  jetzt  in  (11) 
Daraus  folgt 


,  1 

9»  =^-V 


^"  =  cos  2n7tx. 


, ,       sin  2'njex 


9"=1, 


^^- 


und 


/' 


qP2(^)  cos  2n%x  ■  dx  = 


2WÄ 
cos  2n7tx 
(2n»)* 

1 


(2wjt)« 


Man  findet  also 


««  = 


(-1)* 


(2««)" 

Wenn  man  in  der  Entwicklung 

.'Pski^)  °^  ^.^0  +  ^1  cos  2«^;  +  öTg  cos  4Äa;  +  •  •  • 
o;  gleicliNull  setzt,  so  findet  man  die  bemerkenswerte  Formel  von  Euler 
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Diese  Formel  wird  später  benutzt,  um  die  Konvergenz  einiger 
Reihen  zu  untersuchen. 

Aus  den  Werten  für  A^  und  Ä^  folgen  die  Formeln 

1  j_  jL4.i-_i_i__L     -'^ 

■^  "•"  2*  "^  3«  "^  4*  "•"  ■  '  ■         6  ' 
-^  "^  2*  "•"  3*  "'"  4*  "^  90 

Dieses  Kapitel  habe  ich  unter  dem  Einfluß  vonW.G.Imschenetzky 
(Kasaner  Universitäts- Nachrichten  1870  und  französisch  in  J.  Hoüel, 
Cours  de  calcul  infinitesimal,  Paris  1888,  Bd.  1,  S.  476)  und  Herrn 
N.  J.  Sonin  (Warschauer  Universitäts -Nachrichten  1888  und  französisch 
im  Journal  für  Mathematik,  Bd.  116  (1896),  S.  133)  bearbeitet. 
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Eulersche  Summationsformel. 

36.  Eulersche  Formel  für  ganze  Funktionen.     Es  sei 

^(^)  =-^or72^7^;r=l)  +^1  1.2!'..  (n-2)  +  •      +Pn-2^  -^Pn-1' 
Wenn  man  diese  Funktion  von  0  bis  x  summiert,  so  findet  man 


7 

0 


2jf(^)  ==Po<Pni^)  +  Pl<Pn-li^)  +  •  •  •  +  Pn-2'P2i^)  +  ^«-1^' 


Hier  treten  Bernoullische  Funktionen  ein.  Wir  wollen  dieselben 
vollständig  ausdrücken  und  die  Glieder  nach  den  Koeffizienten 
J-i,  J.2,  J-3,  ...  anordnen. 

In  dieser  Weise  entsteht  die  Formel 

X  

0 

+  A(i>0i.2!'!'.(n-l)+i>ll.2!".(n-2)+     -•  +  ^^-2^) 

+  ^  (^0i.2!'".(n-2)+^U.2.''..(n-3)  +  '--+^^-3^) 

-{-Än-2[po^+PlO0J 
+  Ä„-ip^X 
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oder 


(1) 


2 


0 


Das  ist  die  Eulersche  Summationsformel.     Sie  ist  genau  für 
ganze  Funktionen  vom  (n  —  l)*^'^  oder  niedrigeren  Grade. 


37.    Bestglied  der  Eulerschen  Formel.     Wir   wollen  die 
Funktion 


(2) 


X 

t(x)=^f(x)-  ff{x)dx-A,\f{x)-m-\-A,[r(x)-f'{0)-\- 
.•  ■  -An-,[r—'\x)-fi—\0)-\-Ä„-^[f^—'){x)-f(—\(S)-\ 


ausdrücken    in    dem    Falle,    daß    f{x)    keine    ganze    Funktion    vom 
(w  —  ly^^  Grade  ist. 
E  Die  Differenz  von  il){x)  ist 

Ai,{x)  =  f{x)  -  A  jf(x)dx  -  A^t^f{x)  -  A,Af'(x) 

0 

A^_,Afi—^\x)  -  A_iA/'(«-2)(a;). 

Aus  der  Integralrechnung  ist  bekannt,  daß 

AF{x)  =  Fix  +  1)  -Fix)  ^F'ix)  +  ^F"{x)  +  ^-F"'ix)-\- •■^ 


(3) 


+ 


1-2.3  ...m 


i^'W 


W+/r^ 


I^"*+^\x  +  u)du 


ist. 


Nachdem  wir  sukzessive 


X 

F{^)  =fmdx,     f^x\     fix), ...,  r-^)ix) 


0 

m  =  w, 


w  —  1,        n  —  2,  . . . ,         1 


gesetzt  haben,  so  verwandelt  sich  (3)  in 

Seliwanoff,  Differenzenreclmuiig. 
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Ltp{x)  =  f{x) 


1 


0 

1 

0 

1 
-^^"W-  ■  -1.2.  .^(n-2)  /'-"(^)-/ff  :(f_",/'"(^+») 


0 

1 


n-tf(—')(x)-rÄ„^l{l-u)fi-)(x  +  tl) 


oder 

A^(^)  =  -  9'2(l)/"(^)  -  9'3(l)/'"(^)  -  •  •  ■ 

0 

Da  die  Werte  ^'2(1),   9)3(1),...,  qp^(l)   gleich  Null   sind,   so  ist 

1 

AiIj(x)  =  —  j  (pjil—  u)f^''\x  +  W-)<?M. 
0 

Nach  der  Formel  (2)  ist  ^(0)  =  0  und  folglich 
1 

0 
0 

Es  wird  hier  nach  x  summiert  und  nach  u  integriert. 
Wir  werden  n  =  21i  setzen.     Wegen  der  Relation 

9'2i(l-w)  =  gP2^(w) 

ist  dann 

1 


(5) 


^{x)  =  -  f(p^^{u)du^p'^\x  +  u). 


Das  ist  das  Restglied  von  Jacobi  (Journal  für  Math.,  Bd.  12  (1834), 
p.  263  oder  Werke,  Bd.  6,  p.  64). 
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Wenn  man  in  (5)  die-Bernoullische  Funktion  (p^f^{u)  durch 
die  trigonometrisclie  Reihe  (Nr.  35)  ausdrückt,  so  findet  man  das 
Restglied  von  Poisson  (Pariser  Memoires,  Bd.  6  (1826),  p.  580). 
Wir  werden  die  Jacobische  Form  gebrauchen. 

Die  Funktion  9^^.(.w)  behält  ihr  Vorzeichen  im  Intervall  zwischen 
0  und  1.     Deswegen  folgt  aus  (5) 

1 

0  y 


wo  0  <  Ö  <  1  isi 
Da  aber 

so  ist 


und  folglich 

(6) 


1 

J  q>3,{u)du  = 


La* 


^{x)  =  Ä,,^f(^^)(x-[-e). 


Wir    haben    also    bewiesen,    daß    für    jede    Funktion    f(cc),    die 
gewissen  Stetigkeitsbedingungen  genügt,  die  Formel  besteht 


(7) 


X 

=^ß{x)dx  +  A,[f{x)  -  /-(O)]  +  Ä,[f{x)  -  /"(O)]  +  •  •  • 

+  A^Z^ir^'-'^ix)  -  /•(2^-3)(0)]  +  J.2,_i  [/■(2*-2)(a;)  -  /•(>*-2)(0)] 
i+  ^2*[/-(^*>(0)  +  /•(^*)(9  +  1)  +  •  •  •  +  /•'">(0  +  x-  1)]. 
Für  X  =  1  ist 


(8) 


m  =  ff{x)dx + A,[fi\)  -  /-(o)]  4-  A,[f{i)  -  rm + 

^>*-a[/-("-»)(l)  -  /'(«*-»)(0)] 
-f  ^2,_i[/'(2*-2)(l)  _  /•("-8)(0)]  +  ^>*r*)(0). 


Die   Glieder    mit   A^k—x   fallen    weg,    sobald   Ä>1  ist   (Nr.  32, 
Form.  9).     Für  Ä;  =  1  ist 

4* 
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X 

(9)  ■=^fmdx  +  A,\f{x)-f{0)-\ 

1+  Aif(P)  +  /""(ö  + 1) + r(ö  +  2)  + ... + /-"(ö  +  ^  - 1)] 

und 

(10)  m  =ß{x)dx  +  A,\f{l)  -  /-(O)]  +  AJ"ie). 

0 

38.  Andre  Form  des  Bestgliedes.  Wenn  die  Ableitungen 
/'(2*)(^)  und  /■(2i  +  2)(^)  für  Werte  von  t  zwisclien  0  und  x  ihre  Vor- 
zeicten  behalten  und  dieselben  Vorzeichen  haben,  kann  das  Rest- 
glied (6)  eine  andre  Form  annehmen. 

In  der  Formel 

/  (pc)"du  =  (pco'  —  cp' (o  -\-  j  G}(p"du 
setzen  wir 

X 

0 

Dann  wird 

X 

0 

X 

0 

und  folglich  ist 

1 


0 

1 


/»  X 

^,,+,{^)dn-^r''+%x^u). 
0 


0 

, .  Es  besteht  also  die  Relation 


1  1 

X         ^  /^  « 
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Da  aber 
1 


/' 


0 

so  ist 


f^^\x  +  u)du  =  /'(2*-i)(a;  +  1)  -  /•(8*-i)(a;)  =  A/'(2*-i)(a;), 
1 


«     0 


und  folglich 


/• 


0 

1 


Das  erste  und  letzte  Glied  dieser  Formel  haben  verschiedene 
Vorzeichen,  weil,  nach  dem  in  Nr.  33  bewiesenen  Satze,  (p^j.(u)  und 
9)3^  ,  2  (w)  im  Intervalle  (0,  1)  verschiedene  Vorzeichen  haben  und 
nach  unsrer  Voraussetzung  /'(**)(^)  und  f^^^—^)(t)  gleiche  Vorzeichen 
haben. 

Die  gefundene  Relation  hat  die  Form 

a  =  b  -{-  c, 

wo  ac  <.0  ist.     Aus  der  Gleichung 

a^  =  ah  -\-  ac 
folgt 

0<a2<a&,     0<|<1. 

Es  ist  also  ^, 

a  =  üO, 

wo  0  <  6  <  1  ist. 

Nach  dieser  Bemerkung  ist 

(11)  t(x)  =  e'Ä,,[f^^^-^)(x)-f^''-%o)i 

Diese  Formel  kann  statt  (6)  gebraucht  werden,  sobald  f(t)  gewissen 
Bedingungen  genügt.     (Vgl.  Markoff,    Differenzenrechnung,    S.  121.) 

39.  Zweite  Form  der  Eulerschen  Formel.  Wir  werden 
voraussetzen,  daß  für  alle  positiven  Werte  von  t  die  Ableitungen 
f"(t\  P'^if),  /(6)(0, .  . .  dieselben  Vorzeichen  behalten  und  daß  alle 
Ableitungen  /"'(^j,  /"(<), /'"(i), .. .  unendlich  klein  werden,  sobald  t 
unendlich  groß  wird. 
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In  diesem  Falle  hat  das  Integral 
1 

/»  CO 

0 

einen  Sinn.     Es  bedeutet,  daß  das  Integral 

0 

bei  unendlich  wachsendem  x  sich  einer  bestimmten  Grenze  nähert. 

In    der   Tat,    bei    wachsendem    x   vergrößert    sich    der    absolute 
Betrag  des  Integrals,  weil  neue  Elemente  mit  demselben  Vorzeichen 
zutreten.     Sobald  aber  x  hinreichend  groß  ist,  so  ist 
1 


ü 


<  1^2*1(1  +|/(^*-^>(0)i}, 


weil  hier  die  Formel  (11)   anwendbar   ist   und   nach    der    gemachten 
Voraussetzung  \f^^''~'^'^(x)\  beliebig  klein  sein  kann. 

Nach    bekanntem    Prinzip    hat     das    Integral    eine    bestimmte 
Grenze. 

Aus  der  Relation 

1 

0 

1 

/-»  x 

folgt  0  " 

1 

/>  CO 

0  « 

1 

/»  CO 

.  =  ^2.P*-^)(0)  +  /  (p,,  +  ,(u)du^f^^''+^^(t  +  U). 

«/  0 


0 

Da  aber 
1 


hp,,{u)du^fm(t  +  u)     und       U,,^,{u)du^p^+^){t  +  u) 

verschiedene    Vorzeichen    haben,    so    ist    nach    der    Bemerkung    im 
vorigen  Paragraphen 
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1 

wo  0  <  01  <  1  ist. 

Um  diese  Formel  zu  verallgemeinern,  genügt  es  zu  setzen 

wo  V  eine  neue  unabhängige  Variable  ist   und    x   einen   bestimmten 
Wert  hat. 

Nach  der  bewiesenen  Formel  ist 
1 

0 

oder 

1 

0 

wo  0  <  02  <  1  ist. 

Wenn  man  in  der  Euler  sehen  Formel  das  Restglied 
1 

durch 

1  1 

<P,,(u)du^f('^)(t  +  u)  +      cp^,{u)du^f^^^){t  +  u) 

0  «  0  * 

ersetzt,  so  findet  man 

X 
X  ^ 

^m  =  /  fmt  4-  C*  +  A,f{x)  4-  A,f{x)  +  •  •  • 

«        sr 

wo 

C.  =  -A/'(0)-^/"(0)--.. 
1 

/>  00 

0 

Jetzt    wollen    wir   beweisen,    daß    die  Konstante    C*   von   Je   un- 
abhängig ist. 
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Für  Ä;  =  1  ist 

*  1 

Mit  Benutzung  von  (12)  findet  man 

X 
X  >-» 

^m  =  c, + /  mdt + Aj{x)  4-  .4/'(^) + •  • 

nnd  ^  ^ 

2/(0  =  Ci + / /(o^<  +  A/'W  +  KAfip^y 

Der  Yergleicli  der  beiden  Formeln  gibt  die  Relation 

die  für  jeden  ganzzabligen  positiven  Wert  von  x  besteht. 

Bei    unendlich    wachsendem   x    sind   alle    Ableitungen   von   f{x) 
unendlich  klein  und  folglich  ist  Ck  =  C^. 

Es  ist  also  die  Formel  bewiesen 


(13) 


X 
X  ^ 

^f(t)  =  0  +     f(t)dt  +  AJ(x)  +^/''(^)  +... 

-j-A,,_,fi^^-^)(x)  +  eA,kf^''-'K^), 


wo  0  <  0  <  1  ist. 

40.   Der  allgemeinere  Fall  der  Eulerscheu  Formel.    Bis 

jetzt  haben  wir  die  bestimmte  Summe  von  der  Form 

m+m+m+----^f{x-i) 

betrachtet.     Nun  wollen  wir  zur  Summe 

F(a)-\-  F{a  +  h)  +  F{a  +  2h) +■"  +  F(b  -h), 
wo  h  ==  a-\-  xh  ist,  übergehen. 

^''''  F{a^th)^f{ty, 

dann  ist 

F{a)  =  m,     I{a  +  h)=^f(l),    F{a  +  2h)  =  f{2)r" 

F(b  -h)  =  fix  -  1) 
und  folglich 

F{a)  +  F{a-\-h)+-..  +  F(h-h)=m  +  f{l)+---\-f{x-V). 
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Wenn  wir  berücksictitigen,  daß 

so  nehmen  die  Formeln  (7),  (8),  (9)  und  (10)  die  Formen  an 

(F(ia)-\-F{a  +  h)  +  F(a  +  2h)  +  '"  +  F(b-K) 


(14) 


^rF(u)du+Ä,[F{h)-F{a)-]+A^h[F'{b)-F'{a)] 


+ 


(15) 


(16) 


+  ^2,/i2*[F(2*)(a+0/i)+i?'(2*)(a+l+0/^)+. .  .+F^*(h-h+eh)l 

a+h 

F(a)=~  rF(u)du+ÄlF(a-^h)-F(a)]+A^h[F'(a+h)-FXa)]+-  •  • 

+A,k-ih^'-''[F(^'-\a+h)-F(^^-\a)']+Äi,h^'F^''^){a+eh), 
(F(a)  +  F(a  +  Ä)  +  F{a -\- 2h) -\- -- -  +  F(h  -  h) 

b 

=  ^fF(uydu  +Ä,  [F(h)  -  F{a)\ 
+ilÄ'[jP"(a+ÖÄ)+JP"(a+l4-0Ä)  +  -"+i^"(&-Ä  +  ÖÄ)] 


und 

a  +  h 


(17)  F(a)  =  j  JF(u)du  +  ^1  [F{a  +  Ä)  -  i^(a)]  +  ^/^^jjr,,^^  _j_  ^-^^ 

a 

In  allen  diesen  Formeln  ist  }>  =  a  -\-  xh  und  a;  eine  ganze 
positive  Zahl. 

Wenn  die  Ableitungen  F'\i),  F^\t),  F(^){t),  .  .  .  für  t^a  die- 
selben Vorzeichen  behalten  und  die  Ableitungen  -^'(0,  F"(t),  F"'(t),.. . 
unendlich  klein  werden,  sobald  t  unendlich  groß  wird,  dann  besteht 
die  Formel 


(18) 


{F{a)-\-F{a  +  h)  +  F{a-\-2h)-\-'-'  +  F(b-h) 
b 

=  C  +  irF{u)du-\-Ä,F(b)-{-Ä,hF'{b)-\---' 
+  Ä,,Jlh"'-'Fi''-')(b)  +  dÄs,h''-^F^''-^)(Jb), 


wo  die  Konstante  C  von   &  und  Je  unabhängig    ist.     (Vgl.  Markoff, 
Differenzenrechnung,  S.  131 — 133.) 
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Viertes  Kapitel. 
Anwendungen  der  Eulerschen  Formel. 

41.    Entwickluugf  von  — j —  in  eine  Reihe.     Setzt  man  in 

e-1 

der  Formel  (8)  (Nr.  37) 

m  =  e^% 

so  wird 

1 


ei^dx=^^-—^,     /"^(l)  =  ^"»e*,     /■H(0)  =  ^'»,     pK){x)  =  t^^e''' 


ß-   . 

0 

und  folglicli 

1  =  '*-=^  +  A(e'-  1)  +Ä,tie*-  1)  +  ■  • . 


^  - .    Nacli  der  Formel  (12)  (Nr.  35)  ist  das  Restglied 

Bei  wachsendem  Ic  bleibt 
endlich,  weil 

oder  1 

1<M,< V- 


22^-1 

ist.     Daraus  folgt,  daß  lim  Jf^.  =  1  ist. 

Es  ist  also  für  |^|<2;t   das  Restglied    dem    absoluten  Betrage 
nach  beliebig  klein  für  hinreichend  große  Werte  von  Je. 

Folglich  besteht  für  |  i^  |  <i27t  die  konvergente  Reihe 

1  =  ^  +  A  (e'  -  1)  +  AKe'  -  1)  +  ÄJ%e*  -  1)  +  •  •  • 
oder 
(1)  -^  =  |+A  +  A^  +  A^'  +  --- 

42.    Entwicklung' '  von    ctg  x   und   tg  x    in   eine  Reihe. 

In  der  Formel  (8)  (Nr.  37)  setzen  wir 

„    '        .  ,  f(x)  =  cos  (tx). 

Dann  wird  ^ 


1 


/■ 


cos  {tx)dx  =  '^^     f'(t)  =  -  t  sin  (tx),     f"'{t)  =  +  ^^  g^^  ^^^^^ 
.  .  .,  fm(x)  =  (-  If-H^'^  cos  (tx), 
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und  folglich  ist 

1  =  ^-7 — |-  Ä^  (cos  t  —  1)  —  A^t  sin  t  4-  AJ,^  sin  < 

+  A,,-.,{-  l)*-1^2*-8  gui  ^  ^_  ^2^(_  \)k-lt2k  cos  (0^). 

Da  aber 

^         ./        ,       ^v        1  +  cos  *  9  i 

1  —  J-i  (cos  ^  —  1)  =  —^2 —  ^  ^^^  "2  ' 

so  findet  man  nach  der  Division  durch 

Bva.t  =  z  sm  —  cos  — 
die  Formel 

4-  (-  1)^-^.42.-2^^^-«+  (-  l)*-^^2.i^*^^^- 

Ebenso,  wie  in  Nr.  41,  kann  man  beweisen,  daß   das  Restglied 
für  i^l  <  2:7r  und  für  hinreichend  große  Je  beliebig  klein  wird. 

Wenn  man  x  statt  -^  einsetzt,  so  entsteht  die  Entwicklung 

(2)  ctg  :r  =  -^  -  2^A^x  +  2^A^x^  -  2^A^x'>  +  ■■■, 

die  für  \x\<.n  konvergiert. 

Mit  Hilfe  der  Relation 

tg  a;  =  ctg  x  —  2  ctg  2x 
findet  man  die  Reihe 

(3)  iq,x==2^i2^~l)A^x-2\2^-r)A^x^-^2'\2^-l)A^x'' , 

w.enn  \x\  <.  -^  ist. 

43.  Stirlingsche  Reihe.    Wir  wenden  die  Formel  (18)  (Nr.  40) 
an  und  setzen 

F{x)  =  \o^x,     a  =  l,    h  =  l,     h  =  X. 
Man  findet  dann 

I  logxdx  =  xlog  X  —  X,    F'(x)  =  x~\    F"(x)  =  —  x~^,-  ■  ■ 
. . .,  jPM(a;)  =  (-  1)"»-!  1 . 2  •  3  . .  .  (m  -  l)x-"' 
und  folglich 

log  [1  •  2  •  3  . . .  (a;  -  1)]  =  C  +  Ä-  log  a;  -  ic  +  ^1  log  ^ 

+  ^a;-i  +  ^,l-2a;-3  +  -..  +  42,_2l-2-3--.(2Ä-4)ir-(»*-») 

+  6^8*1-2. 3.-. (2Ä;-2)a;-(^*-^>. 
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Indem    man    zu    beiden    Seiten    log  x    addiert    und    A^    durch 
seinen  Wert  (—  — j  ersetzt,  so  erhält  man 

log  (1  •  2 . 3  . . .  :r)  =  0  +  (^  4-  -^)  log  aj  -  :?;  4-  ^2  a;-i 

(4)  +  ^J  .  2x-^  +  •  • .  +  ^2 i_2 1  •  2  •  3  . . .  (27c  -  4)ä;-(2*-3) 

4- 0^2i  1  •  2  •  3  . . .  (2Ä  -  2)^-(2*-i). 

Hier  ist  ä;  >  1  vorausgesetzt.     Für  li  =  \  besteht  die  Formel 

(5)  log  (1  -2.3  ...  r^)  =  C4-  (^  +  1)  log  ri;  -  ^  +  0,  J, 


<— 1 


Es  sind  0  und  0^  Größen  zwischen  0  und  1. 
Es    bleibt    noch    die    Konstante    G  zu   bestimmen.     Zu    diesem 
Zweck  wenden  wir  die  Wallis  sehe  Formel  an: 


71 


=  lim  ^{n)     für     lim  w  =  oo, 


wo 

,  ,   s         2-  2     4-4     6 • 6  2%- 2w 

^^        1-3     3-5     5-7         (2w-l)(2n  +  l) 

ist.     Mittels  einfacher  Umformungen  erhalten  wir 

^   ^        Ll-3-5...(2w-l)J     2w4-l 

_  r2-2-4-4-6-6...2n  •  2w"|2         1 
~  L  l-2-3-4...(2n-l)2«  J  "  2w  +  l 

_     (2.4-6...2W)*  1  2*"  •(1-2 -3...%)* 


(1-2-3-4...2W)*     2w  +  l        (l-2-3...2n)*-2w     .    t  JL 

_  „  ,  "'"2» 

Daraus  folgt 

log  ^(w)  =  4w  log  2  +  4  log  (1  •  2 . 3  . . .  w) 

-  2  log  (1  -2.3  ...  2»*)  -  log  (2w)  -  log  (l  4-  ^) 

oder  mit  Benutzung  von  (5) 

log  ^(w)  =  4w  log  2  4-  4C  4-  (4w  4-  2)  log  w  -  4w  +  4017^^  ^ 

-       log2  -  logw  -log(l  +  ^) 

-(4w4-  l)log2-2(7-(4w4-  l)logw4-  4w  -%A^^ 

=  -21og2  4-2C4-40,^,i-log(l4-^)-02^V 
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Wenn  wir  jetzt  n  unendlich  wachsen  lassen,  so  wird 
lim  log  0(n)  =  log  |,      Hm  [40^  A  ^J  =  0, 

limlog(l4-^)=0,        lim[03^i]  =  O 

und  folglich 

log|  =  - 2  log  2  + 2a 

Daraus  findet  man 

(7  =  log}/23r. 
Es  bestehen  also  die  Formeln 

log  (1  •  2  •  3 ... a;)  =  log  yJjt  +  (x -}- ^]  log  X  —  X  -\-  Ä^x-'^ 

(6)  +^^.2^-3  +  ..   +^2*-2l-2-3...(2Ä-4)^-(2*-8) 

l  +  0^2il  •  2 . 3  . . .  (2/c  -  2)a;-(2*-i) 

und 

(7)  log  (1 . 2  •  3 . . . a;)  =  log  yJn  +(x  +  ^\ogx-x+  Q^A^x-\ 

Da  0  <  01  <  1   und  A^=  -z  ist,  so  folgt  aus  (7) 

log  }/2jc  +  (^  +  I)  log  a;  -  2;  <  log  (1 . 2 . 3  . . .  ic)  <  log  Y2% 

+  (^  +  {)log^-a:  +  ^ 
oder 


(8)         y^nx  '  oif  '  e-*  <  1  •  2  •  3  . . .  ic  <  Y¥Jix  •  x^  •  e-"  ■  e^^''. 
Für  große  Werte  von  x  kann  die  Näherungsformel 


(9)  l'2-3  ...x=y27tx-x''-e-^ 

gebraucht  werden. 

44.  Untersuchung  des  Bestgliedes  der  Stirlingschen 
Reihe.  Die  Gleichung  (6)  ist  die  Stirlingsche  Reihe,  deren  Kon- 
vergenz wir  jetzt  untersuchen  wollen.  Zu  diesem  Zwecke  muß  man 
sehen,  ob  das  Restglied  bei  unendlich  wachsendem  k  unendlich 
klein  wird. 

Nach  der  Formel  (12)  (Nr.  35)  ist 

0^2Ä-l-2-3...(2Ä;-2)a;-(2*-i) 
=  0.(-l)-^2..(,i^y.,|^...^^[l  +  ^  +  ^  +  ^,+  ...} 
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In  Nr.  41   ist  bewiesen,   daß  bei  wachsendem  k  sich  die  Reihe 

1 

der  Grenze  Eins  nähert.    Für  jeden  Wert  von  x  wird  aber  das  Produkt 

1         2         3  2k-2 

■2'7tx     2nx     27tx  2%x 

beliebig  groß  für  hinreichend  große  Werte  Ton  yfc. 

Daher  ist  die  Stirlingsche  Reihe  divergent.  Für  große  Werte 
von  X  und  kleine  Werte  von  k  ist  das  Restglied  dem  absoluten 
Betrage  nach  sehr  klein.  In  diesem  Falle  gibt  die  Stirlingsche 
Reihe  eine  gute  Annäherung  für 

log  (1 .  2  •  3  . . .  :r). 

Z.  B.  für  X  =  100  und  ^  =  3  ist  der  Fehler  <  -^-  (Markoff, 
Dijöferenzenrechnung,  S.  137.) 

Für  einen  bestimmten  Wert  von  x  ist  der  Fehler  dem  absoluten 
Betrage  nach  kleiner  als 

1^2,  |.l-2-3...(2;b-2)ic-(2*-i). 

Die  beste  Annäherung  bekommt  man  für  einen  solchen  Wert 
von  Je,  für  den  der  Ausdruck 

—  9        /    1    Y       2  3  2k-2     r-i    ,    J^    ,    J^    ,  "1 

^''~'^^'\2^x)  '  27tx'  2^x"'    27CX    'L    "^2^*        32*"^         J 

den  kleinsten  Wert  annimmt,  weil  der  absolute  Betrag  des  Restgliedes 
kleiner  als  u^  ist. 
Für  k^6  ist 

-i- +  — 4- —  +  —  +  .•.<■ -^. 

22t   ~  ^2k  ^  ^2k   ^  ^2k   ~  ^  2047 

Deswegen  kann  man  angenähert  setzen 

u   =2x    (   ^   Y      ^ - ^^"^ 

*  \2  7tx/      2%x     2Ttx  2  7tx 

Da  aber  angenähert 

—  9        {    ^    Y       ^  2lc-2     2Jc-l      2k 

Uk+i  —  ^X-  \j^)  ■  2^  ^^^  '  ~2^^  "2nx' 

SO  ist 


%+i  _  2k-l      2k    ^  (  k  \2 
u,.  2itx 


2k         fky- 
\%x        \nx) 
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Die  absoluten  Beträge  der  Glieder  der  Stirlingschen  Reihe 
nehmen  also  ab,  solange  k  <  :jtx  bleibt.  Die  beste  Annäherung  erfolgt 
also,  wenn  h  nahe  zu  tcx  -}-  1  liegt.* 


45.  Augenäherte  Berechnung  der  Summe 

für  grofse  Werte  von  x.     In  der  Formel  (18)  (Nr.  40)  setzen  wir 
F(x)  =  -}     a  =  l,    h  =  l,     h  =  X. 

Dann  finden  wir 

äx 


fi 


=  log  X,    F('»)(a;)  =  (-  1)"»1  •  2  •  3  . . .  wa^-^^+i) 
und  folglich 

^4l-2-3^-* ^2i_2l-2-3...(2Z;-3)a;-(2*-2) 

-  dÄ2kl  -2  ■  ^  . . .  {2k  -  l)x-^K 


(10) 


Die  Konstante  C  läßt  sich  nicht  durch  bekannte  Konstanten  aus- 
drücken, sie  läßt  sich  nur  angenähert  berechnen. 

Wenn  man  x  =  10,  k  =  1  setzt  und  das  Restglied  vernachlässigt, 
so  findet  man 

C  =  0,57  721     56649     01532... 

Das  ist  die  Eulersche  Konstante. 

Ebenso  wie  in  Nr.  43  kann  man  zeigen,  daß  die  Reihe  (10)  für 
jeden  Wert  von  x  divergiert.  Deswegen  darf  zur  Berechnung  der 
betrachteten  Summe  nur  eine  geringe  Anzahl  der  Glieder  der  Reihe 
genommen  werden. 

In  die  Untersuchung,  wie  groß  diese  Anzahl  sein  kann,  werden 
wir  uns  nicht  einlassen. 


*  Vgl.  Encykl.  der  Math.  Wiss.,  Bd.I,  p.  103. 
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Differenzengleichnngeii. 


Erstes  Kapitel. 

Allgemeines  über  Differenzengleicliuiigeii. 

46.   Zwei    Formen    der    Differeuzengleichuugeu.      Eine 
Gleicliiing  von  der  Form 

(1)  ^(^,  y.,  ^y,,  ^\, ■■■,  A-2/J  =  0 

heißt  Differenzengleichung.     Hier  ist  0  eine  gegebene  Funktion, 
y    eine  unbekannte  Funktion. 

Bei  Anwendung  der  Formeln 

^\  =  y,+,  -  ^y,+t  +  y., 


^my       ^y  —my,  _|_    W(W  ^)y  1_    (_    1)«  ^^ 

^x        i^x-\-m  ^x-j-m — 1    '  1  •  2       ^x-f-m — 2  '     V  /    ^a; 

nimmt  die  Gleichung  (1)  die  andre  Form  an 

(2)  t(x,  y^,  2/,+i,  y,_^„  . . . ,  y,^J  =  0. 

Ist  die  Gleichung  in  der  Form  (2)  gegeben,  so  kann  man  sie  in 
die  Form  (1)  überführen.     Es  genügt,  die  Formeln 

"x+l  "a;    '"         "x^ 


'x+: 


2/.+  2AI/  +AV, 


,         .        ,   m(m—l)  J.O       ,  lAm 

yx+m-yx+^^yx+-\T2-^^%+--+^"'yx 

anzuwenden. 

Von  diesen  beiden  Formen  (1)  und  (2)  ist  (2)  die  bequemere. 
Es  ist  einfacher,  eine  gegebene  Funktion  statt  y^  in  die  Gleichung  (2), 
als  in  (1)  einzusetzen. 
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Wenn   bei   der  Transformation  (1)    in    (2)    z.  B.    die    Größen  y 
und  y^ti  wegfallen,   so  wird  die  Gleichung  (2)   dadurch  vereinfacht, 
daß  man  ^/ar+s'^^x  setzt. 

Es  sei  die  Gleichung 

gegeben.     Mit  Anwendung  der  Formeln 

^\  =  y.+s  -  ^y.+,  +  ^y^+i  -  y., 
^%-yx+2-^y.+i  +  y.> 

"x         "x-f-l         "x 

geht  diese  Gleichung  in 

2/x+3-2y,+2=0 

über.     Setzt  man  2/3:4-2'^  ■^x'  ^^  entsteht  die  einfachere  Gleichung 

Aus  diesen  Gründen  werden  wir  voraussetzen,  daß  die  gegebene 
Gleichung  die  Form  (2)  hat  und  daß  in  ihr  y^  enthalten  ist.  Wenn 
außer  y^  noch  y^,^,  aber  kein  y  mit  größerem  Index  vorkommt,  so 
werden  wir  (2)  eine  Gleichung  m*®'  Ordnung  nennen.  Die  eben 
betrachtete  Gleichung 

A\-{-A\-Ay^-y^  =  0 
ist  also  eine  Gleichung  1*°'  Ordnung. 

47.  Partikuläre  und  allgemeine  Lösungen.  Wir  werden 
voraussetzen,  daß  die  vorgelegte  Differenzengleichung  »i*®'  Ordnung 
die  Form  hat 

^^)  Vx+n^  fy"^'  Vxi  yx+v  yx+i'---f  yx+n-iJ- 

Die  Funktion  f  soll  für  alle  betrachteten  Werte  von  x,  y^,  y^t^, 
y   .c'--,V   ,       ,  einen  endlichen  und  bestimmten  Wert  erhalten. 

Jeder  Ausdruck  y^=  (p{x),  der  dieser  Gleichung  genügt,  heißt 
partikuläre  Lösung  von  (3). 

Es  gibt  unendlich  viele  partikuläre  Lösungen.  Um  dies  zu  zeigen, 
wollen  wir  den  Größen  y^,  y^,  y^,..,  y„_i  willkürliche  Werte 
beilegen. 

Aus  (3)  ergibt  sich  für  x  =  0,  1,  2,  3, . . .  das  System  der 
Gleichungen 

Seliwanoff,  Differenzenrechiiang.  6 
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2/„+i  =  /■(!;   2/1;   2/2>   ys,---,   y„), 
2/„+2==/'(2,   y^,   yä7   yi,-'-f   y„+i), 
yn+z^f{^>   2/3.   ^4;   2/5;---;   2/^+2)^ 


(4) 


aus  welchem  y^,  y^_^^,  y^j^^,  y^j^^,-"   sukzessive  berechnet  werden 

können.     Es   hat   also   y^   einen    ganz    bestimmten   Wert   für   ganze 

positive  X,  sobald  y^,  y^,  y^,  '■-,  2/„_i  gegeben  sind.     Man  kann  also 

schreiben 

(5)  y^-F{x,  y^,  y^,  y^,...,  y^_^. 

Die  Wahl  der  Anfangswerte  y^,  y^,  y^,.-.,  y^-x  ist  nur  einer 
Bedingung  unterworfen:  es  sollen  alle  Gleichungen  des  Systems  (4) 
einen  Sinn  haben. 

Es  ist  möglich,  daß  die  Funktion 

f{^>  yx'  yx+v  yx+if-f  yx+n-i) 

einen  bestimmten  Wert   erhält   für   x^a   und   für   beliebige  Werte 

von  y^,  y^^^,  2/^+2;  •••;  Vx+n-r 
Dann  besteht  das  System 


(6) 


Daraus  folgt  die  Relation 

(7)  yx=^{^>  ya'  ya+v  ya+v  ■■■>  ya+n-t) 

für  x^a.  Hier  ist  die  Differenz  x  —  a  gleich  einer  ganzen  posi- 
tiven Zahl. 

In  den  Formeln  (5)  und  (7)  sind  y^,  y^,  y^,  ...,  y^_^  oder 
yay  ya+v  ya+2>  •••'  ^«+«-1   willkürliche  Konstanten. 

Jede  Lösung  der  Grleichung  (3)  von  der  Form 

(8)  yx-t{x,  C„  C„  C„  ...,  (7J, 

die  n  willkürliche  Konstanten  enthält,  heißt  allgemeine  Lösung, 
wenn  es  möglich  ist,  die  Grleichung  (8)  in  die  Form  (5)  oder  (7) 
zu  bringen. 


2/«+„  =/■(«, 

ya> 

ya+v 

ya+v  • 

•''      ya+n  —  lj' 

2/a+«  +  l=/"(«^+l^ 

ya+v 

ya+v 

ya+z^  • 

"}       ya  +  n)> 

2/.+„  +  2  =/"(«+ 2, 

ya  +  2f 

ya+zy 

ya+v  • 

■•>       2/a  +  M  +  l)> 
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Zu  diesem  Zweck  muß  man  aus  dem  System 
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/A  =  ^(l,       C„C„...,C„) 


(9)  oder  (10) 


die  Konstanten  C^,  C^,  .  . .,  C^  durch  y^,  y^,  y^,  . .  .,  «/„_i  oder  durch 

2/a>  2/a+i.  ^«+2.  •'•>  ya+n-i  ausdrücken. 

Man  nennt  die  willkürlichen  Konstanten  C^,  C^,  .  .  .,  C„  wesent- 
lich, wenn  das  System  (9)  oder  (10)  in  bezug  auf  diese  Konstanten 
lösbar  ist. 

Allgemeine  Lösung  von  (3)  ist  also  eine  Lösung,  die  n 
wesentliche  Konstanten  enthält. 

Aus  einer  allgemeinen  Lösung  lassen  sich  viele  andere  ableiten. 
Es  genügt  in  (8) 

C^.=  «iA+«2;,^2  +  --+«,,C„ 
(h=l,  2,  ...,  n) 

zu  setzen.  Die  Koeffizienten  a^^  sind  nur  einer  Bedingung  unter- 
worfen, daß  die  Determinante  1  a^^^  1  (/b  =  1,  2,  . . . ,  n;  7i  =  1,  2,  . . . ,  w) 
nicht  verschwindet.  In  diesem  Falle  lassen  sich  c^,  Cg,  .  .  . ,  c„  durch 
Cj,  Ca,  .  .  . ,  C„  ausdrücken. 

48.    Eigenschaften  der  linearen  Di£ferenzengleichungen. 

Eine  Gleichung  von  der  Form 

(11)  y.+n+Piyx+n-i+p^y^+n-2  +  --+Pn-iy.+i-^Pny.=  Q. 

heißt  lineare  Gleichung.  Die  Koeffizienten  p^,  p^,..-,  i>„_i>  i>„ 
und  das  letzte  Glied  Q^  sind  gegebene  Funktionen  von  x.  Wenn 
wir  Q^  durch  Null  ersetzen,  so  entsteht  die  entsprechende 
homogene  Gleichung 

(12)  y,+„ + Ä2/.+„_i + P2y,+n-2 + •  •  +  p»-iy.+i + Pny. = ^^ 

Die  Gleichung  (11)  heißt  eine  vollständige,  wenn  Q^  nicht 
gleich  Null  ist. 

Es  sei  Yx  eine  partikuläre  Lösung  von  (11).  Wenn  man  in 
die  linke  Seite  dieser  Gleichung 

(13)  y.-Y.+  ^. 

einsetzt,  so  findet  man 

5* 
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+  ^.+n    -^Pl^+n-l  +  P2^.+n-2-^-'-  +  Pn-l^.+l-^Pn^.- 

Die  erste  Zeile  ist  nacli  der  Voraussetzung  gleich  Q  .  Daraus 
folgt,  daß  der  Ausdruck  (13)  der  vollständigen  Gleichung  genügt, 
wenn  ß^  eine  Lösung  der  entsprechenden  homogenen  Gleichung  (12)  ist. 

allgemeine  Lösung  von  (12)  ist,  dann  ist 

y.=  ^a.-^ "pi^,  ^i>  c's, ...,  a) 

allgemeine  Lösung  der  vollständigen  Gleichung  (11). 

Lineare  homogene  Gleichungen  besitzen  folgende  bemerkenswerte 
Eigenschaften : 

1.  Wenn   0^   eine  Lösung   von    (12)    ist,    so    ist    auch  C0 
eine  Lösung.     Dies  folgt  aus  der  Relation 

00  ,    -\-p.Gz  j_     ^■{■PiCz  ,      „H \- p     .Gz  ..-{-pCs 

=  G(z  ,    +  »1^  ,      ,  +  «90  ,      „H \-pz\ 

\  x-\-n    '    -'^1    x-\-n  —  1    '    -^'5    x-f-n  —  2    '  '    -t^n   x) 

Hier  ist  G  eine  willkürliche  Konstante. 

2.  Wenn  ^^^  und  0^^)  Lösungen  von  (12)  sind,  so  ist  auch 
^f(^)  +  0^^)  eine  Lösung,  weil 

^a.+  ^a«  +  i'i(4';«-x  +  0f|_,)  +i^,(0i!{._,  +  0f) ,_,)  +  •  •  • 

'     V  x-\-n    '    J^l    x-\-n  —  1    '    -t^2    x-\-n — 2    '  '    J^n    x   } 

ist. 

Aus  diesen  Sätzen  folgt  der  wichtige  Satz: 

3.  Sind    ^W,    ^(2),    ....   ^W    verschiedene    Lösungen     der 

X    '         x    '  '         x  O 

linearen  homogenen  Gleichung,  so  ist  auch 

(14)  2,^=0,4i)  +  C,.f  4----+Ce) 

eine  Lösung.  Hier  sind  0^,  Cg,  . .  .,  Gn  willkürliche  Konstanten. 
Wenn  dieselben  wesentlich  verschieden  sind,  dann  ist  der  Ausdruck 
(14)  allgemeine  Lösung  von  (12). 

In  Nr.  47  haben  wir  gesehen,  daß  eine  Differenzengleichung  eine 
Funktion  definiert,  wenn  die  Anfangswerte  gegeben  sind.  Um  n 
partikuläre  Lösungen  zu  bilden,  genügt  es,  n  Systeme  von  Anfangs- 
werten anzugeben 
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(15) 


*1  ' 


*0   '         1   f 


2     > 


71  —  1' 

n — 1 ' 


«(»)         «(w)         *(«) 

*0    '  1    '       *2    ?  ' 


n  —  1 


Wenn  die  Determinante,  die  aus  diesen  Elementen  gebildet  ist, 
nicht  verschwindet,  dann  sind  im  Ausdruck  (14)  die  Konstanten 
wesentlich  verschieden  und  die  Formel  (14)  ist  wirklich  allgemeine 
Lösung  der  Gleichung  (12). 

Wir  bekommen  also  den  Satz: 

Die  allgemeine  Lösung  einer  linearen  homogenen 
Gleichung  ist  eine  lineare  homogene  Funktion  der  will- 
kürlichen Konstanten. 

Es  war  bei  dem  Beweise  dieses  Satzes  vorausgesetzt,  daß  die 
Koeffizienten  der  gegebenen  homogenen  Gleichung  einen  Sinn  haben 
für  ganz3  positive  x. 

Wenn  aber  wir  nur  behaupten  können,  daß  Pi,  P2,  ■  •  ■  >  P„  end- 
liche und  bestimmte  Werte  haben  für  x'^a,  dann  ist  das  System  (15) 
durch  das  folgende 


(16) 


«(1) 


a-\-n  —  1" 


a    f  a-\-l'  '         a'-\-n — 1 


ZU  ersetzen. 

Nach    diesen    allgemeinen    Betrachtungen    gehen    wir    jetzt    zu 
speziellen  Gleichungen  über. 

Zweites  Kapitel, 
Lineare  Düferenzengleichungen  erster  Ordnung. 

49.    Homogene  G-leichung.     Um  die  allgemeine  Lösung  der 
Gleichung 

zu   bilden,  geben    wir    der  Größe    i/o    einen    beliebigen    Wert.     Aus 
(1)  folgt  dann 

Vi^^ -^2/2  =  -A-^i-^yo)  ■  •  • 


y*=  AA^  ...  Ä^_^tj^. 
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Weil  2/0  willkürlich  ist,  kann  man  schreiben 

wo  C  die  willkürliche  Konstante  ist.     Das  ist  die   gesuchte  Lösung. 

Wir  werden  voraussetzen,  daß  in  (2)  keine  der  Größen  Äq,  ä^, 
A^,  ...,  Ax—i  gleich  Null  ist.  Sonst  hätten  wir  die  evidente 
Lösung  y^=  0,  die  kein  Interesse  bietet. 

Wenn  für  alle  x^a  die  Funktion  Äx  einen  endlichen,  be- 
stimmten und  von  Null  verschiedenen  Wert  hat,  dann  ist  die 
allgemeine  Lösung  unserer  Gleichung 

(3)  y^^C-Aa-Äa+l-Äa+S   .  .  .  Ä^-l. 

50.  Vollständige  Gleichung.     Um  die  Gleichung 

(4)  y.+i-Ay.=  ^. 

aufzulösen,  setzen  wir 

(5)  2/;,  =  «**  •  ^x, 

wo  Ux  und  Vx  neue  unbekannte  Funktionen  sind;  eine  derselben 
wählen  wir  so,  daß  die  Gleichung  für  die  andere  Funktion  möglich 
einfach  sei. 

Vx+l  =  Vx+£^Vx 

ist,  so  nimmt  die  Gleichung 

Ux+l  •  Vx+1  —  AxÜxVx  =  IBx 

die  Form  an 

W*+l(t>«  +  ^V^  —  AxUxVx  =  Bx 

oder 

(6)  (Wx+l  —  AxUx)  Vx  +  Ux+1  ■Avx  =  Bx. 

Wenn  Ux  eine  partikuläre  Lösung  der  entsprechenden  homogenen 
Gleichung  (1)  ist,  dann  erhalten  wir  die  Gleichung 

Ux+l''^Vx  =  Bx, 

die  Vx  bestimmt.     Es  ist 
und  folglich 

Aus  der  Betrachtung  in  Nr.  49  folgt,  daß  man 

Ux  =  AqAj^A^  . . .  Ax—i     oder    Ux  =  AaAa+i  .  .  .  Ax-i 
setzen  kann.     Li  diesen  partikulären  Lösungen  ist  Uq  =  1  oder  w«  =  1. 
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51.    Beispiele.     1.  Es  sei  die  Gleichung 

gegeben.     Hier  ist  ^a,  =  3,  Bx  =  a", 
Die  Formel 

hat  verschiedene  Werte,  je  nachdem  a  ==  3  oder  a  ^  3. 
Wenn  a  von  3  verschieden  ist,  dann  ist 

(-Y 


und  folglich 


2/  -C.3*  + 


a-* 


a  — : 
Wenn  aber  a  =  3  ist,  so  ist 


^i^-i^i^^  +  c 


und  ^ 

t/,=  C.3-  +  |^.3^ 

2.  Wenn  wir  die  Gleichung 

3x4-1 « 

^*+i~^  3x  +  7^*""  (3a;  +  4)(3a;  +  7) 
nehmen,  so  ist 

A    =       3a;  +  l        ^  _  x 

*  3x4-7'  *       (3a;  +  4)(3a;+7)' 

^*  ~  ^~    -^    ■  Y  ■  lö  '  13  ■  16  ■  ■  ■  Sx-j-l  *  3a; 4-4 
1-4 


-  (- 1) 


(3a;-|-l)(3x  +  4) 
1 


*x-i-l 


(-l)-a;, 


2;(-  1)-  •  o;  =  2:^  .  A  ^^  =  -  ia;(-  1)--  l2;(-  1)- 

Folglich  ist  die  allgemeine  Lösung  dieser  Gleichung 

^«  ""  ^~  ^'    '  (3a;-f  l)(3a;-{-4)  '  ^  +  (2"^  "~  4  j  "  (3x+l[)(3x  +  4)" 
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3.  Als  letztes  Beispiel  wollen  wir  noch  die  Gleichung 
_  2a;-f  7     3a;  4-1       _       1 
yx+i       2a; +  1  *  3a;  +  4^a;~  3a; +  4 
auflösen.     Es  ist 

_  1    1    11        2a;  +  5 
^^~  1  '  3  'T"  2^31l- 

1.1    J-         3a;-2  _  (2x+l)(2x-\-S){2x-\-5) 
.      4  ■  7  *  10  ■■■  3a;+l  ~  15(3a;4-l)  ' 

^x    ^  15 

M^+1  ~  (2a;  +  3)(2a;  +  5)(2a;+y)' 

^(2a;  +  3)(2a;-f  5)(2a;+7)  ^  "sX/ 


(^+l)(^+t)(^+l) 

1    /nr J^  J^ L.  r' 


''(^  +  l)(^  +  l)  4(2:^  +  3)(2a;  +  5) 

_  (2a;+l)(2a;  +  3)(2a;  +  5)    ^_1     2x-\-l 
^*  3a;  +  l  4  '  3a;+l' 

(Vgl.  Markoff,  A.A.  Dijfferenzenrechnung,  S.  156  —  158.) 
Jetzt  wollen  wir  die  Anwendungen  der  Gleichungen  erster 
Ordnung  zeigen.  Die  Koeffizienten  verschiedener  Entwicklungen 
genügen  Differenzengleichungen.  Die  Auflösung  derselben  führt 
zur  Bestimmung  der  Koeffizienten.  Als  Beispiel  nehmen  wir  die 
binomische  Entwicklung  und  die  Entwicklung  von  cos  xt  nach  Potenzen 
von  cos  t. 

52.  Die  binomische  Entwicklung.  Das  Produkt  von  x  Faktoren 

hat  die  Form  (1  +  0  (1  +  0  •  •  •  (1  +  0 

(1  +  ty=  l-{-A,t  +  BJ^+  (7,^3+  ...  +  H,t^-^  +  KJ-. 
Der  letzte  Koeffizient  ist  gleich  Eins.     Die  Funktion 
Äx  hat  einen  Sinn  für  x'^1, 

•^x      }}  „  ,)  )i      X  ^  J, 

^x      ,f  „  „  },      X  ^  O, 


Die  Anfangswerte  sind 

Nach  Multiplikation  mit  1  +  ^  findet  man 

(1  +  ty+i  =  1  +  (1  +  Äx)t  +  (Äx + B^y + {Bx  +  Gxy + 

■■■  +  {Hx+Kx)t^+Kxt-^\ 
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Anderseits  ist 

(1  +  ^)x+i  _  1  4.  A,+J  +  B,+it'  +  C,+,t'  +  ■  •  • 

•  •  •  +  H,+,t—-'-\-K,+,t-+  L,+J-+K 

Daraus  folgt 

A^+i=  Ä^+  1, 

Bx+l  =  J?a,  +  Äx, 
Cx+1  =  Cx  -\-  JBx, 

Die  Auflösimg  dieser  Gleicilungen  gibt 

Ax=  21    =  X  -\-  c;     es  ist  c  =  0,  weil  J.^  =  1  ist. 

Bx=2]Äx=Ux  =  '^^f^  +  c. 
Aus  der  Bedingung  B^^  1  folgt  c  =  0  und 

-^*~      1-2      ■ 

Die  willkürliche  Konstante  c  ist  =  0,  weil  Cg  ==  1  ist. 
In  dieser  Weise  erhalten  wir  die  binomische  Entwicklung 

xix-l)ioi:-2)...[x-{x-l)-\ 
"^  1-2-3...X 

Der  letzte  Koeffizient  ist  gleich  Eins,  wie  es  auch  sein  sollte. 

53.  Entwicklung  von  cos  x  t  nach  Potenzen  von  cos  t. 

Nach  der  Moivreschen  Formel  ist 

(cos  ^  +  i  sin  ^)^  =  cos  xt  -\-  i  sin  xt. 

Der  Vergleich  der  reellen  Teile  beider  Seiten   dieser  Gleichung 
ergibt 

cos  xt  =  cos't  —  ^  f  ~     sin*^  •  cos*— ^^ 

,    x(x—l)(x—2)(x  —  i)     ...  ^     ., 

+  -!^ ^^^.s.^ ^^  ^ '  cos*-^^ 

Wenn  man  die  geraden  Potenzen  von  sin  t  durch  cos  t  ausdrückt, 
so  findet  man  eine  Gleichung  von  der  Form 

cos  xt  =  Äx  COB't  +  Bx  COS^-^t  +  Ca:C0S*-*^  +  •'♦ 
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Es  sollen  die  Koeffizienten  Äx,  JBx,  Cx,  ...  bestimmt  werden. 
Der  Kürze  halber  wollen  wir   cos  t  mit  u  bezeichnen.      In  der 
Entwicklung 

cos  xt  =  ÄxW'  -f-  BxU^-^  +  CxW"-^  4-  DxU^-^  +  •  •  • 

hat  die  Funktion 

Ax  einen  Sinn  für  a:;  >  1, 

J>x        ff  j7  )f      OC  ^  Zf 

^x       ),         }}        «     ^  ^  4, 

■Lfx  7J  };  ,f         OC  ^  O, 


Aus    der  Formel   cos  t  =  A^u   folgt,    daß  A^=  1    ist.      Um    die 
übrigen  Anfangswerte  zu  bestimmen,  genügt  es,  in  den  Formeln 

cos  2t  =  A^u^  +  B^, 

cos  4:t  =  A^u^  +  JB^u^  +  O4, 

cos  6^  =  Aqu"^  +  BqU^  +  OßW^  4-  A 


^  =  —  oder  w  =  0  zu  setzen.     Man  findet 

B,=  -l,     C,=  l,    J>6=-1,--- 
Aus  der  Relation 

cos  (x  +  1)^  +  cos  (x  —  l)t  =  2  cos  t  •  cos  xt 
folgt  die  Identität 

Ax+iu-+^  4-  (5^+1  +  ^,_i)m— 1  +  (Cx+i  +  Bx-x)u''-'  4-  •  •  ■ 

=  2^^w^+i  +  2BxU^-^  4-  2a,M^-3  +  •  • . 
Der  Vergleich  der  Koeffizienten  führt  zu  den  Differenzengleichungen 
Ax+i—2Ax==0, 
Bx+i  —  2Bx  =  —  Ax—i, 
Cx+i  —  2Cx  =  —  Bx—Xf 

-Z/ar-fl —  2Dx=  —  Cx  —  lf 


Diese  Gleichungen  vereinfachen  sich,  wenn  man  neue  Funktionen 
einführt. 

Wenn  wir  setzen 

A-2*-i-a^,    Bx==-2^-'-hx,    Cx=-2—^-Cx,...f 


so  finden  wir 


üx+l —  «a;=  0, 
l>x+X —  ix  =  C^x  —  l} 
Ca;  4-1  ^x  ^^  ^x — 1> 

U/x-{-l        (^x  ^^  Cx — 1, 
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Die  Anfangswerte  sind 

«1=1,     62  =  2,     64=2,     d^=2,... 

Aus  der   ersten  Gleichung  folgt,   daß   a^  konstant   ist;   da   aber 
«j  =  1  ist,  so  ist 

a^=l     und     J.x-2*-^ 

Die  zweite  Gleichung  gibt 

hx=  .21  =  X  -{-  c. 

Weil  h^  =  2  ist,  so  ist  c  =  0, 

Ix^^  X        und     5x=  — 2^-8-ä;. 

Die  weiteren  Rechnungen  sind 
Cx=  2;(a;  -  1)  =  ^ ^h +  ^' 

3-2    ,  . 

C4=]--^  +  C.       C  =  -l, 

(a;-l)(a;-2)        ^ 
C^-  1.2  ■"' 

^^^  ^p-2)(p3)  _  ij  ,  (.-2)(.-8)(.-4)  _  ^  _^  ^^ 

2  =  ^^-6  +  «=,    ^  =  4, 

_  (a;  -  3)  (x  -  4)  (o;  -  5)  (a;  -  6)  _  (a;-5)(a;-6) 

""  1.2-3-4  1-2  "^^ 

^"1.2.3.4        1-2  +  ^'      ^        ^' 

_  (J?-3)(a;-4)(a;-6)(a;-6)       (a!-6)(a;-6) 
^*~  1.2-3 -4  l-,2  ' 

„  _  "^  rix  -  4)  (a;  -  5)  (a;  -  6)  (a;  -  7)  _  {x  -  6)  (a;  -  7)-| 
'*~^L  1-2-3-4  1-2  J 

_  (g - 4)  (x - 5)  (a; - 6)  (x - 7)  {x - 8)  _  (g-6)  (x-7)  (x -8) 

1.2-3-4-5  1-2-3  "*"' 

g        65-43-2        4-3-2  ^^ 

^~l-2.3.4.6        1.2.3  +  ^'      ^~^' 

^  _  (x-4)(x-6)(x-6)(x-7)(x-8)        (x  -  6)  (x  -  7)  (x  -  8) 
'*—  1-2-3-4-5  1-23  ' 

Die    gefundenen  Ausdrücke   für   Cx,    dx,   ex,   fx  lassen    sich    ver- 
einfachen.    Es  ist 
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x{x  —  3) 


c.=  (l-f)(l-|)-l  =  ^ 


2 


d:r 


x{x  —  i)(x  —  5) 


(.-4).[(l-|)(l-|)-l]  =  --fc:,.3      , 

{x  —  6){x  —  6)     r/^        ^\ /i        ^\       -il       x{x  —  5){x  —  &){x  —  l) 


[(i-|)(i-f)-i]  = 


1-2  L\  3/  \  4/  J  1-2-3-4 

{x  —  &)(x—l){x  —  8)     r/^       a;\  /.       a;\       .,1      x(x—&){x—l){x—8){x~9) 


A~'-''><-:7.^f-^>.[(i-l)(i-l)-i]= 


1 • 2-3-4-5 


Die  gesucliteii  Koeffizienten  haben  also  die  Werte 

Tj  9a;_7    a?(a;-4)(a;-5)         ^        „       g    a;(a;  -  5)  (a;  -  6)  (a;  -  7) 

^  ^  _  g^_ji    a;(a;-6)(a;-7)(a;-8)(a;-9) 
*  ■  1 -2. 3-4-5  '  ■ 

P.  L.  Tscliebysclieff  hat  bemerkt,  daß  die  ganze  Funktion 
x^^^  Grades 

/"W  =  ^:?r  =  w*  +  i?i  w^~^  4-  p^u''-^  + H  P^_^u  +  i?^ 

im  Intervall    (—  1,  +  1)    weniger    von   Null    abweicht,     als    jede 
andre  ganze  Funktion  von  derselben  Form: 

cp{u)  ==  Wx  +  g'iW^-^  +  q^u''-^  -{ h  q^_^u  +  q^. 

(Oeuvres  de  P.  L.  Tchebychef,  t.  I,  p.  295—301.  Vgl.  J.  Ber- 
trand, Traite  de  calcul  differentiel  et  de  calcul  integral,  t.  I, 
p.  512  — 521.) 

Die  Abweichung  von  Null  ist  der  größte  absolute  Betrag 
einer  Funktion,     Wenn  z.  B. 

1  F{u)  \^A    für     a  <  M  <  &, 

so  ist  A  die  Abweichung  von  Null  der  Funktion  F(u)  im  Intervall  {a,  h). 
Die  Abweichung  von  Null  der  Funktion 

f(u)  =  -j—j:  cos  (x  arc  cos  u) 

im  Intervall  (—  1,  +  1)  ist  gleich    ^_^  • 

Wir  wollen  voraussetzen,  daß  im  Intervall  (—1,  +1)'  die 
Funktion  (f{u)  weniger  von  Null  abweicht,  als  fiy). 
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Für 


ist 


^  n  2  7t  {x  —  l)it 

'  X              X                               X  * 

^                In                2«  (x  —  l)n  ^ 

U  =  \,  cos—;  COS >•••;       COS  ^^ —)        —  1, 

'                    X                        X  X  ' 


f(u)  =  -^, ^—,       -J^,  .-.,        (_  1)«-1  .  -J—,       (_l)a:_J_. 

Wenn  wir  der  Kürze  wegen  die  Bezeiclinung  einführen 
Uk=coa-^,   (k  =  0,  1,  2,  ..  .,  x), 
so  bestehen  nach  unsrer  Voraussetzung  die  Ungleichheiten 

fM  -  9^M  >  0,     /"(%)  -  9p(wi)  <  0,     /"(wa)  -  (p(u^)  >0,  ... 
In  jedem  der  Intervalle 

(^0)    Wi),       (Wj,    U^),   .  .  .,      (Wx-1,    Wx) 

hat  die  Gleichung  f{u)  —  (p(u)  =  0  mindestens  eine  Wurzel  und  also 
im  Intervall  (—  1 ,  +1)  mindestens  x  verschiedene  Wurzeln.  Dies  ist 
aber  unmöglich,  weil  die  ganze  Funktion 

f{u)  -(p(u)  =  (j>^-  g'i)w^-i  +  (p^  -  <ii)u='-^  + Vijpx-  gx) 

den  Grad  x  —  1  nicht  übersteigen  kann. 

Damit  ist  also  der  Satz  von  Tschebyscheff  bewiesen.  Diesen 
Beweis  verdanke  ich  einer  mündlichen  Mitteilung  von  Herrn 
A.  A.  Markoff. 

(Vgl.  Lacroix,  Traite  du  calcul  differentiel  et  du  calcul  integral. 
Nr.  1052  — 1054.) 
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54.    Partikuläre   Lösungen    der    homogenen   Gleichung. 

Um  die  Lösungen  der  Gleichung  mit  konstanten  Koeffizienten 

(1)  2/.+„  +  i>i2/,+„_i  +  i>22/:,+„_2  +  *-+i>„2/.=  0 

zu  finden,   werden  wir  in  die   linke  Seite  dieser  Gleichung  einsetzen 

(2)  y^=-a'U:„ 

wo  Ux  eine  ganze  Funktion  »w**°  Grades  ist. 
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Mit  Benutzung  der  Formeln 

2/^+3=  <*''+'K+  3Aw^+  3A2m^  + A^M^), 

2/.+,  =  «"+'*  (w.  4- »*  A  «^,  +  ^?^^^  A2w,  +  . . .  +  A«M,) 
findet  man  als  Resultat  dieser  Substitution  die  Identität 

+  a^+1  [wa«-i  +  (w  -  l)i)ia'»-2  +  (w  -  2)^2 «^«-3  +  . .  •  +  i>„_  J  A  W;, 


+  a^+«  •  w(w  -  1)  («*  -  2)  •  ■  •  2  •  1  •  j 


2..  .w 

Zur  Abkürzung  wollen  wir  die  Bezeichnung  einfübren 
(3)  f{z)  =  0-  +  pi^«-i  -f  p^^-^  +---  +  P,. 

Die  gefundene  Identität  nimmt  dann  die  Form  an 


(4) 


=  a-  ■  f{a)  •  w,  +  a-+i .  /-'(a)  •  Am,  +  a-+2  •  /•"(«)  •  ^  + 


Unter  «/^  ist  hier  a'-Ux  zu  verstehen. 

Die  rechte  Seite  von  (4)  verschwindet,  wenn  a  eine  Wurzel  der 
Grleichung  f{ß)  =  0  und  Ux  =  1  ist.  Eine  partikuläre  Lösung  von  (1) 
ist  dann  a". 

Wir  wollen  voraussetzen,  daß  a  eine  )fc- fache  Wurzel  der  Grleichung 
f(/)  =  0  ist.     Die  rechte  Seite  von  (4)  ist  dann 

„.+..^«(„) .  _iV_  +  «.+»+x.^<*+.,(«) .  _^!^ +... 

A"w^ 
'       ^  ^     1 • 2  . .  .n 

Dieser  Ausdruck  ist  gleich  NuU,  sobald  Ux  eine  ganze  Funktion 
vom  (Je  —  1)*®'^  oder  niedrigeren  Grade  ist.     Wir  setzen 
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M«  =  1; 

Uo  =  0,     Aii^=l; 


vk-l. 


Uq  =  0,    Auq  =  0,     AhiQ  =  0,  . . . ,     ^.g      ^^1^^  =  1. 

Dann  ist 

Ux=l,    Ux  =  oo,     Ux  =  x(x— 1),  . . . 

Ux  =  x{x  -  l){x  -2)  ...  (x- Je +  2)', 
wir  haben  also  Je  partikuläre  Lösungen  von  (1) 

a^,     a^'x,     a^ 'x{x— 1),  .  . .,     a""  ■  x{x  —  1)  .  . .  (x —  Je -\- 2). 

In  Nr,  48  war  es  schon  gezeigt,  daß  den  partikulären  Lösungen 
einer  linearen  homogenen  Gleichung  ein  konstanter  Faktor  hinzu- 
gefügt werden  kann.  Deswegen  kann  man  die  gefundenen  partikulären 
Lösungen  durch  die  folgenden  ersetzen 

a^,    xa'^-^    x{x  -  l)a*-2,  .  .  . ,     x{x  -  l)  .  . .  {x -Je -\-  2)a^-*-^ 

Um  diese  Lösungen  zu  bilden,  genügt  es,  die  Funktion  t'  zu 
nehmen  und  (Je  —  1)  sukzessive  Ableitungen  nach  t  zu  bilden  und 
in  den  Resultaten  t  durch  a  zu  ersetzen. 

55.    Allgemeine  Lösung  der  homogenen  Gleichung.    Bei 

der  Auflösung  der  Gleichung 

(1)  !/.+„  +  i^i2/x+.-i+i^2^x+«-2+--+i'«2/.=  0 

werden  wir  zwei  Fälle  unterscheiden:  1.  die  entsprechende  algebraische 
Gleichung 

hat  verschiedene  Wurzeln  und  2.  die  Gleichung  ({s)  =  0  hat  gleiche 
Wurzeln. 

Wenn  wir  im  ersten  Falle  die  Wurzeln  mit 

bezeichnen,  dann  sind  nach  Nr.  54 

partikuläre  Lösungen  von  (1)  und  nach  Nr.  48  genügt 
(5)  2/,  =  C-ia-  +  C^al  +  Cs«^  + .  •  •  +  C««^ 

der  Gleichung  bei  willkürlichen  Werten  der  Konstanten 

^v  ^2»  ^3>  •  •  •>  ^«- 
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Um  zu  beweisen,  daß  der  Ausdruck  (5)  die  allgemeine  Lösung 
ist,  genügt  es  zu  zeigen,  daß  Oj,  O2,  ...,  C„  sich  durch  willkürlich 
gegebene  Werte  «/g,  «z^,  y^,  .  .  .,  y^_i  ausdrücken  lassen. 

Um  das  System 

Vi  =  C^a^      -\-  C^a^     +  Ce^a^      H-  •   •  +  Cnf^n, 


aufzulösen,  werden  wir  die  Gleichungen  des  Systems  mit  unbestimmten 
Faktoren  g_^_^,  ö'„_2>  Q-n—v  •  •  • '  So  multiplizieren  und  dann  addieren. 
Das  Resultat  nimmt  die  Form  an 

(7)  \ 

1=  Cj^(p{a^)  +  C^(pM  +  C^<p{a^  + f-  C„g){an), 

wenn  man  die  Bezeichnung 

(8)  9'(^)  =  2o^-'+2i^''-'  +  --  +  g„_2^  +  ?._i 
einführt. 

Um  Cj  zu  bestimmen,  werden  wir  qQ,  q^,  .  .  .  q.n—\  ^^  wählen,  daß 

wird.     Es  ist  in  diesem  Falle 

^  ^^^  ~  («1  -  «2)  («1  -  »s)  •  •  •  («1  -  ««) 
und  die  Relation  (7)  gibt 

In  derselben  Weise  lassen  sich  C^,  ^3,  ...,  Cn  durch  i/q,  «z^, 
Vi,  ■■•,  Vn-x  ausdrücken. 

Damit  ist  also  bewiesen,  daß  der  Ausdruck  (5)  wirklich  q.11- 
gemeine  Lösung  Ton  (1)  ist. 

Jetzt  wollen  wir  voraussetzen,  daß 

gleiche  Wurzeln  besitzt. 

In  allen  einzelnen  Fällen  ist  die  Beweisführung  dieselbe.  Des- 
wegen genügt  es,  nur  einen  Fall  zu  betrachten.  Wir  werden  voraus- 
setzen, daß 

f{3)  =  {3-  a^f{z  -  a^{z  -a^)  ...  (z-  an) 


(10) 
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ist     und    daß     die    Zahlen     a^,    a^,    a^,    ...,   a„    voneinander    ver- 
schieden sind. 

Da   wir   n   partikuläre   Lösungen   von  (1)    kennen    (Nr.  54),    so 
genügt  der  Ausdruck 

{y.=  ^i«T  +  CgO^api-l-  G,x(x  -  l)a^-2 
I      +  C^a?  +  C,a^  +  •  ■  •+  C  «^ 

l  '         4    4'         55'  '         n    n 

unsrer  Di£ferenzengleichung.     Es  ist  noch  zu  zeigen,  daß  das  System 


\     yo  =  C,l     +6V0  +03-0  4-04-1     +--+C„-1 

in  bezug  auf  C^,  C^,  .  •  .,  C„  lösbar  ist. 

Wenn  man   die   Gleichungen   (10)   mit    q^,    q.^,  ...,   q^_2)  Un—i 
multipliziert  und  zusammenaddiert,  so  findet  man 

j    ^o^„_i  +  ^i^„_2  + •■•  +  3„_2^i  +  g„_i2/o 

1=  Gi<pM  +  G29>'M  +  C39P"K)  +  C'49  W  +  •  •  •  +  Cn(p(a„). 

Hier  ist 

(p{z)  =  q^z»--"  +  q,z»-^  +  . . .  4-  ^^_^^  +  q^_^. 

Um  Ci  zu  bestimmen,  genügt  es,    die  Koeffizienten  qQ,  q^,  .  .  ., 
q„_2f  Qn—i  ^^  ^^  wählen,  daß 
(p(a,)  =  1,     (p'(a^)  =  0,     <p"(a^)  =  0,     (p{a^)  =  0,  .  .  . ,     (jp(a„)  =  0. 
Um  die  Gleichungen 

ffM  =  0>     <pM  =  0,  .  .  . ,     (p{an)=-0 
zu  befriedigen,  genügt  es,  zu  setzen 
(12)  (p{3)  =  {z-  a^{z  -  «5)  ...  (^  -  an)il>{2), 

wo  '^{z)  eine  willkürliche  Funktion  ist. 
Aus  (12)  und  aus  den  Bedingungen 

9)(ai)  =  l,     9)'(aJ  =  0,     ^"(aj  =  0 

folgen    die    Werte    für    '4>{a^,   ri^'ia^),   t^"(«i).     Die    Funktion    vom 
niedrigsten  Grade,  die  allen  diesen  Forderungen  genügt,  ist 

qp(^)  ^{z-  a^{z  -  ag)  ...(z-an)  [^(aj  +  (^  -  ai)t^'(aO 
4.(£Z^V"(%)]. 

Seli-wanoff,  Differenzenrechnnng.  R 
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In  derselben  Weise  lassen  sich  die  Konstanten  C^  und  C^  durcli 
yo>  yii  •-•>  Vn-x  ausdrücken. 

Um  C^  zu  bestimmen,  genügt  es,  zu  setzen 

^^^^  ~  («4  -  «i)'(a4  -  «5)  («4  -  «e)  •  •  •  («4-  ««)" 

Ebenso  lassen  sich  die  übrigen  Konstanten  Cg,  Cg,  ...,  6'„ 
bestimmen. 

Damit  ist  also  bewiesen,  daß  (9)  wirklich  eine  allgemeine 
Lösung  der  vorgelegten  Dijfferenzengleichung  ist. 

Durch  die  Einführung  andrer  willkürlicher  Konstanten  kann 
diese  Lösung  eine  einfachere  Form  annehmen.     Da  der  Ausdruck 

Ol  4-  C^a~'^x  4-  C^aY'^x{x  —  1) 

eine  willkürliche  Funktion   2*®^  Grades   ist,   so  kann   man  ihn  durch 
A^  -\-  A^x  +  A^x^  ersetzen. 

Wir  bekommen  die  folgende  allgemeine  Lösung  unsrer  Gleichung 


(13) 


y,-{A  +  A^^  A^>1  +  c,a-  4-  •  •  •  +  C^al- 


Li  allen  Fällen  gibt  es  also  solche  n  Lösungen  u^^\  u^^\ 
der  gegebenen  homogenen  Gleichung,  daß  man  das  System 


,  m(") 


^n — 1  1    n  —  1    '         2    n  —  1    '  '         n    n  —  1 

in    bezug    auf   C^,    Cg,  ...,    0„   auflösen    kann.     Deswegen    ist    die 
Determinante 


(14) 


n  —  1  n — 1 


«*L"^  1 

n  —  1 


von  Null  verschieden. 

Als  Beispiele  wollen  wir  einige  Gleichungen  mit  ihren  allgemeinen 
Lösungen  angeben. 
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56.    Olelchuugeu   mit   reellen  Koeffizienten.     Wenn   alle 

Koeffizienten  der  Gleichung 

y,+„ + Piy,+„-i  +  p,y,+„.,  +  ■■■+ pjj.  =  ^ 

reell  sind,  so  kann  man  verlangen,  daß  in  der   allgemeinen  Lösung 
keine  imaginären  Größen  auftreten. 

Wir   werden    voraussetzen,    daß    die    entsprechende    algebraische 
Gleichung 

eine    dreifache   komplexe  Wurzel   hat,    die   wir   mit    «i    bezeichnen. 
Nach  einem   Satze  der  Algebra   hat    diese   Gleichung    eine    dreifache 
konjugierte  Wurzel,  die  wir  mit  a^  bezeichnen.    Die  übrigen  Wurzeln 
a-j,  a^,  .  .  .,  a„  soUen  reell  und  voneinander  verschieden  sein. 
Die  allgemeine  Lösung  unsrer  Gleichung  ist 

y.=  (0,+  C,x  +  C,x')al  +  (C,+  C,x  +  C\x')a2 

Wenn 

«1  =  ()  (cos  a  +  *  sin  a) 
gesetzt  wird,  so  ist 

a^—  Q (cos  a  —  «  sin  a), 

a*  =  ()'(cos  XU  -f  i  sin  xa), 

a*  ==  9^  (cos  xcc  —  i  sin  xcc). 
Deswegen  ist 

(Ol  +  C^x  +  Ogo;»)«*  +  (C^+  C^x  +  C,x')al 
=  [(Ol  -\-C^x-\-  C^x^)     +  (0,  +  C,x  +  C,x^)]q'  cos  ax 
+  [(G,  +  C^x  +  C,x')  i   -  {C,  +  C,x  +  Oe^r»)  i]  q-  sin  ax 
Wenn  man  neue  Konstanten  einführt 

Ä  =  Ol  +  Oj,        -^1  =  Oj  +  Cß,        -^8  =  Oj  4"  Oj, 
B  =  (Oi-0,)^      ^1=  (0,-05)*,      J5,  =  (03-OeK 
so  erscheint   die   allgemeine  Lösung   unsrer  Gleichung  in  der  Form 
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Vx^Q^'  M  +  A^  +  A^')  cos  xa  +  {B  +  B^x  +  B^x^)  sin  xu'\ 
+  C,a^,  +  --.  +  C^al, 

die  keine  imaginären  Größen  enthält. 
Beispiele: 

2/^=  2-.  [^  cos  (^^)  +  5sin  (^^)]. 

y^-  Ci  cos  (xj'j  +  C2  8in(a;|)  +  O3COS  (x^^  +  C^4  sm(;r  ^)- 

57.    Vollständige    Gleichung.      Um    die    Lösung    der    voll- 
ständigen Gleichung 

(15)  y,+n  +  Piy:c  +  n-l+P2y.+  n-2-^-+Pny.=   Q. 

ZU  suchen,  werden  wir  voraussetzen,  daß  die  Lösung  der  entsprechenden 
homogenen  Gleichung 

schon  bekannt  ist.     Es  sei   die   allgemeine  Lösung  dieser  Gleichung 

Statt    der    willkürlichen    Konstanten    werden    wir    unbekannte 
Funktionen  von  x  einführen,  und  der  Ausdruck 

(16)  y  =  OWw(i)  +  0(2)^(2)  +..._!-  C^»)u'n) 

\         /  '^X  X         X        '  X         X        '  'XX 

soll  der  vollständigen  Gleichung  genügen.  Eine  Bedingung  haben 
wir  zu  erfüllen  und  n  Funktionen  zur  Verfügung.  Deswegen  können 
wir  n  —  1  neue  Bedingungen  einführen.  Lagrange  wählt  dieselben 
so,  daß  y^,^,  y^_L.^,  •  •  •;  y^-i-n—x  dieselbe  Form  haben  sollen,  als  ob 
die  Cx  Konstante  wären. 

Da  Cx-{-x  =  Oj:  -f  AOx  ist,  so  finden  wir  das  System  von  Gleichungen 

yx+,   =q^^-^iV    +GT'^f+.   +-+G^:'-^i\v 

^x-\-i  X  x-\-2  '        X  x-\-2  '  '        X  x-\-2'  ^ 

y  =OW-mW  J_0(2).m(2)  +...  +  0W-w(1        ,.* 

''«+« — 1  X  x-\-n  —  1    '        X  x-\-n  —  1    '  '        x  x-\-n  —  1 
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Diese  Formeln  gelten,  wenn  die  Funktionen  Cx  den  Bedingungen 
genügen 


(17) 


/.(i) 


Es  ist 

^       =  Cd)  •  w(ii  +  C(2) .  w(2)     4- .  ■ .  +  OW  •  wW 

^x-\-n  X  x-f-n    '        X  x-\-n    '  '        x  x+n 

+  wj).  „ •  AC'(i)  +  w(f)  „  •  AOf  +  . . .  +  ,,w    .  AC(«), 

*       x-j-n  X  x-f-w  X  x-\-n  x    ' 

und  es   kann  nicht   mehr  «   ,      dieselbe  Form   haben,   als   ob   die  Gx 

»'x  +  n  '  ■^ 

konstant  wären. 

Wenn  wir  die  Werte  y^,  y^_^^,  y^^^,  ...,  y^_^^  in  die  gegebene 
vollständige  Gleichung  einsetzen,  so  fallen  die  Glieder  mit  C^\ 
Cf\  .  .  . ,  G^^  weg  und  es  bleibt  die  Gleichung 

(18)         w(i)     •  AC(i)  +  ^t(2)     .  AOW  +  . . .  +  mW    .  aOW  =  ^  . 

V       /  x-|-n  X       '        x  +  w  '''  X+'*  ^  ^ 

Aus  den  m  Gleichungen  (17)  und  (18)  mit  den  w  Unbekannten  AO^^^, 
AC^^),  . . . ,  AO^")  lassen  sich  dieselben  bestimmen,  wenn  die  Determinante 

X       '  '  T.  ' 


F{x) 


^%.     ^.-■-    ^%. 


x  +  2 


mW  , 

x-\-% 


w(i) 

x-\-n 


«f+. 


x-j-w 


von  Null  verschieden  ist. 

Wir  haben  schon  gesehen  (Nr.  55),  daß  die  partikulären  Lösungen 
**x  ^  w^^),  .  .  . ,  M^")  der  homogenen  Gleichung  so  gewählt  werden  können, 
daß  die  Determinante 


<) 

<)  .. 

•    wW 

F{-  1)  = 

Mf) 

lef     .  . 

.     MJ") 

ei 

<^l.i  •  • 

•     <ll 

nicht  gleich  Null  ist. 
Aus  der  Relation 

x+n    '    J^i    a:+«  — 1-^2    *+n  —  2    '  '    ^n  —  1     x+1    '    -^n     x 

folgt,  daß 
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F{x)  = 


X-\-i. 

x-f2 


**«+«— 1 


a;+«  —  1 


mW 

X-\-} 


—  p  u^^^     —  p  w(^) 

-'^«    X  -^n    X 


P  d"^ 

-t    M        X 


weil  man  zu  den  Elementen  der  letzten  Zeile  die  Elemente  der  vorigen 
Zeilen,  resp.  multipliziert  mit  p^^  p^,  •  •  •>  l>„_i>  addieren  kann. 
Es  besteht  also  die  Relation 

Wäre  F(x)  =  0,  so  müßte  auch  F(x  —  1)  ==  0  sein,  weil  p^  nicht 
Null  ist. 

Aus  denselben  Gründen  wäre  auch 

F(x-2)  =  0,    i^(a;-3)  =  0,  ...,    JP(- 1)  =  0, 

was  aber  nicht  der  Fall  ist. 

Es  ist  also  bewiesen,  daß  F(x)  nicht  Null  ist.  Die  Gleichungen 
(17)  und  (18)  sind  lösbar.     Man  findet 

A  C(i)  =  «;W,     A  0(2)  =  «;(2),  . . . ,     A  0(«)  =  «;(«) 

X  X    '  X  X    '  '  X  X    • 

Die  Summation  gibt  Resultate  von  der  Form 

Demnach  ist  die  Lösung  der  vollständigen  Gleichung 

2/.=  C>W+02Mf+---  +  0„M« 

-f  m(1>  W;(1)  +  t((2)  ^(2)  _| [_  u{n)^{,n) 

'  X  X         '  X  X         '  '  X  X 

Es  ist  die  allgemeine  Lösung,  weil  die  Determinante  (14)  von 
Null  verschieden  ist. 

Diese  Methode  von  Lagrange  heißt  Variation  der  will- 
kürlichen Konstanten.  (Lagrange,  Oeuvres, t. 4, p.  156.  Recherches 
sur  les  series  recurrentes.) 

Als  Beispiel  nehmen  wir  die  Gleichung 


Wir  setzen 
Es  ist 


wenn 


yx+^-^Vx+i  +  ^^x-^' 

i^  X  X  '  X 

•'a+l  X       '        X  ' 

AO(i)+6^+i-AO(2)=0 


ist. 
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Es  wird  y    der  Gleichung  genügen,  wenn  außerdem 

X  X 

ist.     Daraus  folgt 

X  5     '  *  30        \6/ 

Die  Summation  gibt 

Die  allgemeine  Lösung  unsrer  Gleichung  ist  also 

S8.  Man  kann  noch  in  einer  andern  Weise  von  der  Lösung 
der  homogenen  Gleichung  zur  vollständigen  (15)  übergehen.  Es 
genügt,  eine  partikuläre  Lösung  der  vollständigen  Gleichung  zu 
kennen  (Nr.  48).     Wenn  das  letzte  Glied 

ist,  wo  Vg  eine  ganze  Funktion  ist,  dann  ist  es  leicht,  eine  solche 
Lösung  zu  finden. 

Wir  setzen  „ 

wo  Ux  eine  ganze  Funktion  ist,  und  führen  diese  Größe  statt  y^  in 
die  vollständige  Gleichung  ein.  Vermöge  der  Identität  (4)  findet 
man  die  Gleichung 

a-fia)  .  u,+  a-+y-  f{a)  •  Aw,+  a-+2  ./""(a)  •  -^  +  .  • . 

. . .  +  a- + «  . /•  W  (a)  .  ^^^  =  a- .  i;„ 

die  zur  Bestimmung  von  Ux  dient. 

Der  Faktor  a*  kann  gehoben  werden  und  es  wird  diese 
Gleichung  die  Form  annehmen 

A'w  A"w 

(19)  f{a)  .  u,  +  af'(a)-Au,  +  ar'ia)  ■j7i+--  -i-a^f^^Ka)  ■  ^-^-f-  =  t;,. 

Wenn  a  keine  Wurzel  der  Gleichung 
ist,  dann  ist  Ux  von  demselben  Grade  wie  Vx- 


Ux  —   (^1 
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Es  sei 

x(x—l)(x  —  2)...(x  —  m-\-l)  x{x  —  l)...{x  —  m-\-2) 


1  •  2  •  3  . . .  m  "^  ^1       1  •  2  •  3  .  . .  (m  -  1) 


x(x—i)...{x  —  m-\-3)  x{x  —  l)    , 

"^  "2 1  .  2  .  .  .  (m  -  2) ' ^  ««-2      1.2      +  CCm-lOS  +  ««; 

^•*  ~  Po  1 .  2  •  3  . .  .  m  +  ^1       1.2.3...(m-l)"~ 

+  P2  1  .  2  .  .  .  (m  -  2) r-  •  •  •  +  Pm-2      ^     g H  Pm-X^  -r  Pm- 

Die  Zahlen  /3  sind  gegeben  und  die  a  werden  gesucht. 
Aus  (19)  findet  man  durch  den  Vergleich  der  Koeffizienten 

aj(a)  +  a^afia)  =  ß^, 
Da  /■(«)  nicht  Null  ist,  so  bestimmen  sich  sukzessive  die  Zahlen 

«Q,    ^IJ    ^2?    ^3>   •  •  * 

Wenn  a   eine   j*!; fache  Wurzel   von   f{d)  =  0   ist,   so  nimmt   die 
Gleichung  (19)  die  Form  an 

a^f^Ha)  ■  ^^±-,  +  «^+^  •  P+^)(a)  •  ,.,...(4,)  +  •  •  ■ 

A"w^ 

Die    linke    Seite    dieser    Gleichung    ist    eine    ganze    Funktion 
m^^^  Grades,  wenn  Ux  vom  (Ji  +  w)*®'^  Grade  ist. 
Man  kann  setzen 


Ux=  a, 


x{x—l)...{x  —  'k  —  m-\-l)  x(x  —  1) . . .  {x  —  k  —  m  -{■  2)    . 

'        1  •  2  •  3  . .  .  (Ä  +  m)  ^  "1       l-2-3...(fc  +  TO-l)         *" 


,  a;(a;— 1) . . .  (a;  — Ä:)  a;(a;  — 1) . . .  te  — ä;  + 1) 

"^  "»"-^   1.2.3...(Ä;+1)    "^  ^"^  1  ■  2  ■  3  . .  .  fc  ' 

und   zur  Bestimmung   der  Koeffizienten  findet   man   die  Gleichungen 

0      l-2...Jc~  P^' 

"1^   1.2...ä;+''o^  1.2.3...(ä;  +  1)~^1' 

«  gl.    n^m      ,   ^  .,,+1        /•(^+i)(a)  ,  +  ,        /■(^  +  2)(a) 

*^2«  1.2.../c  +  "i^  1.2...(Ä  +  l)  +  '^o^  1-2. ..(Ä  +  2)        '^^^ 
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Diese  Gleichungen   sind  lösbar,   weil  a  und  p^'^{a)  nicht  gleich 
Null  sind. 

Z.  B.  wollen  wir  eine  partikuläre  Lösung  der  Gleichung 

suchen.     Hier   ist  a  =  1    und   dieser  Wert   ist   eine  einfache  Wurzel 

der  Gleichung 

^2-7^  +  6  =  0. 

Deswegen  gibt  es  eine  partikuläre  Lösung  von  der  Form 

Li  die  gegebene  Gleichung  setzen  wir  ein 

Yx==bx  -\-  cx^ 
und  wir  finden 

^  =  5Ö'      ^  =  -10- 

Nach  dem   bekannten  Satze   (Nr.  48)   ist   die  allgemeine  Lösung 
dieser  Gleichung 


,,^C,  +  G,-&'+~x-i.x\ 


Z       _   1^ 
6Ö  ^        10 ' 


Die  Gleichung 
hat  eine  partikuläre  Lösung  von  der  Form 

Aus  der  Gleichung 

«0—  5a(,+  6a(j=  1 

folgt,  daß  «0  =  ^  is*-      Deswegen   ist   die    allgemeine  Lösung  unsrer 
Gleichung 

Wir  wollen  jetzt  einige  Anwendungen  der  Differenzengleichungen 
höherer  Ordnung  zeigen. 

69.  Rekurrente  Reihen.     Dies  sind  solche  Eeihen 

Vat  Vij  Vif  •  •  -y  Vx'  Vx-^-v  '  '  '> 
für  die  je  n  aufeinander  folgende  Glieder  einer  Relation 

Vx^n^  PxVx+n-l-^  P^yx+n-2+  '     '  +  PnVx^  ^ 

mit  konstanten  Koeffizienten  genügen. 

Seliwanoff,  Differenzenreolmung.  6* 
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Die  Auflösung  dieser  Differenzenglei ctung  dient  dazu,    um  das 
allgemeine  Glied  y^  als  Funktion  von  x  darzustellen. 
Als  Beispiel  wollen  wir  die  Reihe  betrachten 

0,  1,  1,  2,  3,  5,  8,..., 

deren  jedes  Grlied  die  Summe  der  beiden  vorhergebenden  Glieder  ist.  Je 
drei  nacheinander  folgende  Glieder  sind  durch  die  Relation  verbunden 

Die  Anfangswerte  sind 

Vo  =0,     2/i  =  1. 
Die  Wurzeln  der  quadratischen  Gleichung 

^2-^-1=0 

sind  ^  ^ 

Deswegen   ist   die    allgemeine   Lösung    der   Differenzengleichung 

Wenn    C^  und  C^   durch   die  Anfangswerte  «/^  und  y^  bestimmt 
werden,  findet  man,  daß  das  allgemeine  Glied  unsrer  Reihe 

ist. 

Die  Potenzreihe 

F(t)  =  2/0  +  yJ  +  y,t'  +  y,t'+-+  yj^  +  •  •  •. 

deren  Koeffizienten  eine  rekurrente  Reihe  bilden,  ist  gleich  einer 
rationalen  Funktion  von  t.  Die  Relation  zwischen  den  n  Koeffizienten 
der  Reihe  sei 

y,+„ + Piy.+n-i  +  p,y.+„-,  +  •  •  • + i>.?/.  =  o. 

Nach  Multiplikation  von  F(t)  mit 
findet  man 

yo  +  (yi-^Piyo)i + (2/2  -^Ptyi+p^yo^ + (2/3  ^Piy^  -^p^yx  +Psyo)i^ + 

Die  Koeffizienten   der   folgenden   Potenzen   sind   die  Werte    der 
Funktion 

y,  +  n  +  P,y,+n-l+P2yx+n^S-^--+Pnyx 

für  x  =  0,  1,  2, .  .  .    Deswegen  sind  alle  diese  Koeffizienten  gleich  Null. 
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Umgekehrt,  die  rationale  Funktion 

läßt  sich  in  die  Reihe 

entwickeln,  für  die  je  n  aufeinander  folgende  Koeffizienten  durch  die 
Relation 

verbunden  sind.     Die  n  ersten  Koeffizienten  bestimmen  sich   aus  den 
Gleichungen 

y2+Piyi+P2yo  =  <h> 


y„-l  +  P^y„_2  +  P2yn-B  +  '-'+Pn-iyO=^n-.l' 

Also  können  die  Difi'erenzengleichungen  dazu  dienen,  um  die 
Koeffizienten  der  Entwicklungen  rationaler  Funktionen  in  Potenz- 
reihen zu  bestimmen. 


60.    Tschebyscheffscher  Beweis  eines  Satzes  von  Lam6. 

In  seinen  Vorlesungen  zeigte  Tschebyscheff  die  folgende  An- 
wendung der  Differenzengleichungen. 

Man  soll  bestimmen,  wie  groß  die  Anzahl  der  sukzessiven 
Divisionen  sein  kann,  um  den  größten  gemeinschaftlichen  Teiler 
zweier  Zahlen  a  und  h  zu  finden. 

Es  sei 

«  =  &9'i  +  **!>     ^  =  ^iQ2+^i>  •  ",     ^n-2  =  r„-iq^  +  r„,     r,,  =  0. 

Die  Anzahl  der  Divisionen  ist  hier  n. 

Weil  alle  Quotienten  größer  oder  gleich  Eins  sind,  so  ist 

a^l  +  r^,     h^r^  +  r^,     r^^  r^  +  r^,  .  .  . , 

. .  .  »-„-s  ^  r„_2  +  r„_i,     r„_2^r„_i,    r„_i  ^  1. 

Daraus  folgen  die  Beziehungen 

r„=0      r«_i^l,     rn_2^1,     r„_8^2, 
rn-4  ^  3,    r„_5  ^  5,    r„_6  ^8,  ... 
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Wenu  wir  die  Zahlen 

2/0  =  0,    y,=  l,    2/2  =  1,    2/^  =  2,...,    y,=  y,_,  +  y,_„  ... 

einfüliren;  so  ist 

^n  =  Vq,     r^-i  ^  «/i,     r«_2  ^  2/27  •  •  •  ;     ^  ^  «/„• 
Wenn  für  eine. ganze  positjye  Zahl  m  die  Ungleichheit 

besteht,    dann    ist   notwendig   y^  >  y^    und    m>  n.      Die    gesuchte 
Anzahl  der  Divisionen  ist  also  kleiner  als  eine  solche  Zahl  m,  für  die 

y   >b 

ist.     Den  Ausdruck   für  y^  kennen    wir    schon   aus   Nr.  59.      Es    ist 
deswegen  die  Ungleichheit 


aufzulösen. 
Da 


/5 


y5  =  2,2...;     1+J^==1,6...;     ilJ^  =  -  0,  6  . 
so  können  wir  y    annähernd  in  der  Weise  ausdrücken 


«=Ä(^T±«- 


WO  0  zwischen  0  und  1  ist. 
Wenn 

ist,  dann  ist 

2/^  ^  &  +  (1  ±  0)  und  folglich  y^  >  K 
Es  genügt 

log[2,2.(t+l)] 

=  log  1,6 

ZU  nehmen.     Wenn  die  Zahl  h  aus    Zj  Ziffern  besteht,    dann    ist   der 
gemeine  Logarithmus  von  2,  2-(&  +  1)  ungefähr  gleich  h,  log  1,  6 

ißt  beinahe  —  und 

Daraus  folgt  der  Satz  von  G.  Lame:  Die  Anzahl  der  Divi- 
sionen, die  auszuführen  ist,  um  den  größten  gemeinschaft- 
lichen Teiler  der  Zahlen  a  und  &  (a  >  6)  zu  finden,  ist  kleiner 
als  die  fünffache  Anzahl  der  Ziffern  der  Zahl  &. 

Lame  hat  diesen  Satz  auf  elementarem  Wege  bewiesen.  (Paris. 
Comptes  rendus,  19.  [1844],  p.  867.) 
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Encyklopädie  der  Elementar^Mathematik. 

Ein  Handbucb  fllr  Lehrer  und  Studierende. 

Von 

Heinrich  Weber,       und         Josef  Wellstein, 

Professor  in  StraBburg  Professor  in  GieBen. 

In  3  Bänden.  Band  I  (XIV  u.  446  S.)  gr.  8.   1903.  In  Leinw.  geb.  n.  ^  8.— 

Das  Werk,  dessen  erster  Band  soeben  erschienen  ist,  richtet  sich  in  erster 
Linie  an  die  Lehrer,  die  darin  Anregung  finden  sollen,  ihren  Unterrichtsstoff 
auszuwählen  und  namentlich  in  den  höheren  Klassen  zu  vertiefen,  sodann  aber 
auch  an  Studierende,  die  eine  Anlehnung  an  die  Elemente  und  Auffrischung  und 
Ergänzung  früher  erworbener  Kenntnisse  suchen. 

Durch  das  Zusammenwirken  mehrerer  Gelehrter  hoffen  die  Herausgeber,  die 
möglichste  Vollständigkeit  zu  erreichen.  Der  erste  Band  umfaßt  den  algebraisch- 
analytischen Teil.  Der  zweite  Band,  der  Ostern  1904  erscheint,  wird  die  Geometrie 
nach  ihren  verschiedenen  Seiten  behandeln.  Ein  dritter  Teil,  der  im  Herbst  1904 
zur  Ausgabe  gelangen  soll,  wird  die  Anwendungen  bringen,  deren  Stoff  aus  der 
darstellenden  Geometrie,  der  Mechanik,  Physik  und  Wahrscheinlichkeitsrechnung 
entnommen  ist.  Die  Vorarbeiten  sind  so  weit  gediehen,  daß  die  Vollendung  des 
ganzen  Werkes  im  nächsten  Jahre  erwartet  werden  darf. 


Yerlag  Yon  B.  Gl.  Teubner  in  Leipzig. 

Markoff,  A.  A.,  o.  Prof.  an  der  Kaiserl.  Universität  zu  St.  Petersburg 
usw.,  Differenzenrechnung.  Autorisierte  deutsche  Übersetzung  von 
Th.  Friesendorff  und  E.  Prümm.  Mit  einem  Vorwort  von  R.  Mehmke, 
o.  Prof.  an  der  Kgl.  Techn.  Hochschule  zu  Stuttgart.  [VI  u.  194  S.] 
gr.  8.     1896.     geh.  n.  M  7.— 

Das  Mark  off  sehe  Werk,  das  in  Rußland  sehr  verbreitet  ist  und  hier  in  deutscher  Über- 
setzung vorgelegt  wird,  kann  einerseits  als  ein  Lehrbuch  dienen,  anderseits  durch  seine  Vielseitig- 
keit ein  allgemeines  Interesse  erwecken,  besonders  da  die  Differenzenrechnung  in  letzter  Zeit  in 
Deutschland  wenig  bearbeitet  worden  ist.    Das  Buch  enthält  zwei  Teile. 

Der  erste  Teil  (der  im  Jahre  1889  erschienen  ist)  behandelt  die  Interpolation  und  ihre 
Anwendungen.  Nachdem  der  Verfasser  in  Kap.  I  die  Hauptinterpolationsformeln  abgeleitet  hat,  mit 
Berücksichtigung  des  Eestgliedes,  geht  er  in  Kap.  II  zu  den  allgemeinen  Begriffen  der  Differenzen 
verschiedener  Ordnungen  liber.  In  Kap.  IH  wird  der  Zusammenhang  zwischen  Differenzen  und 
Differentialquotienten  auseinandergesetzt.  In  Kap.  FV  werden  Herstellung  und  Benutzung  mathe- 
matischer Tabellen  besprochen  und  an  einigen  Beispielen  erläutert.  Die  Kap.  V,  VI  und  VII  ent- 
halten Anwendungen  der  Interpolation  auf  die  Berechnung  bestimmter  Integrale,  einige  Eigen- 
schaften der  Degendreschen  Funktionen,  Bntwickelungen  einiger  Integrale  in  Kettenbrüohe  u.  dgl. 

Der   zweite    Teil   (der   im  Jahre  1891   erschienen  ist)   behandelt   die  Summation  und  die 

a-\-hn 

Theorie   der  Differenzengleiohungen.    In  Kap.  I  wird  der  Begriff  der  Summe  2_\  /(^)   gegeben  und 

a 
die  Summation  einiger  Ausdrücke  angeführt.    In  Kap.  II  und  III  wird  die  Eulersche  Formel 
h  h 

2/(^)  =  t/^(^)  ^^  -  Y  i  ^(*)  -/(«)}  +  ^  { /  (6)  -f{a)\.... 
a  a 

abgeleitet  und  interpretiert  und  ihre  Anwendung  bei  der  Berechnung  einiger  Integrale  und  Summen 
auseinandergesetzt,  sowie  aus  ihr  die  Stirlingsche  Formel: 

lg  1.2.3...  a:  =  -i-lg  27r -f  ^x-f -^  lg  a:  -  X  +  A -^  +  ^.  ~  +  ••■ 

abgeleitet.  In  Kap.  IV  werden  Begriffe  und  Theorie  der  Differenzengleichungen  gegeben.  In  Kap.  V 
und  VI  werden  lineare  Differenzengleichungen  erster  Ordnung  besprochen.  Endlich  Kap.  VII  enthält 
die  Anwendung  der  Doppelsummen  auf  die  Transformation  der  Reihen,  an  einigen  Beispielen  erläutert. 
Es  ist  zu  bemerken,  daß  der  Verfasser  ganz  besonders  die  Abschätzung  der  Fehler  der 
Näherungsformeln  hervorgehoben  hat,  die  z.  B.  im  Buche  von  Schlömilch  „Theorie  der  Differenzen 
und  Summen"  nicht  so  ausführlich  behandelt  wird.  Außerdem  folgt  nach  jedem  Kapitel  ein  aus- 
führliches Literaturverzeichnis,  das  natürlich  eine  große  Erleichterung  für  den  Leser  bietet. 

Klein,  Geh.  Reg.- Rat  Dr.  Felix,  Professor  an  der  Universität  Göttingen, 
Anwendung  der  Differential-  und  Integralrechnung  auf  Geo- 
metrie, eine  Revision  der  Prinzipien.  Vorlesung,  gehalten  während 
des  Sommersemesters  1901.  Ausgearbeitet  von  Conrad  Müller.  (Der 
autographierten  Vorlesungen  VIT.  Heft.)    [VIII  u.  468  S.]   4.    1902.   geh. 

n.  .//  10.— 

Der  Verfasser  tut  mit  der  Herausgabe  einen  neuen  Schritt  in  der  Richtung,  der  seine  Be- 
mühungen in  den  letztvergangenen  Jahren  vorwiegend  gegolten  haben:  die  mathematische  Wissen- 
schaft als  ein  zusammengehöriges  Ganzes  nach  allen  Seiten  wieder  so  zur  Geltung  zu  bringen,  wie 
es  früher,  als  es  noch  keinen  Spezialismus  gab,  als  selbstverständlich  galt,  insbesondere  aber  zu 
erreichen,  daß  zwischen  den  Vertretern  der  abstrakten  uud  der  angewandten  Mathematik  wieder 
mehr  gegenseitiges  Verständnis  Platz  greife,  als  augenblicklich  besteht.  Er  hat  schon  wiederholt 
darauf  hingewiesen,  daß  in  letzterer  Hinsicht  eine  klare  Erfassung  des  Unterschiedes  und  der  gegen- 
seitigen Beziehungen  derjenigen  beiden  Teile  der  Mathematik  von  zentraler  Bedeutung  sein  dürfte, 
die  im  folgenden  als  Präzisionsmathematik  und  Approximationsmathematik  bezeichnet 
werden.  Ähnliche  Gedanken  sind  in  den  letzten  Jahren  von  den  Herren  Burkhardt  und  Heun 
geäußert  worden ;  um  so  mehr  darf  vielleicht  eine  ausführliche  Auseinanderlegung  der  in  Betracht 
kommenden  Gegensätze,  tvie  sie  die  folgende  Vorlesung  für  das  Gebiet  der  Geometrie  gibt,  nun- 
mehr auf  Interesse  und  Varständnis  in  weiteren  Kreisen  rechnen. 

Mit  dieser  Vorlesung  vervollständigt  der  Verfasser  die  Auseinandersetzung,  welche  er  ver- 
schiedentlich über  die  Methoden  des  mathematischen  Unterrichtes ,  insbesondere  des  mathematischen 
Hochschulunterrichtes,  gegeben  hat.  Seine  Ansicht  ist  nach  wie  vor,  daß  der  Unterricht  der  An- 
fänger, sowie  derjenigen  Zuhörer,  welche  die  Mathematik  nur  als  ein  Hilfsmittel  für  anderweitige 
Studien  gebrauchen  wollen,  von  den  anschauungsmäßigen  Momenten  einen  naiven  Gebrauch  machen 
soll;  die  Überzeugung,  daß  dies  aus  pädagogischeu  Gründen,  mit  Rücksicht  auf  die  Veranlagung 
der  Mehrzahl  der  Studierenden,  notwendig  sei,  kommt  in  den  letzten  Jahren  sichtlich  immer  mehr, 
insbesondere  auch  im  Auslande,  zur  Geltung.  Aber  nicht  minder  ist  seine  Überzeugung,  daß  ein 
solcher  Unterricht,  entsprechend  der  heutigen  Entwickelung  der  Wissenschaft,  bei  der  Ausbildung 
der  Faohmathematiker  auf  der  Oberstufo  nicht  genügt,  daß  hier  vielmehr  neben  den  Tatsachen 
der  Anschauung  die  zentrale  Bedeutung  des  modernen  Zahlbegriffes  und  der  weitgehenden ,  •  mit 
ihm  zusammenhängenden  Entwickelungen  hervortreten  muß.  Er  vermißt  in  den  Lehrbüchern  und 
Vorlesungen,  welche  er  kennt,  die  Überleitung  von  der  einen  Auffassung  zur  anderen.  Hier  will 
sich  die  folgende  Darstellung  als  eine  Ergänzung  einschieben ;  ihr  bestes  Ziel  wird  erreicht  sein, 
wenn  sie  sich  eines  Tages  als  überflüssig  erweisen  sollte,  weil  die  Überlegungen,  die  sie  bringt,  zu 
selbstverständlichen  Bestandteilen  jedes  höheren  mathematischen  Unterrichtes  geworden  sind. 


(^  Selivanov,  Dmitrii  Fedorovich 
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